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Abstract. We analyze the causal action principle for a system of relativistic fer- 
mions composed of massive Dirac particles and neutrinos. In the continuum limit, 
we obtain an effective interaction described by a left-handed, massive SU(2)-gauge 
field and a gravitational field. The off-diagonal gauge potentials involve a unitary 
mixing matrix, which is similar to the Maki-Nakagawa-Sakata matrix in the standard 
model. 
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1. Introduction 

In [3j it was proposed to formulate physics based on a new action principle in space- 
time. In the paper [7], this action principle was worked out in detail in the so-called 
continuum limit for a simple model involving several generations of massive Dirac 
particles. In the present article, we extend this analysis to a model which includes 
neutrinos. In the continuum limit, we obtain an effective interaction described by a 
left-handed massive SU(2) gauge field and a gravitational field. 

More specifically, we again consider the causal action principle introduced in [4]. 
Thus we define the causal Lagrangian by 

C\A xy \ = |^-xy I ~~ g l^icyl 2 , (1.1) 

where A xy = P(x, y) P(y, x) denotes the closed chain corresponding to the fermionic 
projector P(x,y), and \A\ = Ya=x M ^ s ^ ne spectral weight (where Aj are the eigen- 
values of A counted with algebraic multiplicities). We introduce the action S and the 
constraint 7~ by 

S[P] = [[ £[A xy ]d 4 xd 4 y, T[P] = 11 \A xy \ 2 d 4 x d 4 y , (1.2) 
JJmxm JJmxm 

where (M, (., .)) denotes Minkowski space. The causal action principle is to 

minimize S for fixed T ■ (1-3) 

This action principle is given a rigorous meaning in Section 2]. Every minimizer is 
a critical point of the so-called auxiliary action 

S„[P] = [[ C^A xy ] d 4 x d 4 y , Cfj[A xy \ = \A 2 \ - y, \A xy \ 2 , (1.4) 

J JMxM 

which involves a Lagrange multiplier fj, € R. 

We model the configuration of the fermions by a system consisting of a doublet of 
two sectors, each composed of three generations. Thus we describe the vacuum by the 
fermionic projector 

P(x,y) = P N (x,y)®P c (x,y), (1.5) 
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where the charged sector P c is formed exactly as the fermionic projector in [7] as a 
sum of Dirac seas, i.e. 

3 

P C (x,y) =J2 P m (x,y), (1-6) 

/3=1 

where mp are the masses of the fermions and P m is the distribution 

Pm(x, y) = I (# + m) 5(k 2 - m 2 ) 0(-k°) e~ lk ^ . (1.7) 

For the neutrino sector P N we consider two different ansatze. The first ansatz of chiral 
neutrinos is to take a sum of left-handed, massless Dirac seas, 

3 

P N (x,y) = Y J XLP (x,y). (1.8) 

/3=1 

The configuration of Dirac seas (|1.5p . fjl .6|) and (jl.8p models precisely the leptons in 
the standard model. It was considered earlier in [H §5.1]. The chiral ansatz (I1.8P 
has the shortcoming that the neutrinos are necessarily massless, in contradiction to 
experimental observations. In order to describe massive neutrinos, we proceed as 
follows. As the mass mixes the left- and right-handed spinor components in the Dirac 
equation, for massive Dirac particles it is impossible to to restrict attention to one 
chirality. This leads us to the ansatz of massive neutrinos 

3 

P N (x,y) = ^Prr l ,(x,y). (1.9) 

/3=1 

Here the neutrino masses rhp > will in general be different from the masses mp 
in the charged sector. Except for the different masses, the ansatze (|1.6p and (|1.9p 
are exactly the same. In particular, it might seem surprising that (|1.9p does not 
distinguish the left- or right-handed component, in contrast to the observation that 
neutrinos are always left-handed. In order to obtain consistency with experiments, 
if working with (|1.9p we need to make sure that the interaction distinguishes one 
chirality. For example, if we described massive neutrinos by f 1 1 . 9 [) and found that the 
neutrinos only couple to left-handed gauge fields, then the right-handed neutrinos, 
although being present in (|1.9p . would not be observable. With this in mind, working 
with (ll.9p seems a possible approach, provided that we find a way to break the chiral 
symmetry in the interaction. It is a major goal of this paper to work out how this can 
be accomplished. 

Working out the continuum limit for the above systems gives the following results. 
First, we rule out the chiral ansatz (jl.8D by showing that it does not admit a global 
minimizer of the causal action principle. Thus in the fermionic projector approach, we 
must necessarily work with the massive ansatz (]1.9p . We find that at least one of the 
neutrino masses fhp must be strictly positive. In order to break the chiral symmetry, 
we introduce additional right-handed states into the neutrino sector. It is a delicate 
question how this should be done. We discuss different approaches, in particular the 
so-called shear states and general surface states. The conclusion is that if the right- 
handed states and the regularization are introduced suitably, then the continuum limit 
is well-defined. Moreover, the structure of the effective interaction in the continuum 
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limit is described as follows. The fermions satisfy the Dirac equation coupled to a 
left-handed SU(2)-gauge potential Al, 

f = 0, (1.10) 

where we used a block matrix notation (where the matrix entries are 3 x 3-matrices). 
Here mY is a diagonal matrix composed of the fermion masses, 

mY = diag(mi, 7712, m^, mi, 1112,1713) , 

and C/mns is a unitary 3 x 3-matrix. In analogy to the standard model, we refer 
to Uuns as the Maki-Nakagawa-Sakata (MNS) matrix. The gauge potentials Al satisfy 
a classical Yang-Mills equation coupled to the fermions. More precisely, writing the 
isospin dependence of the gauge potentials according to Al = X^a=i ^L^ a m terms of 
Pauli matrices, we obtain the field equations 

d\{ATf - U(At) k - Ml {A%) k = c a *{ XLl k * a )* , (1.11) 

valid for a = 1,2, 3. Here M a are the bosonic masses and c a the corresponding coupling 
constants. The masses and coupling constants of the two off-diagonal components are 
equal, i.e. M\ = M2 and c\ = C2, but they may be different from the mass and coupling 
constant of the diagonal component a = 3. 

Moreover, our model involves a gravitational field described by the Einstein equa- 
tions 

Rjh~ 2 R #/fc + ^9jk = KT jk , (1.12) 

where Rjk denotes the Ricci tensor, R is scalar curvature, and Tjk is the energy- 
momentum tensor of the Dirac field. Moreover, n and A denote the gravitational and 
the cosmological constants, respectively. We find that the gravitational constant scales 
like k ~ 5~ 2 , where 5 is the length scale on which the shear and general surface states 
become relevant. The dynamics in the continuum limit is described by the coupled 
Dirac- Yang/Mills-Einstein equations (jl.lOp . (jl.lip and (|1.12p . These equations are 
of variational form, meaning that they can be recovered as Euler-Lagrange equations 
corresponding to an "effective action." The effective continuum theory is manifestly 
covariant under general coordinate transformations. 

For ease in notation, the field equations (and similarly the Einstein equa- 

tions (I1.12p ) were written only for one fermionic wave function But clearly, the 
equations hold similarly for many- fermion systems (see Theorem 17. ip . In this context, 
it is worth noting that, although the states of the Dirac sea are explicitly taken into 
account in our analysis, they do not enter the Einstein equations. Thus the naive 
"infinite negative energy density" of the sea drops out of the field equations, making 
it unnecessary to subtract any counter terms. 

Similar as explained in [7] for an axial field, we again obtain corrections to the field 
equations which are nonlocal and violate causality in the sense that the future may 
influence the past. Moreover, for a given regularization one can compute the coupling 
constant, the bosonic mass, and the gravitational constant. 

We note that in this paper, we restrict attention to explaining our computations 
and results; for all conceptual issues and more references we refer to [7] and the survey 
article [9]. 



/ 4 1 (Afu Mm y\ 

' + Xfi 21 11 — m ^ 

\4l Umns ~4l J 



THE CONTINUUM LIMIT OF A FERMION SYSTEM INVOLVING NEUTRINOS 



5 



2. Regularizing the Neutrino Sector 

In this section, we explain how the neutrino sector is to be regularized. We begin 
in E j2.ll by reviewing the regularization method used in [3] . Then we give an argument 
why this method is not sufficient for our purposes (see ^2.2[) . This leads us to extending 
our methods (see ^2.3j) . and we will explain why these methods only work for the 
ansatz of massive neutrinos (see £ j2.4p . In £ 12.51 we introduce the resulting general 
regularization scheme for the vacuum neutrino sector. In §2.61 we explain how to 
introduce an interaction, relying for the more technical aspects on Appendix[A] Finally, 
in Q2.7\ we introduce a modification of the formalism of the continuum limit which 
makes some computations more transparent. 

2.1. A Naive Regularization of the Neutrino Sector. As in [7, Section 3] we 
denote the regularized fermionic projector of the vacuum by P £ , where the parameter e 
is the length scale of the regularization. This regularization length can be thought of as 
the Planck length, but it could be even smaller. Here we shall always assume that P £ 
is homogeneous, meaning that it depends only on the difference vector £ := y — x. This 
is a natural physical assumption as the vacuum state should not distinguish a specific 
point in space-time. The simplest regularization method for the vacuum neutrino 
sector is to replace the above distribution P N (x, y) (see fjl .8[) ) by a function which 
is again left-handed, 

P?{x,y) = X L9M)l j ■ (2-1) 
Such a regularization, in what follows referred to as a naive regularization, was used 
in [4] (see [H eq. (5.3.1)]). It has the effect that the corresponding closed chain vanishes 
due to so-called chiral cancellations (see [H eq. 5.3.2]), 

A xy-=P^{x,y)Pe{y,x) = XL^{x,y)xL^{v,x) = xl xr <${ x , y) <$(y, x ) = 0. 

Regularizing the charged sector as explained in [4, Chapter 4] or [7j, the closed chain 
of the regularized fermionic projector P £ of the whole system is of the form 

A xy = P £ (x, y) P £ (y, x) = 0®Ag. 

Hence the closed chain has the eigenvalue zero with multiplicity four as well as the 
non-trivial eigenvalues A + and A_, both with multiplicity two (see [H §5.3]). Let us 
recall from [H Chapter 5] how by a specific choice of the Lagrange multiplier \i we 
can arrange that the EL equations are satisfied: The operator Q corresponding to the 
action f j 1 . 4 [) is computed by (see [H §3.5] or [7J Section 6]) 

Q(x, y) = {-An) [(1 - 4/x)] T ' F < f )• 
In order for the operator Q to vanish on the charged sector, we must choose 

M = \ ■ (2-2) 

Then 

Q( x , y) = -J2 ^ y) © °> 

s=± 

and multiplying by P(y, z), we again get chiral cancellations to obtain 
Q(x, y)P(y, z) = -Y, XL 4{x, z) X l j(z, y) = . 
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Similarly, the pointwise product P(x,y) Q(y, z) also vanishes, showing that the EL 
equations [P, Q] are indeed satisfied in the vacuum. 

Before going on, we remark for clarity that in [3], the chiral regularization ansatz (|2.ip 
was overridden on the large scale in order to arrange a suitable normalization of the 
chiral fermionic states (see [U Appendix C]). More precisely, P^ was constructed by 
projecting out half of the states of a Dirac sea of mass m. The formula (|2.ip was 
recovered in the limit m \ 0. In this so-called singular mass limit, the normaliza- 
tion integrals did not converge, making it possible to arrange a proper normalization, 
although for the limit (|2,ip the normalization integral would vanish due to chiral can- 
cellations. However, in [JJ §C.l] it was explained that the formalism of the continuum 
limit is well-behaved in the singular mass limit, thus justifying why we were allowed 
to describe the regularized chiral Dirac seas by (|2.ip . 

2.2. Instability of the Naively Regularized Neutrino Sector. We now give an 
argument which shows that if the neutrino sector is regularized in the neutrino sector 
according to (12. ip . the system (II. 5p cannot be an absolute minimum of the causal 
action principle (11.3p . Suppose conversely that a fermionic projector P £ , which in the 
neutrino sector is regularized according to (|2.ip . is an absolute minimum of the action 
principle (|1.3p . Then any variation of the fermionic projector can only increase the 
action. Evaluating this condition for specific variations leads to the notion of state 
stability, which we now recall (for details see [U §5.6] or |12|). This notion makes it 
necessary to assume that our regularization is macroscopic away from the light cone, 
meaning that the difference P £ (x, y) — P(x, y) should be small pointwise except if the 
vector y — x is close to the light cone (see [H §5.6]). This condition seems to be fulfilled 
for any reasonable regularization, and thus we shall always assume it from now on. 
Suppose that the state ^ is occupied by a particle (i.e. that ^ lies in the image of 
the operator P £ ), whereas the state $ is not occupied. We assume that ^/ and <3? are 
suitably normalized and negative definite with respect to the indefinite inner product 

<^\(f» = / il)[x)<j>(x) d A x . (2.3) 
Jm 

Then the ansatz 

8P £ {x, y) = ^(x)^(y) - (2.4) 

describes an admissible perturbation of P £ . Since the number of occupied states is 
very large, 5P £ is a very small perturbation (which even becomes infinitesimally small 
in the infinite volume limit). Thus we may consider 5P as a first order variation and 
treat the constraint in (jl.3p with a Lagrange multiplier. We point out that the set 
of possible variations 5P £ does not form a vector space, because it is restricted by 
additional conditions. This is seen most easily from the fact that —5P £ is not an 
admissible variation, as it does not preserve the rank of P £ . The fact that possible 
variations 5P £ are restricted has the consequence that we merely get the variational 
inequality 

5 M [P £ + SP E ] > <S M [P e ] , (2.5) 

valid for all admissible variations of the form (12. 4D . 

Next, we consider variations which are homogeneous, meaning that Vl/ and <£ are 
plane waves of momenta k respectively q, 

y(x) = * e~ ikx , $(x) = $ e~ iqx . (2.6) 
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Then both P £ and the variation 5P depend only on the difference vector £ = y — x. 
Thus after carrying out one integral in (jl.4p . we obtain a constant, so that the second 
integral diverges. Thinking of the infinite volume limit of a system in finite 4-volume, 
we can remove this divergence simply by omitting the second integral. Then (|2.5p 
simplifies to the state stability condition 

[ 5^[A(O]d 4 Z>0. (2.7) 

JM 

In order to analyze state stability for our system (jl.5p , we first choose the Lagrange 
multiplier according to (|2.2|) . Moreover, we assume that ^ is a state of the charged 
sector, whereas <3? is in the neutrino sector, 

* = oe$ c , |> = d Ar ©0. (2.8) 

Since ^ should be an occupied state, it must clearly be a solution of one of the Dirac 
equations — m a )^ = with a G {1, 2, 3}. The state 3>, on the other hand, should 
be unoccupied; we assume for simplicity that its momentum q is outside the support 
ofPf, 

q supp g (2.9) 

(where g is the Fourier transform of the vector field g in (|2.ip ). Thus our variation 
removes a state from a Dirac sea in the charged sector and occupies instead an unoc- 
cupied state in the neutrino sector with arbitrary momentum q (in particular, $ does 
not need to satisfy any Dirac equation). Let us compute the corresponding variation 
of the Lagrangian. First, using that the spectral weight is additive on direct sums, we 
find that 

6Ci = 5{\A 2 \-\\A\ 2 )=5\A 2 \-\\A\ 6\A\ 

= 5\(A C ) 2 \+5\(A N ) 2 \-±(\A C \ + \A N \) (^l + ^l) • 

This formula simplifies if we use that vanishes due to chiral cancellations. More- 
over, the first order variation of (A ) 2 vanishes because 

5 ((A N ) 2 ) = (5A N )A N + A N (SA N ) = . 

Finally, 8\A N \ = \(A N + 5A N )\ -\A N \ = \SA N \. This gives 

8Ci = 5(\(A C ) 2 \ - \ \A C \ 2 ) - I \A C \ \SA N \ . (2.10) 
4 V 4 / 2 

Note that only affects the first term, whereas $ influences only the second term. In 
the first term the neutrino sector does not appear, and thus the state stability analysis 
for one sector as carried out in [U §5.6] and |12j applies. From this analysis, we know 
that the charged sector should be regularized in compliance with the condition of a 
distributional A^P-product (see also [6]). Then the first term in (|2.10p leads to a finite 
variation of our action. The point is that the second term in (I2.10p is negative. In the 
next lemma we show that it is even unbounded below, proving that our system indeed 
violates the state stability condition (|2.7p . 

Lemma 2.1. Suppose that P £ is a regularization of the distribution (jl.5p which is 
macroscopic away from the light cone and which in the neutrino sector is of the 
form (|2.ip . Then for any constant C > there is a properly normalized, negative 
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definite wave function $ satisfying (|2.6p . (|2.8p and (|2.9p smc/i i/iot i/ie corresponding 
variation of the fermionic projector 

5P e {x,y) = -<S>{x)My) (2.11) 

satisfies the inequality 

I \A C \ \5A N \ d A i>C. 
Jm 

Proof. For convenience, we occupy two fermionic states of the same momentum q such 
that 

5P E N (x,y) = (jf + ro) e-^-^ , (2.12) 

where p is a vector on the lower hyperboloid H m := {p \ p 2 = m 2 and p° < 0}, and m 
is a positive parameter which involves the normalization constant. For this simple 
ansatz one easily verifies that the image of 5P N is indeed two-dimensional and negative 
definite. By occupying the two states in two separate steps, one can decompose (I2.12P 
into two variations of the required form ()2.11|) . Therefore, it suffices to prove the 
lemma for the variation fl2. 12f) . 

Using (|2.ip and ([2.12p . the variation of A N is computed to be 

SA N = xl^(x, y)(j> + m) e iq t + XR ($ + m ) <j( y , x ) e** . 

To simplify the notation, we omit the arguments x and y and write g(£) = g(x,y). 
Then g is a complex vector field with g(£) = g(y,x). Using that our regularization 
is macroscopic away from the light cone, there clearly is a set !] C M of positive 
Lebesgue measure such that both the vector field g and the function l^^l are non-zero 
for all £ G f2. Then we can choose a past directed null vector n such that (n,g) is 
non-zero on a set 0,' C O, again of positive measure. We now consider a sequence of 
vectors p\ £ 7i m which converge to the ray M + n in the sense that there are coefficients q 
with 

pi — ci n — > and q — > oo . 

Then on il', the inner product (pi,g) diverges as / — > oo. A short computation shows 
that in this limit, the eigenvalues of the matrix 5Af also diverge. Computing these 
eigenvalues asymptotically, one finds that 

\5A?\>A\(p u g)\+0(P). 

Hence for large /, 

/ \A C \\8A?\> [ \A c \\( Pl ,g)\ ^ oo, 
Jm Jn' 

completing the proof. □ 

It is remarkable that the above argument applies independent of any regularization 
details. We learn that regularizing the neutrino sector by a left-handed function ()2.ip 
necessarily leads to an instability of the vacuum. The only way to avoid this instability 
is to consider more general regularizations where also involves a right-handed 
component. 
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Figure 1. Plots of qr exemplifying different regularization mecha- 
nisms in the neutrino sector. 



2.3. Regularizing the Vacuum Neutrino Sector — Introductory Discussion. 

We begin by explaining our regularization method for one massless left-handed Dirac 
sea, 

P{x,y) = xl Po(x,y) 

(several seas and massive neutrinos will be considered later in this section). Working 
with a left-handed Dirac sea is motivated by the fact that right-handed neutrinos have 
never been observed in nature. To be precise, this physical observation only tells us 
that there should be no right-handed neutrinos in the low-energy regime. However, 
on the regularization scale e^ 1 , which is at least as large as the Planck energy Ep 
and therefore clearly inaccessible to experiments, there might well be right-handed 
neutrinos. Thus it seems physically admissible to regularize P by 

P e (x,y) = XL<i L (x,y) + Xr$r{x,v) , ( 2 -!3) 

provided that the Fourier transform gn(k) vanishes if \k° \ + \k\ <C e _1 . 

In order to explain the effect of such a right-handed high-energy component, we begin 
with the simplest example where gn is supported on the lower mass cone, 

|#) = 8i 2 PM^ 2 ), (2.14) 

where u = k , and the non-negative function h is supported in the high-energy re- 
gion to ~ e^ 1 (see Figure Q] (A)). We compute the Fourier integrals by 



(27T) 



2^ e iK 



-2i% 
2^ 



2tt K ' 



f 

Jo 



(27T)' 



h(u) 5{k 



2\ e iki 



where we set t - 
simple formula 



du 
2^ 
doj 
2^ 



h(u) 



iut 



JjJlt 



p 2 dp S(oj 2 — p 2 ) 



p dp 8{<jj 2 — p 2 )(e 



dcos$e- iprcos ^ 



-ipr 



') 



|£| and chose polar coordinates (p 



\k\, </?). This gives the 



where h is the one-dimensional Fourier transform of h. Under the natural assumption 
that the derivatives of h scale in powers of e, the function h decays rapidly on the 
regularization scale. Then ^ R vanishes except if £ is close to the light cone, so that 
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the regularization is again macroscopic away from the light cone. But the contribu- 
tion (|2.14p does affect the singularities on the light cone, and it is thus of importance 
in the continuum limit. More specifically, on the upper light cone away from the 
origin t ~ r 3> e, we obtain the contribution 

^ r -17 7 J r +7 r2 (215) 
+ (rapid decay in r) , 

where we set 7 r = (£7)/?". This contribution is compatible with the formalism of the 
continuum limit, because it has a similar structure and the same scaling as correspond- 
ing contributions by a regularized Dirac sea (see [6] , where the same notation and sign 
conventions are used). 

Regularizing the neutrino sector of our fermionic projector (jl.5p using a right-handed 
high-energy component has the consequence that no chiral cancellations occur. Hence 
the EL equations become 



where i labels the eigenvalues of A xy . For these equations to be satisfied, we must 
choose 

» = (2-17) 

and furthermore we must impose that the eigenvalues of A xy all have the same absolute 
values in the sense that 

(M - l A jl)rrT F i p (x,v) = for a11 ho- 

I Ail 

In simple terms, the matrix A^ must have the same spectral properties as A . 

This consideration points to a shortcoming of the regularization (|2.14p . Namely, 
the expression (|2.15p does not involve a mass parameter, and thus the corresponding 
contribution to the closed chain A N cannot have the same spectral properties as A , 
which has a non-trivial mass expansion. A possible solution to this problem is to 
consider states on a more general hypersurface, as we now explain again in the example 
of a spherically symmetric regularization. We choose 



1R 



(k) = -4vr 2 ( 7 ° + 7 fc ) h(u) 8(\k\ - K(u)) , (2.18) 



where j k = kj/k, and h is chosen as in (|2.14p . We again assume that § is supported 
in the high-energy region, meaning that 

h(cj) = if \u\ < e' 1 . (2.19) 
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Setting K = -co, we get back to (|2.14p ; but now the function K gives a more general 
dispersion relation (see Figure Q] (B)). Carrying out the Fourier integrals, we obtain 



-4tt z 



d 4 k 

W)~ 4 



h(u))$(\k\-K(u))) e ifc « 



A, , 

— h(co) e iut 

-™ 27T 



p 2 d P 5(p- K(u)) / dcosde~ ipr 



• cos i? 



2^ 



h(u))e wt / pdp5(p-K(u)) (e 



-ipr 



tpr\ 







2tt 



(19R)(0 = "4vr 2 (17V) 



iuit I giKr 



-iKr 



d 4 k 



-7V 

T 

r 

Y 



—00 




(2vr)4 

° du hi \ 
n(u!) e 



h(w)5(k- K(u)) 



2tt 



iuit I ^J-Kr 



\k\ 
iKr 



duj 

27T 



iuit I ^iKr _|_ ^—iKr 



giA"r ^—iKr 



Evaluating as in (|2.15p on the upper light cone away from the origin, we conclude that 



-00 

+ (rapid decay in r) . 



i^(w) + \ ) e i < w * +Kr > 



(2.20) 



For ease in notation, from now on we will omit the rapidly decaying error term. Re- 
arranging the exponentials, we obtain 



IR 



(0 



i{ui+K)r 



2tt 



c"~ 1 • '■ h(ui) 



Now the mass expansion can be performed by expanding the factor exp(i(w + K)r), 



h{u) (u> + K) n [i 



n=0 

r° du 
2vr 



7 - 7 



K{u) + ^ ) e iaj(t " r) (2.21) 



7 r K7 + a; 7 r 
+ i h 

^1 ^ ^» 



Jui(t-r) 



We conclude that the general ansatz (|2. 18|) gives rise to a mass expansion which is 
similar to that for a massive Dirac sea (see [H Chaper 4]). By modifying the geometry 
of the hypersurface {\k\ = K(lv)}, we have a lot of freedom to modify the contributions 
to the mass expansion. We point out that, in contrast to the mass expansion for a 
massive Dirac sea, the mass expansion in (|2.2ip involves no logarithmic poles. This is 
because here we only consider high-energy states (|2.19p . whereas the logarithmic poles 
are a consequence of the low-frequency behavior of the massive Dirac seas (for details 
see the discussion of the logarithmic mass problem in [H §2.5 and §4.3]). 
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We now come to another regularization effect. The regularizations (|2.14p and ()2.18p 
considered so far have the property that $ R is a multiple of the matrix xl(7° + 7 fc )> 
as is indicated in Figure ([1]) (B) by the arrows (to avoid confusion with the signs, we 
note that on the lower mass shell, k 1 = W7 — kj = uj (7 + 7 fc )). Clearly, we could also 
have flipped the sign of 7 fc , i.e. instead of (|2.18p . 

§ R {k) = -47T 2 (7 - 7") h(u) ~ K(u)) (2.22) 

(see Figure Q] (C)). In order to explain the consequence of this sign change in the 
simplest possible case, we consider the two functions 

$ ± (k) = 87r 2 u J ( 1 ± 1 k )h(u J )5(k 2 ), 

whose Fourier transforms are given in analogy to (|2.15p on the upper light cone by 

When multiplying ^ + by itself, the identity (7 + 7 r ) 2 = gives rise to a cancellation. 
For example, in the expression 



\ Tr W + (^ + R)«) = 2Re( "' (t -; )M '- r)> - (2-24) 

the term ~ r~ 2 has dropped out. The situation is different if we multiply ^ + by <fi_. 
For example, in 



1 T (A (A J mt^K 2ilm(h'(t-r)h{t-r)) \h(t - r)\ 2 

1 Tr(^ + (0^_(0 ) = ~ 2 ~ 3 + ~ 4 (2-25) 

no cancellation occurs, so that the term ~ r -2 is present. From this consideration we 
learn that by flipping the sign of 7 r as in (|2.22p . we can generate terms in the closed 
chain which have a different scaling behavior in the radius. 

In order to clarify the last construction, it is helpful to describe the situation in terms 
of the general notions introduced in [4, §4.4]. The fact that the leading term in ([2.15P 
is proportional to (7 — r y r ) can be expressed by saying that the vector component 
is null on the light cone. When forming the closed chain, the term quadratic in the 
leading terms drops out, implying that A xy ~ r~ 3 . In momentum space, this situation 
corresponds to the fact that the vector g(k) points almost in the same direction as k. 
In other words, the shear of the surface states is small. Thus in (I2.14p and ()2.18p as 
well as in g + , the shear is small, implying that the vector component is null on the 
light cone, explaining the cancellation of the term ~ r~ 2 in (I2.24H . The states in (I2.22p 
and however, have a large shear. Thus the corresponding vector component is not 
null on the light cone, explaining the term ~ r~ 2 in (I2.25p . We point out that states of 
large shear have never been considered before, as in [4] we always assumed the shear 
to be small. For simplicity, we refer to the states in (I2.22P and g- as shear states. 

We next outline how the above considerations can be adapted to the general an- 
satze (jl.8p and (jl.9p . In order to describe several chiral Dirac seas, one simply adds 
regularized Dirac seas, each of which might involve a right-handed high-energy compo- 
nent and/or shear states. In other words, in the chiral ansatz (I1.8P one replaces each 
summand by a Dirac sea regularized as described above. In the massive ansatz (ll.9p . 
we regularize every massive Dirac sea exactly as described in [H Chapter 4]. More- 
over, in order to distinguish the neutrino sector from a massive sector, we add one or 
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several right-handed high-energy contributions. In this way, the regularization breaks 
the chiral symmetry. 

We finally make a few remarks which clarify our considerations and bring them into 
the context of previous work. 

Remark 2.2. (1) We point out that the above assumption of spherical symmetry 
was merely a technical simplification. But this assumption is not crucial for 
the arguments, and indeed it will be relaxed in £ 12.51 We also point out that 
in all previous regularizations, the occupied states formed a hypersurface in 
momentum space. In this paper, we will always restrict attention to such 
surface states (see [H §4.3]). The underlying guiding principle is that one 
should try to build up the regularized fermionic projector with as few occupied 
states as possible. This can be understood from the general framework of causal 
variational principles as introduced in [51 [8]. Namely, in this framework the 
minimum of the action decreases if the number of particles gets large iQ. Thus 
to construct minimizers, one should always keep the number of particles fixed. 
Conversely, one could also construct minimizers by keeping the action fixed 
and decreasing the number of particles. With this in mind, a regularization 
involving fewer particles corresponds to a smaller action and is thus preferable. 

(2) It is worth mentioning that in all the above regularizations we worked with 
null states, meaning that for every k, the image of the operator P(k) is null 
with respect to the spin scalar product. Such null states can be obtained from 
properly normalized negative definite states by taking a singular mass limit, 
similar as worked out in [H Appendix C]. 

(3) At first sight, our procedure for regularizing might seem very special and ad- 
hoc. However, it catches all essential effects of more general regularizations, 
as we now outline. First, states of large shear could be used just as well for 
the regularization of massive Dirac seas, also in the charged sector. However, 
our analysis in Section [6] will reveal that the EL equations will only involve the 
difference in the regularization used in the charged sector compared to that in 
the neutrino sector. Thus it is no loss in generality to regularize the charged 
sector simply according to [U Chapter 4], and to account for shear states only 
in the neutrino sector. Next, in the high-energy region one could also work 
with massive states. In order to break the chiral symmetry, one could project 
out one spin state with the ansatz 

j(p) = ±(l-pj)W + m)h(k) (2.26) 

with p 2 = m 2 , q 2 = —1 and {q, k) = (see [H eq. (C.1.5)], where a correspond- 
ing Dirac sea is considered before taking the singular mass limit). However, 
this procedure would have two disadvantages. First, massive states would yield 
additional contributions to the fermionic projector, whereas (|2.26p even gives 
rise to bilinear and pseudoscalar contributions, which would all cause technical 

lr To be precise, this results holds for operators in the class (see [5] Def. 2.7]) if the fermionic 
operator is rescaled such that its trace is independent of /. In the formulation with local correlation 
matrices (see [8] Section 3.2]) and under the trace constraint, the canonical embedding > C^ + 1 
allows one to regard a system of / particles as a special system of / + 1 particles. Since varying within 
the set of / + 1-particle systems gives more freedom, it is obvious that the action decreases if / gets 
larger. 
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complications. Secondly, massive states involve both left- and right-handed 
components, which are coupled together in such a way that it would be more 
difficult to introduce a general interaction. Apart from these disadvantages, 
working with massive states does not seem to lead to any interesting effects. 
This is why we decided not to consider them in this paper. 
(4) We mention that for a fully convincing justification of the vacuum fermionic 
projector (jl.5p and of our regularization method, one should extend the state 
stability analysis from [12] to a system of a charged sector and a neutrino sec- 
tor. Since this analysis only takes into account the behavior of the fermionic 
projector away from the light cone, the high-energy behavior of P £ plays no 
role, so that one could simply work with the explicit formula for the unreg- 
ularized fermionic projector (jl.5p . Then the methods of |12| apply to each 
of the sectors. However, the two sectors are coupled by the term \A\ 2 in the 
Lagrangian. The results of this analysis will depend on the value of the La- 
grange multiplier (|2.17p as well as on the choice of all lepton masses (including 
the neutrino masses). Clearly, the details of this analysis are too involved for 
predicting results. For the moment, all one can say is that there is no general 
counter argument (in the spirit of §2. 2D which might prevent state stability. 

2.4. Ruling out the Chiral Neutrino Ansatz. In this section, we give an argument 
which shows that for chiral neutrinos there is no regularization which gives rise to a 
stable minimum of the causal action principle. More precisely, we will show that even 
taking into account the regularization effects discussed in the previous section, it is 
impossible to arrange that the vacuum satisfies the EL equations in the continuum 
limit (|2.16p and (|2.17p . Our argument applies in such generality (i.e. without any 
specific assumptions on the regularization) that it will lead us to drop the ansatz of 
chiral neutrinos (jl.8p , leaving us with the ansatz of massive neutrinos (jl.9p . 

Considering massive neutrinos is clearly consistent with the experimental observa- 
tion of neutrino oscillations. Based on these experimental findings, we could also have 
restricted attention to the ansatz (jl.9p right away. On the other hand, considering also 
chiral neutrinos (jl.8p has the advantage that we can conclude that massive neutrinos 
are needed even for mathematical consistency. This conclusion is of particular interest 
because in the neutrino experiments, the mass of the neutrinos is observed indirectly 
from the fact that different generations of neutrinos are converted into each other. 
This leaves the possibility that neutrinos might be massless, and that the neutrino 
oscillations can be explained instead by modifying the weak interaction. The following 
argument rules out this possibility by giving an independent reason why there must 
be massive neutrinos. 

Recall that the Dirac seas in the charged sector P c , (|1.6p . can be written as 





where T m 2 is the Fourier transform of the lower mass shell 
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Computing this Fourier integral and expanding the resulting Bessel functions gives the 
expansion in position space 

1 /PP 

T m i(x,y) = -^-3 + m<5(£ 2 ) £ (£°) 

+ ^ (log \m 2 e\ + c + in 9(f) e(£ )) + 0(f log^ 2 )) . 

(see [U §2.5] or [7, §4.4]). The point for what follows is that the light-cone expan- 
sion of P m (x,y) involves a logarithmic pole ~ log(^ 2 ). As a consequence, in the EL 
equations (|2. 16|) we get contributions to ()2. 16j) which involve the logarithm of the ra- 
dius \x — y\ (for details see § 5.1] or the weak evaluation formula (|2.32|) below). In 
order to satisfy the EL equations, these logarithmic contributions in the charged sec- 
tor must be compensated by corresponding logarithmic contributions in the neutrino 
sector. 

Now assume that we consider the chiral neutrino ansatz fjl .8|) . Then the light-cone 
expansion of T/v does not involve logarithmic poles (indeed, the distribution Pq can be 
given explicitly in position space by taking the limit m \ in (|2.27j) and (|2.28p ). Thus 
the logarithmic contributions in the radius must come from the high-energy component 
to the fermionic projector. However, as one sees explicitly from the formulas (|2.2ip 
and (I2.23p . the high-energy component is a Laurent series in the radius and does not 
involve any logarithms. This explains why with chiral neutrinos alone it is impossible 
to satisfy the EL equations. 

This problem can also be understood in more general terms as follows. The log- 
arithmic poles of P m (x,y) are an infrared effect related to the fact that the square 
root is not an analytic function (see the discussion of the so-called logarithmic mass 
problem in [H §2.5 and §4.5]). Thus in order to arrange logarithmic contributions 
in the high-energy region, one would have to work with states on a surface with a 
singularity. Then the logarithm in the radius would show up in the next-to leading 
order on the light cone. Thus in order to compensate the logarithms in (I2.28p . the 
contribution by the high-energy states would be just as singular on the light cone as 
the contribution by the highest pole in (12.28|) . Apart from the fact that it seems dif- 
ficult to construct such high-energy contributions, such constructions could no longer 
be regarded as regularizations of Dirac sea structures. Instead, one would have to put 
in specific additional structures ad hoc, in contrast to the concept behind the method 
of variable regularization (see [3J §4.1]). 

The above arguments show that at least one generation of neutrinos must be massive. 
In particular, we must give up the ansatz (|1.8p of chiral neutrinos. Instead, we shall 
always work with massive neutrinos (jl.9p . and we need to assume that at least one of 
the masses fhn is non-zero. 

For clarity, we finally remark that our arguments also leave the possibility to choose 
another ansatz which involves a combination of both chiral and massive neutrinos, i.e. 

A) 3 

P N (x,y) = ^ XL P (x,y)+ Yl P m^,v) with A) €{1,2}. (2.29) 

0=1 /3=A)+1 

The only reason why we do not consider this ansatz here is that it seems more natural to 
describe all neutrino generations in the same way. All our methods could be extended 
in a straightforward way to the ansatz (|2,29D . 
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2.5. A Formalism for the Regularized Vacuum Fermionic Projector. In the 

following sections §2.51 and §2.61 we incorporate the regularization effects discussed 
in §2.31 to the formalism of the continuum limit. Beginning with the vacuum, we re- 
call that in [4, §4.5] we described the regularization by complex factors and 
(see also §5.1]). The upper index n tells about the order of the singularity on 
the light cone, whereas the lower index keeps track of the orders in a mass expan- 
sion. In §2.3[ we considered a chiral decomposition ()2. 13j) and chose the left- and 
right-handed components independently. This can be indicated in our formalism by 
a chiral index c 6 {L,R}, which we insert into the subscript. Thus we write the 
regularization (|2. 13[) and (|2.14p symbolically as 

P £ ^y) = \(xLtT[-2 + XBtT{ R 2 

If the regularization effects of the previous section are not used in the left- or right- 
handed component, we simply omit the chiral index. Thus if we work with general 
surface states or shear states only in the right-handed component, we leave out the 
left-handed chiral index, 

* / — (-1) , „ irri-l) 



P s (x,y) = -(xL$T{->+ XR $T [R0] 
When using the same notation as in the charged sector, we always indicate that we 

(n) 

assume the corresponding regularizations to be compatible. Thus for factors T 
without a chiral index, we shall use the same calculation rules in the neutrino and in 
the charged sector. This will also make it possible to introduce an interaction between 
these sectors (for details see §2.61 and Appendix [A}. If we consider a sector of massive 
neutrinos ([1.9 j) . we first perform the mass expansion of every Dirac sea 



pe _ H V- ^ T {-l+n) A m 2n+1 (n) , * 

r m~ 2 „| J [2n] + „] 1 [2n+l] ^ Z - 6U ) 



2 ^ n\ L ""' J n\ 

n=0 n=0 

and then add the chiral index to the massless component, 

Ad 00 ™2n , , , 00 rT7 2n+l , ^ ( 2 -31) 
, H_ ™ T (-l+n) , \ ^ ™ /„) 

+ 2 n! I 2 "! ^ n! \ 2n +^\ ' 

n=l n=0 

Now the regularization effects of the previous section can be incorporated by in- 
troducing more general factors and t|"\ and by imposing suitable computation 
rules. Before beginning, we point out that the more general factors should all comply 
with our weak evaluation rule 



dt r,(t,0 tl U Io lo = = 77(|^,0 l ° g }J l) , (2.32) 

m-e T (^) ... T (^) T (di) ... T (ds) (i\C\) L £ 

which holds up to 

(higher orders in e/£ ma cro and e/|£|) • (2.33) 

(n) 

Here L is the degree defined by deg T Q = 1 — n, and c reg is a so-called regularization 
parameter (for details see again [H §4.5] or [7, §5.1]). The quotient of products of 
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factors and in (|2.32|) is referred to as a simple fraction. In order to take into 
account the mass expansion (|2.2ip . we replace every factor T^ c ^ by the formal series 

oo 

V — J_T ( - l+n) (2 34) 

n=0 

This notation has the advantage that it resembles the even part of the standard mass 
expansion (|2.30p . In order to get the scaling dimensions right, we inserted a factor 5~ 2n , 
where the parameter 5 has the dimension of a length. The scaling of 5 will be specified 
later (see (|4.2ip . §4.61 and Section [8]). For the moment, in order to make sense of the 
mass expansion, we only need to assume that the 

length scale 5 » e . (2.35) 

But 5 could be much smaller than the Compton wave length of the fermions of the 
system. It could even be on the same scale as the regularization length e. We thus 
replace (I2.3ip by 

pe ( x V ) = y r l l T ( - 1} +y4V- 2_ T (~l+n) 

( 2 - 36 ) 

i$^m 2n (-i+n) ^ m 2n+1 T (n) 
^ 2 ^ n\ I 2 ™! ^ n\ i 2n + 1 \ ' 

n=l n=0 

The effect of large shear can be incorporated in our contraction rules, as we now 
explain. Recall that our usual contraction rules read 

fd n ')\ - 1 (M , >') 



(t$Y u$f)j = 2 (ii + %']) + ( hi § her orders in e /i*D ( 2 - 37 ) 

ii^S^-^-T^+T^ (2.38) 



(and similarly for the complex conjugates, cf. [4, §4.5] or [7J §5.1]). We extend the first 
rule in the obvious way by inserting lower chiral indices. In the second rule we insert 
a factor 5~ 2 , 

*Si T S = - 4 (« t m 1} + ^ ?tl 2) ) • ( 2 - 39 ) 

The factor 5~ 2 has the advantage that it ensures that the factors with square and curly 
brackets have the same scaling dimension (as one sees by comparing (|2,39p with (I2.34p 
or (|2,30p ; we remark that this point was not taken care of in [1] and [7], simply because 

the factors with curly brackets played no role). The term 5~ 2 T^^ can be associated 
precisely to the shear states. For example, in the expression 

1 „ (-i) \ (0 )„(-i) , „,(-lW0) ^(lM-l) 1 ^(-1)^(1) 



g X1 ^? J [0] ^ ^ [fl,0]^ ~ 1 [0] 1 [R,0] ~t~ 1 [0] 1 [R,0] 1 {0} 1 [R,0] §2 [0] 1 {R,0} ' 

the last summand involves an additional scaling factor of r and can thus be used to 
describe the effect observed in (|2.25p . Using again ()2.35p , we can reproduce the scaling 
of the first summand in (|2,25p . 

In the weak evaluation formula (|2,32p . one can integrate by parts. This gives rise 

(n) 

to the following integration-by-parts rules. On the factors T we introduce a deriva- 
tion V by 

vrj n) = T^~ l) . 
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Extending this derivation with the Leibniz and quotient rules, the integration- by-parts 
rules states that 



17-1(01) rpifla) rp(bl) rp(bfj) 

V I — - ° 1=0. (2.40) 



rp(ci) rp(<>y) rp(d\) rp(d & ) 

J- o ' ' ' -Lq 1 o " ' -io 



(a) 



As shown in [4, Appendix E], there are no further relations between the factors T Q 

(n) (n) 

We finally point out that the chiral factors Tj c \ and T) c p y were introduced in such 
a way that the weak evaluation formula (|2.32p remains valid. However, one should 
keep in mind that these chiral factors do not have logarithmic singularities on the 
light cone, which implies that they have no influence on the power k in (|2.32p . This 
follows from the fact that the chiral factors only describe high-energy effects, whereas 
the logarithmic poles are a consequence of the low-frequency behavior of the massive 
Dirac seas (see also the explicit example fl2.21f) and the explanation thereafter). 

2.6. Interacting Systems, Regularization of the Light-Cone Expansion. We 

now extend the previous formalism such as to include a general interaction; for the 
derivation see Appendix [A] For simplicity, we restrict attention to the system (j 1 . 5 1) 
with massive neutrinos (|1.9p and a non-trivial regularization of the neutrino sector by 
right-handed high-energy states. But our methods apply to more general systems as 
well (see Remark 12.31 below) . In preparation, as in [H §2.3] and §4.1] it is helpful to 
introduce the auxiliary fermionic projector as the direct sum of all Dirac seas. In order 
to allow the interaction to be as general as possible, it is preferable to describe the 
right-handed high-energy states by a separate component of the auxiliary fermionic 
projector. Thus we set 



p aux = p£L®p!L, (2-41) 



where 



P iux= (0^)0 and pf ux = 0P m/3 . (2.42) 

0=1 p=i 

Note that p aux is composed of seven direct summands, four in the neutrino and three 
in the charged sector. As the fourth component of the neutrino sector is reserved 
for right-handed high-energy neutrinos (possibly occupying shear or general surface 
states), the corresponding component vanishes without regularization ()2.42p . 

In order to recover P aux from a solution of the Dirac equation, we introduce the 
chiral asymmetry matrix X by 

X = (1 C3 r rcg X r) e 1 C 3 . (2.43) 

Here dimensionless parameter, which we always assume to take values in the 

range 

< T reg < 1 . 

It has two purposes. First, it indicates that the corresponding direct summand involves 
a non-trivial regularization. This will be useful below when we derive constraints for 
the interaction. Second, it can be used to modify the amplitude of the regularization 
effects. In the limit r reg \ 0, the general surface states and shear states are absent, 
whereas in the = 1, they have the same order of magnitude as the regular 

states. 
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Next, we introduce the mass matrix Y by 

Y = — (mi, 7fi2, rhs, 0, mi, m,2, ms) (2.44) 
m 

(here m is an arbitrary mass parameter which makes Y dimensionless and is useful for 
the mass expansion; see also [H §2.3] or [7, §4.1]). In the limiting case r reg \ 0, we 
can then write P aux as 

7 

paux =X t = tX* with t := PmYp ■ (2.45) 

0=1 

In the case r reg > 0, the fourth direct summand will contain additional states. We 
here model these states by a massless Dirac sea (the shear, and general surface states 
will be obtained later from these massless Dirac states by building in a non-trivial 
regularization) . Thus we also use the ansatz (|2.45p in the case r reg > 0. Since t is 
composed of Dirac seas, it is a solution of the Dirac equation 

(»0 - mY) t = 0. (2.46) 

In order to introduce the interaction, we insert an operator B into the Dirac equa- 
tion, 

(i# + B - mY) t = . (2.47) 

Just as explained in [H §2.2] and [TO], the causal perturbation theory defines t in terms of 
a unique perturbation series. The light-cone expansion (see [H §2.5] and the references 
therein) is a method for analyzing the singularities of t near the light cone. This gives 
a representation of t of the form 

oo 

t(x,y) = ^2 rn Pk (nested bounded line integrals) x T^ n \x,y) 

n=— 1 k 

+ P lc (x,y) + P hc (x,y), (2.48) 

where P le (x,y) and P he (x,y) are smooth to every order in perturbation theory. The 
remaining problem is to insert the chiral asymmetry matrix X into the perturbation 
series to obtain the auxiliary fermionic projector with interaction _p aux . As is shown 
in Appendix El the operator p aux can be uniquely defined in full generality, without 
any assumptions on B. However, for the resulting light-cone expansion to involve only 
bounded line integrals, we need to assume the causality compatibility condition 

(i@ + B - mY) X = X* (i@ + B - mY) for all r reg G (0, 1] . (2.49) 

A similar condition is considered in [U Def. 2.3.2]. Here the additional parameter T reg 
entails the further constraint that the right-handed neutrino states must not interact 
with the regular sea states. This constraint can be understood from the fact that 
gauge fields or gravitational fields should change space-time only on the macroscopic 
scale, but they should leave the microscopic space-time structure unchanged. This 
gives rise to conditions for the admissible interactions of the high-energy states. As is 
worked out in Appendix El the gauge fields and the gravitational field must not lead 
to a "mixing" of the right-handed high-energy states with other states. 

Assuming that the causality compatibility condition holds, the auxiliary fermionic 
projector of the sea states p sea is obtained similar to (12.45f) by multiplication with the 
chiral asymmetry matrix. Incorporating the mass expansion similar to (|2.34p leads to 
the following formalism. We multiply the formulas of the light-cone expansion by X 
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from the left or by X* from the right (which as a consequence of f|2.49|) gives the same 
result). The regularization is built in by the formal replacements 

i p T^^m p T^ ] , (2.50) 



rrv 



r T(»)a T V^J_J_ T (fc+n) r „ Bn 

fc=0 

Next, we introduce particles and anti-particles by occupying additional states or re- 
moving states from the sea, i.e. 



P anx (x,y) = P sea (x,y) - —2^*k(x)* k (y) + — ^i(x)^(y) . (2.52) 

k=l 1=1 

For the normalization of the particle and anti-particle states we refer to [H §2.8] and 
§4.3]. Finally, we introduce the regularized fermionic projector P by taking the partial 
trace (see also @J §2.3] or [7] eq. (4.4)]), 

{p) i = £(j3««)gj0 > (2.53) 

a,/3 

where i,j G {1 5 2} is the sector index, whereas the indices a and /3 run over the 
corresponding generations (i.e., a £ {1, . . . 4} if i = 1 and a 6 {1,2,3} if i = 2). We 
again indicate the partial trace of the mass matrices by accents (see [U §7.1] or 
eq. (5.2)]), 

a «,/3,7i,...,7 p _i 

Remark 2.3. (Regularizing general systems with interaction) We now outline how 
the above construction fits into a general framework for describing interacting fermion 
system with chiral asymmetry. Suppose we consider a system which in the vacuum 
is composed of a direct sum of sums of Dirac seas, some of which involve non-trivial 
regularizations composed of right- or left-handed high-energy shear or general surface 
states. Then the interaction can be introduced as follows: To obtain the auxiliary 
fermionic projector, we replace the sums by direct sums. For each Dirac sea which 
should involve a non-trivial regularization, we add a direct summand involving a left- 
or right-handed massless Dirac sea. After reordering the direct summands, we thus 
obtain 

©Pm i )©( XL^o)©( XrPo) (2.55) 
=1 ' l=t\+i e=e 2 +i 
with parameters 1 < i\ < £2 < f max . In order to keep track of which direct summand 
belongs to which sector, we form a partition L\, . . . ,Ln of {1, . . . ,4nax} such that Li 
contains all the seas in the i th sector. Then the fermionic projector of the vacuum is 
obtained by taking the partial trace as follows, 

P l= E E ( paUX )^' 1 V - ( 2 -56) 

The next step is to specify the intended form of the regularization by parame- 
ters r[ eg , . . . , Tp 6g with p £ The rule is that to every left- or right-handed massless 
Dirac sea which corresponds to a non-trivial regularization we associate a parame- 
ter rl eg . Regularizations which we consider to be identical are associated the same 
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parameter; for different regularizations we take different parameters. Introducing the 
chiral asymmetry matrix X, the mass matrix Y, and the distribution t by 

mY = (mi, . . . , m tl ) (0, . . . , 0) (0, . . . , 0) (2.57) 
X = (l,...,l)ex i (l ) ---,l ! r^,...r^)©XK(l,--- ! l,r fe 7 +1 ,...r^) (2.58) 

*=(©^)®( P o)e( Pa), (2.59) 
e=i ' i=h+i £=£ 2 +i 
the interaction can again be described by inserting an operator 23 into the Dirac equa- 
tion (|2.47p . Now the causality compatibility condition (|2.49p must hold for all values 
of the regularization parameters r[ eg , . . . ,T^ eg , thus allowing for an interaction only 
between seas with identical regularization. Using the causal perturbation expansion 
and the light-cone expansion, we can again represent t in the form (|2.48p . The reg- 
ularization is again introduced by setting _p sea = tX* and applying the replacement 
rules (|2.50p as well as 

oo 

reg rpin) , rog 1 1 rp(k+n) 

XL/R Tj 1 -> XL/R Tj l^tf-pk 1 [R/L,2n,j] ' 

k=0 

where the additional index j in the subscript [B/L, 2n,j] indicates that the factors 
corresponding to different parameters r,- must be treated as different functions. This 
means that the basic fractions formed of these functions are all linearly independent 
in the sense made precise in [U Appendix E]. Finally, we introduce particles and 
anti-particles again by (I2.52P and obtain the fermionic projector by taking the partial 
trace (p3Hjh 



2.7. The t-Formalism. In the formalism of the continuum limit reviewed in §2.51 
the regularization is described in terms of contraction rules. While this formulation 
is most convenient for most computations, it has the disadvantage that the effect of 
the regularization on the inner factors is not explicit. The t- formalism remedies 
this shortcoming by providing more detailed formulas for the regularized fermionic 
projector in position space. The formalism will be used in §3.21 §4.61 and §5.21 It will 
also be important for the derivation of the Einstein equations in Section Here we 
introduce the formalism and illustrate its usefulness in simple examples. 

We begin for clarity with one Dirac sea in the charged sector. Then the mass 
expansion gives (cf. (|2.30p ; see also [H §4.5]) 



A Hfl a T (-i+«) . m2n+1 T M 

2 Z-> n \ * 1 [In] ^ n ! 1 [2n+l] 

n=0 ' n=0 



We choose a vector £ which is real-valued, lightlike and approximates £, i.e. 

£ 2 = 0, | = | and £ = £ + (higher orders in e/\£\) . (2.60) 
Replacing all factors £ in P^ by £, we obtain the function P^, 



,2n , „ , . 00 _2n+l 



r m — 2 2_^ , « J [2n] n \ 1 [2n+l]- 

n=0 ' n=0 

Clearly, this function differs from P^ by vectorial contributions. We now want to 
determine these additional contributions by using that the contraction rules (|2.37|) 
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and (|2.38|) hold. It is most convenient to denote the involved vectors by i^) , which we 
always normalize such that 

(|,4^) = 1. (2.61) 
Then the contraction rules (|2.37[) and (|2.38p are satisfied by the ansatz 

00 2r> 



n! 

n=0 



as is verified by a straightforward calculation. To explain the essence of this compu- 
tation, let us consider only the leading contribution in the mass expansion, 



P £ = \ ^ + 1$ + (deg < -1) + Q(m) . (2.63) 

Taking the square, we obtain 

(P £ f = '(I 4 1} > T^T® - <^% 1} > T$T$ + (deg < -2) + 0(m) 

= -^^S? - i^^) T tSi T m + ( de § < - 2 ) + °H • 

The first summand reproduces the contraction rules ()2.37|) and ()2.38p . Compared to 
this first summand, the second summand is of higher order in e/|£|. It is thus omitted 
in the formalism of the continuum limit, where only the leading contribution in e/|£| 
is taken into account (for details see [H Chapter 4]). More generally, when forming 
composite expressions of ()2.62p in the formalism of the continuum limit, only the mixed 

products (£0 , to ) need to be taken into account, whereas the products (l^,^^) 
involving two factors t may be disregarded. With this in mind, one easily sees 
that the ansatz (|2.62p indeed incorporates the contraction rules (|2.37p and (|2.38p . 
Concerning the uniqueness of the representation (|2.62p , there is clearly the freedom to 

change the vectors t n \ as long as the relations (I2.6ip are respected. Apart from this 
obvious arbitrariness, the representation (I2.62P is unique up to contributions of higher 
order in which can be neglected in a weak evaluation on the light cone. 

In order to extend the above formalism to include the regularization effects in the 
neutrino sector, we define by replacing all factors £ in (12.36P by £. Writing 

I*(x, y) = P e m - iXL^ (~T$ + 2$) (2.64) 
" ™ E^ ^ ((* " 1) T $ M + ^ (2-65) 



n=0 

oc 

i 



00 2n 

E G +n) - D + t^) , (2.66) 



n! 

n=l 



a straightforward calculation shows that the contraction rules (|2.37p , ()2.38j) and (|2.39p 
are indeed respected. 

Clearly, the i-formalism is equivalent to the standard formalism of ^2.51 However, 
it makes some computations more transparent, as we now explain. For simplicity, we 
again consider the leading order in the mass expansion (|2,63p and omit all correction 
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terms, i.e. 

p £ ^y) = \i^ l) + ^ 1)T n , , 

1 . __ (2.67) 

P%y,x) = P%x,y)* = -- ir{- l) - T$ . 

Suppose we want to compute the eigenvalues of the closed chain. As we already saw in 
the example (|2.63p . contractions between two factors to are of higher order in e/|£|. 
Thus, in view of the relations (12.60p . it suffices to take into account the mixed terms, 
i.e. 

= \ T^T^ + I i ^ T^T$ + (higher orders in e/\$) . (2.68) 

When taking powers of A xy , any product of the first summand in (|2.68p with the 
second summand in (|2.68p vanishes, because we get two adjacent factors Similarly, 
we also get zero when the second summand is multiplied by the first summand, because 
in this case we get two adjacent factors /. We thus obtain 

where we again omitted the higher orders in e/|£|. Moreover, powers of products of £ 
and ( can be simplified using the anti-commutation relations; for example, 

2 



and applying (|2.6ip together with the fact that ^ is real, we obtain 



This shows that the Dirac matrices in (|2.69p in the first and second summand in (|2.69p 
both have the eigenvalues two and zero. Prom this fact we can immediately read off 
the eigenvalues of (|2.68p to be 



_ (0) (-1) , , _ (-1) (0) 

A + - 1 [o] 1 [o] and [0] 1 [o] • 

Clearly, these formulas were obtained earlier in the usual formalism (for details see [U 
§5.3 and §6.1] or [7J Section 6.1]). But the above consideration gives a more direct 
understanding for how these formulas come about. 

Another advantage is that it becomes clearer how different contributions to the 
fermionic projector influence the eigenvalues. We explain this in the example of a 
left-handed contribution of the form 

P{x,y)~ X Li. (2.70) 

The corresponding contribution to the left-handed component of the closed chain is 
given by 

If we substitute P £ (y,x) according to ()2.67|) . the factor i will be contracted in any 
composite expression either with u or with another factor l. In both cases, we get 
contributions of higher order in e/|£|. Hence we can disregard the factor l, 
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When multiplying with (|2.68p . the product with the second summand vanishes. Even 
more, using the anti-commutation relations, one finds that 

(A xy f i i (A xy f = <«, |) (1 ffii T^T^) p+q . 

This implies that only the eigenvalue \l+ is influenced; more precisely, 



A L+ x - - Uj ? Tf 0] ij and A L _ x . 

Of course, this result is consistent with earlier computations (see for example [7J Proof 
of Lemma 7.4 in Appendix B]). 

3. The Euler-Lagrange Equations to Degree Five 

Before entering the analysis of the EL equations, we briefly recall the basics. Count- 
ing with algebraic multiplicities, the closed chain A xy has eight eigenvalues, which we 
denote by A^L, where n E {1,2}, c E {L,i?} and s E {+,—}. The corresponding 
spectral projectors are denoted by Fncs- In case of degeneracies, we usually omit the 
lower indices on which the eigenvalues do not depend. For example, in the case of the 
four- fold degeneracy Ail + = X2L+ = ^\R+ = ^2R+-, we simply denote the correspond- 
ing eigenvalue by A+ and the spectral projector onto the four-dimensional eigenspace 
by /••. . 

The considerations in the previous section led us to choosing the Lagrange multi- 
plier /j = I (see (|2.17p ). and thus a minimizer P is a critical point of the auxiliary 
action 

S[P) = [[ £[A xy ]d 4 xd 4 y 

J J MxM 

with C according to (II. 

4^*1 = £ w&i 2 - g iW) = ^£ £ (kli-ia^i) 2 . 

n,c,s \n,c,s / n,c,s n',c',s' 

Considering first order variations of P, one gets the EL equations (see [U §3.5] or for 
more details |7J eq. (5.20)]) 

[P,Q]=0, (3.1) 
where the operator Q has the integral kernel (see [U §3.5 and §5.4]) 

1 , dC 

Q(x, y) = 2 £ QA^r C ^(z, y) 



Ei\^ i_I V 1^ 1 
l A ncsl o ^ I'Vc's'l 



I Anc. 



' F%>,P{x,y). (3.2) 



By testing on null lines (see [7J Section 5.2 and Appendix A], one sees that the com- 
mutator (|3.ip vanishes if and only if Q itself is zero. We thus obtain the EL equations 
in the continuum limit 



Q(x, y) = if evaluated weakly on the light cone . (3-3) 



By relating the spectral decomposition of A xy to that of A yx (see [H Lemma 3.5.1]), 
one sees that the operator Q is symmetric, meaning that 

Q(x,y)* = Q(y,x). (3.4) 
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As in [7] we shall analyze the EL equations (|3.3p degree by degree on the light cone. 
In this section, we consider the leading degree five, both in the vacuum and in the 
presence of gauge potentials. In Section 0] we then consider the next degree four. 

3.1. The Vacuum. Applying the formalism of £ j2.5l and £ j2.6l to the ansatz (|1 .5|) . (jl.6p 
and (jl.9p and taking the partial trace, we obtain according to ()2.41|) and (|2.42p for the 
vacuum fermionic projector the expression 

i / 3 £ T rn7 1} + XR Tree $T^, \ 

^> = 5 p°> * S fljM>) +«*<*>. (3.5, 

where we used a matrix notation in the isospin index. Thus 



XlA xv = — XL 



(-lW(-l) 



3 £7^7] 

+ i (deg < 3) + (deg < 2) , 

and the right-handed component is obtained by taking the adjoint. The eigenvalues 
can be computed in the charged and neutrino sectors exactly as in [3 Section 6.1] to 
obtain 



A 2 +L = X2+R = \2-R = A 2 _ L = 9 T^T^ 1} + (deg < 3) (3.6) 



and 



Ai+l = Xi-r = 31$ (Slf 0] 1} + TregT^J}) + (deg < 3) 

Xi+r = ATZ = (sif§ + r reg Tg 0] ) 3^ + (deg < 3) . 

The corresponding spectral projectors can be computed exactly as in [H §5.3 and §6.1] 
or §6] to 



where F± are given by 



F± := \ ( 1 ± ^4 ) + £(deg < 0) + (deg < 0) . (3.8) 



2 V z-z, 

Here the omitted indices of the factors £, z and their complex conjugates are to be 
chosen in accordance with the corresponding factors Tj ^ and To ^ , respectively. In 
the charged sector, this simply amounts to adding indices to all such factors. In 
the neutrino sector, however, one must keep in mind the contributions involving r reg , 

(n) 

making it necessary to keep the factors r^'i, More precisely, setting 

L w= T w )+ r- r S]' (3 - 9) 
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we obtain 

2 X rF± = 1± 

2 XL F± = 1± 



/(-d t (-i) , 1 _ 

'[0] 1 [o] + 3 Tre gMfi,o]- 1 [K,o]' 



(-l) r (-l) rf(-l) 



L 1] ± l) 



4L 



(0) 



/(-I) 3 /- 1 ) T (- 1) + T 



/(-i) T (-i) 

re g*[R,0]- z [R,0] 



with the error terms as in (|3.7p . Moreover, a direct computation shows that (cf. [U 
eq. (5.3.23)]) 



F nc+ P(x,y) = (deg<0) 
'1 



Fi c _ P(x, y) = Xc 







P(x, y) + (deg < 0) 



F 2c „ P(x, y) = Xc (I fj P(x, y) + (deg < 0) 



(3.10) 
(3.11) 

(3.12) 



Evaluating the EL equations (j3.3j) by substituting the above formulas into (|3.2p . we 
obtain the three conditions 



T (-i) T (o) 



L (-i) T (o) 



T (-l)r(0) 



,(-l) T (0) 



| T (-i) T (o)i 



T, 



(-1) 



L (-i) T (o) 



r (-i) r (o) 



T (-i) L (o) 



L (-i) T (o) 



T (-i) L (o) 



T (-i) T (o) 



L (-i) r (o) 



I r (-l) (0) I 

r[o] J [o] I 



D L (o) 



(-1) 



| T (-i) L (o)i 



T, 



(-1) 



(3.13) 
(3.14) 
(3.15) 



These three equations must be satisfied in a weak evaluation on the light cone. 

To summarize, evaluating the EL equations for the fermionic projector of the vac- 
uum (|3.5p . we obtain a finite hierarchy of equations to be satisfied in a weak evaluation 
on the light cone. As the detailed form of these equations is quite lengthy and will not 
be needed later on, we omit the explicit formulas. 



3.2. The Gauge Phases. Let us introduce chiral gauge potentials. As the auxiliary 
fermionic projector ()2.41|) has seven components, the most general ansatz for chiral 
potentials would correspond to the local gauge group U(7)l x \J(7)r. However, the 
causality compatibility conditions (|2.49p reduce the local gauge group to 

U(6)ixU(6) r xU(1)h, (3.16) 

where U(6)l and \J(6)r act on the first and last three components, whereas the 
group U(l)# acts on the fourth component. Similar as in [3 Section 6.2], to degree 
five the gauge potentials describe phase transformations of the left- and right-handed 
components of the fermionic projector, 

P*"*(x,y) ^ UL ^ y ) + XR U R (x,y))P aux (x,y) + (deg < 2) . (3.17) 
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However, as the gauge group (|3.16p is non-abelian, the unitary operators Ujjh now 
involve the ordered exponential (for details see [H §2.5] or [3, Section 2.2]) 

U L/R = Pexp (-if A{ /R . (3.18) 



Substituting the gauge potentials corresponding to the gauge group (|3.16j) and taking 
the partial trace, we obtain 



XlP(x, v) = Xl % \ K U A + (deg < 2) 



XRP(x,y) = XR 



T (-D uk 1 u R 2 \ 



(3.19) 



+ (deg < 2) 



where 



U 



m = {S R 3) e U(6) 



^eU(l), 



and the hat denotes the partial trace, 



(3.20) 



a,0=l 



a,/3=l 



At this point it is important to observe that our notation in f)3. 19|) is oversimplified 
because it does not make manifest that the four matrices Ufj R and Ufj R on the block 
diagonal describe a mixing of three regularized Dirac seas. Thus when the partial 
trace is taken, one gets new linear combinations of the regularized Dirac seas, which 
are then described effectively by the factor TjL l \ The analysis in [6] gives a strong 
indication that an admissible regularization can be obtained only by taking a sum 
of several Dirac seas and by delicately adjusting their regularizations (more precisely, 
the property of a distributional .MP-product can be arranged only for a sum of at 
least three Dirac seas). This means that if we take a different linear combination of 
our three regularized Dirac seas, we cannot expect that the resulting regularization is 
still admissible. In order to avoid this subtle but important problem, we must impose 
that each of the four matrices Uj}, U 22 , U R and U R 2 is a multiple of the identity 
matrix, because only in this case we get up to a constant the same linear combination 
of regularized Dirac seas as in the vacuum (for more details and similar considerations 
see [H Remark 6.2.3] and [7, §9.3]). This argument shows that the matrices U^ 1 , U 22 , 
U R l and U R 2 must be multiples of the identity matrix. The following lemma tells us 
what these conditions mean for Ul and Ur. 



Lemma 3.1. Suppose that 9 C U(6) is a Lie subgroup such that in the standard 
representation on C 6 ; every g £ 9 is of the form 



9 



ah 



* 

cl C 3 



with a,c G 



(3.21) 



where we used a (3 x 3) block matrix notation, and the stars stand for arbitrary 3x3- 
matrices. Then there is a matrix U £ U(3) such that every g G 9 has the representation 



al c s 
bU 



bU* 

Cl C 3 



with 



e U(2) . 



(3.22) 
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In particular, S is isomorphic to a Lie subgroup o/U(2). 

Proof. For any A G T e 9 we consider the one-parameter subgroup V(r) = e %rA (r G 
Evaluating (|3.2ip to first order in r, we find that 

'al c3 Z* 

Z C1(P3 



,4 



with a 3 x 3-matrix Z. Considering (|3.2ip for the quadratic terms in r, we find that 
the matrices ZZ* and Z*Z are multiples of the identity matrix. Taking the polar 
decomposition of Z, we find that 

A = ( a }^ hT J* ] with a, c G R and b G C . (3.23) 
V 6/7 cl C 3y 

Exponentiating, one finds that V(r) is of the required form (|3.22p . but with U de- 
pending on A. 

We next choose two matrices A, A G T e S and represent them in the form (|3.23p 
(where tildes always refer to A) . It remains that to show that U and U coincide up to 
a phase, 

U = e ilp U with p G R . (3.24) 

To this end, we consider the one-parameter subgroup V(r) = e lT ( A + A ) . Evaluat- 
ing (I3.2ip to second order in r, we obtain the condition 

{A,A}= ( dlc3 * ) withd,eGM. 

Writing out this condition using (I3.23p . we find that 

ad + bU* bU = dl c s . (3.25) 

Let us show that there is a parameter tp G R such that (|3.24p holds. If b or b vanish, 
there is nothing to prove. Otherwise, we know from (|3.25p that the matrix U*U is a 
multiple of the identity matrix. Since this matrix is unitary, it follows that U*U = 
e^l C 3, proving flgJH) . □ 

We point out that the matrix U G U(3) is the same for all g G S; this means that U 
will be a constant matrix in space-time. 

Using the representation (I3.22P in (I3.19p . the left-handed component of the fermionic 
projector becomes 

XlP(x, V) = XL% ( "f U \lt S ) + (deg < 2) , (3.26) 

where Ul G U(2), and U M ns £ U(3) is a constant matrix. The matrix U M ns can be 
identified with the MNS matrix in the electroweak theory. In (|3,26p . we still need to 
make sense of the expressions 

Z7 M ns ^ and J/^ns t \o] 1 ■ (3.27) 

Again, the matrix Umns describes a mixing of regularized Dirac seas, now even com- 
bining the seas with different isospin. Since C/ MN s is constant, one can take the point 
of view that we should adjust the regularizations of all six Dirac seas in such a way 
that the expressions in (|3.27j) are admissible (in the sense that the fermionic projector 
has the property of a distributional TWP-product; see [6]). 
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For the right-handed component, the high-energy component 7^1' makes the ar- 
gument a bit more involved. Applying Lemma 13. II to the right-handed component, we 
obtain a representation of the form 



XrP(x,v) = XR y 



+ (deg < 2) 



with (U R ,V) G 17(2) x 17(1) and a fixed matrix U G 17(3). As explained after (f3T20]) . 
our notation is again a bit too simple in that it does not make manifest that the three 
Dirac seas and the right-handed high-energy states will in general all be regularized 
differently, and that only their linear combination is described effectively by the fac- 
tors To ^ . With this in mind, we can repeat the argument after (|3.20p to conclude 
that the relative prefactor of the regularization functions in the upper left matrix entry 
should not be affected by the gauge potentials, i.e. 

U R lT [o] 1] + ^ V 2*$ = K (T^ + r rcg T^JJ) with « G c . 

In particular, one sees that U R l must be a phase factor, and this implies that Ur must 
be a diagonal matrix. Moreover, we find that V = U R l . 

Putting our results together, we conclude that the admissible local gauge group is 

9 = U(2) L x U(l) fl x U(l) fl . (3.28) 

Choosing a corresponding potential {Al, A r , A r ) G u(2) x u(1) x u(1), the interaction 
is described by the operator 

m \ Al Al ^mns i i ( 4-r \ , . 

S = XR U'f/M NS 4f J +Xi U (3 ' 29) 

Thus the U(l)-potentials A R and A^ couple to the right-handed component of the 
two isospin components. The U(2)-potential Al, on the other hand, acts on the left- 
handed components, mixing the two isospin components. The L^MNS-matrix describe 
a mixing of the generations in the off-diagonal isospin components of Al. 

In order to analyze the EL equations to degree five in the presence of the above gauge 
potentials, we need to compute the eigenvalues of the closed chain (see (j3.3|) and (|3.2p ). 
Combining ()3.19j) with the form of the gauge potentials as specified in ()3.28j) and (|3.29p . 
we obtain 



XL P(x, y) = \ X L it T[ ( 0] 1} (J^ ) + (deg < 2) (3.30) 



o fyN jr,(-l) o \ 

XR P(x, y) = - X Rit[ R [°1 + (deg < 2) (3.31) 

Z V V R T [0] J 

with Ul G U{2) and V R ,V R G U{1), where we again used the notation (|3.9|) and 
introduced the complex number 

c = — Umns ■ 
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It follows for the closed chain that 
XLA xy = 9 - XL ("J; 21 -cU?)(\n " )l> l« T-l CU2, 




+ $ (deg < 3) + (deg < 2 



When diagonalizing the matrix (|3.32|) . the factor causes major difficulties 

because it leads to microscopic oscillations in the eigenvectors. Let us explain this 
problem in detail. First, it is convenient to use the i-formalism, because then, similar 
as explained after (|2.68p . the contributions ~ /| and ~ act on different subspaces. 
Thus it remains to diagonalize the 2 x 2-matrices 



E# cU?\ V»T$ L ^) 0_ 



and 



The characteristic polynomial involves square roots of linear combinations of the inner 
matrix elements, describing non-trivial fluctuations of the eigenvalues on the regular- 
ization scale e. Such expressions are ill-defined in the formalism of the continuum 
limit. A first idea for overcoming this idea would be to extend the formalism such as 
to include square roots of linear combinations of simple fractions. However, even if one 
succeeded in extending the continuum limit in this way, it would be unclear how the 
resulting square root expressions after weak evaluation would depend on the smooth 
parameters Uj? and V^ 1 . The basic difficulty is that integrating over the microscopic 
oscillations will in general not preserve the square root structure (as a simple example, 
an integral of the form J y/a + x f{x) dx cannot in general be written again as a 
square root of say the form \Jb + cx). This is the reason why the complications related 
to the factor in (pT32|) seem to be of principal nature. 

In order to bypass this difficulty, we must restrict attention to a parameter range 
where the eigenvalues of the above matrices can be computed perturbatively. In order 
to make the scaling precise, we write 

T reg = (me) Plcs with < p Icg < 2 . (3.33) 

Under this assumption, we know that that the relation 

T^> = L$ (1 + 0((m £ )^)) holds pointwise (3.34) 

(by "holds pointwise" we mean that if we multiply T^) — L^) by any simple frac- 
tion and evaluate weakly (|2.32p . we get zero up to an error of the specified order). 
Making T re g small in this sense does not necessarily imply that the above matrices 
can be diagonalized perturbatively, because we need to compare r reg to the size of the 
off-diagonal matrix elements U]^ and XJ\ ~. As they are given as line integrals over the 
chiral potentials (cf. (|3.18p ). their size is described by 

||Ai 2 ||.|H and H^Hll 
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(where |.| is a Euclidean norm defined in the same reference frame as £). This leads 
us to the following two cases: 

(i) |g|>> PFiFPH ! |g|<< PFFPH " (3 ' 35) 

In fact, the computations are tractable in both cases, as we now explain. 

Case (i). We expand in powers of r reg . We begin with the case r reg = 0. Then in 
the vacuum, (|3.34p implies that the relations (|3.13p - (|3.15p are trivially satisfied. If 
gauge potentials are present, in the above matrices we can factor out the scalar func- 
tions rlyTjL ^ and tL Tjv, respectively. Thus it remains to compute the eigenvalues 
and spectral projectors of the 2 x 2-matrix 

cuf u?){o v£)- (3 - 36) 

Lemma 3.2. The matrix in (|3.36p is normal (i.e. it commutes with its adjoint). 
Moreover, its eigenvalues have the same absolute value. 

Proof. We denote the matrix in ()3.36p by B and write the two factors in (|3.36p in 
terms of Pauli matrices as 

B = (al + iva) e itp (bl + iwa) 

with a,b,ip £ M. and v, w G K 3 . Using the multiplication rules of Pauli matrices, one 
finds that 

e -itp B = (^-^l + ^aw + W + tJAw)^. (3.37) 

A short calculation shows that this matrix is normal. Moreover, the eigenvalues of B 
are computed by 

e tip ((a6 — vw) ± % law + bv + v A w\) . 
Obviously, these eigenvalues have the same absolute value. □ 

We denote the eigenvalues and corresponding spectral projectors of the matrix 
in (I3.36P by v n L and I n . Then, according to the above lemma, 

I fix I = \"2l\ and I* = I n . (3.38) 

For the left-handed component of the closed chain (|3.32p we thus obtain the eigenval- 
ues X n Ls and spectral projectors F n L S given by 

K.Ls = V n L X s , F n L s = XL In F S j (3.39) 

where X± and F s are given by (cf. (I3.6P and (I3.8p . 



A + = 9 rg } T t ( 0] 1} + (deg < 3) , A_ = 9 1} + (deg < 3) (3.40) 
F± = I ( 1 ± i^z) + ^(deg < 0) + (deg < 0) . (3.41) 



2 V z-z 

The spectral decomposition of xrA xv is obtained by complex conjugation 



X n R± = v n R, A± = A n L T = u n L X± , F n L± = F* Rt . (3.42) 

Combining these relations with (|3.38p and (13.40f) , we conclude that all the eigenvalues 
of the closed chain have the same absolute value. Thus in view of (|3.2p . the EL 
equations are indeed satisfied for T reg = 0. In order to treat the higher orders in T reg , 
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one performs a power expansion up to the required order in the Planck length. The 
EL equations can be satisfied to every order in r reg by imposing suitable conditions 
on the regularization functions. Thus one gets a finite hierarchy of equations to be 
satisfied in a weak evaluation on the light cone. 

Case (ii). We perform a perturbation expansion in the off-diagonal elements U^ 1 
and U^ 2 . If we set these matrix elements to zero, we again get a spectral representation 
of the form ([QjjHpnSj) . but now with 

*-(;$). *-(»!) (3-43, 

and 

v 1 = U 1 L 1 Vg, u 2 = U%Vg. 

Since the diagonal elements of any J7(2)-matrix have the same absolute value, it follows 
that ()3.38p again holds. Hence the EL equations are again satisfied in the case U^ 1 = 
= Uj 2 . Expanding in powers of Uf} and U^ 2 again gives a finite hierarchy of equations 
to be evaluated weakly on the light cone, which can again be satisfied by imposing 
suitable conditions on the regularization functions. 

We conclude that to degree five on the light cone, the EL equations can be satisfied 
by a suitable choice of the regularization functions, whenever the EL equations have 
a well-defined continuum limit. Clearly, the detailed computation of admissible reg- 
ularizations is rather involved. Fortunately, we do not need to work out the details, 
because they will not be needed later on. 



4. The Euler-Lagrange Equations to Degree Four 

We now come to the analysis of the EL equations to degree four on the light cone. 
Before beginning, we clarify our scalings. Recall that the mass expansion increases the 

(n) 

upper index of the factors T and thus decreases the degree on the light cone. In view 
of the weak evaluation formula (|2.32j) . the mass expansion gives scaling factors m 2 
Moreover, the parameter T reg gives scaling factors (me) Plcg (see ()3.33[) ) . Unless stated 
otherwise, we shall only consider the leading order in (me) Plcg , meaning that we allow 
for an error term of the form 

(l + 0((me) p ^)) . (4.1) 

Finally, the weak evaluation formulas involve error terms of the form (|2.33|) . Since 
the contributions to the EL equations of degree four on the light cone involve at least 
one scaling factor m 2 e|£| (from the mass expansion) or a factor with the similar scal- 
ing £|£|/^ macro (f rom the light-cone expansion), the factors e/|£| (which arise from 
the regularization expansion) give rise to at least one factor m 2 e 2 , which can be ab- 
sorbed into the error term (I4.ip . Hence, unless stated otherwise, in all the subsequent 
calculations we neglect the 

(higher orders in e/l maao and (me) Prcg ) . 

For ease in notation, in most computations we omit to write out the corresponding 
error term (1 + 0(e/4iacro) + 0((me) Pl ^)). 
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4.1. General Structural Results. We again denote the eigenvalues of the closed 
chain A xy by Ancs- These eigenvalues will be obtained by perturbing the eigenvalues 
with gauge phases as given in (|3.39[) and (|3.42p . As a consequence, they will again 
form complex conjugate pairs, i.e. 



K y R± = KU ■ ( 4 - 2 ) 

As the unperturbed eigenvalues all have the same absolute value (see (|3.39[) . (|3.38p 
and (I3.40P ). to degree four we only need to take into account the perturbation of the 
square bracket in (|3.2p . Thus the EL equations reduce to the condition 



= AQ(*,y):= Wa|A&|1 £ A |A^vl 



A 



A 



xy 

^i^ s P(x,y), (4.3) 

ncs 



where we again evaluate weakly on the light cone and consider the perturbation of 
the eigenvalues to degree two (also, the superscript xy clarifies the dependence of the 
eigenvalues on the space-time points). 

Here the unperturbed spectral projectors F ncs are given explicitly by (13.39P and (13.411) . 
Moreover, the relations (13.10p - (13.12p can be written in the shorter form 

F?t = (deg < 0) , F X J( = i + (deg < 0) . (4.4) 

Combining these relations with the explicit formulas for the corresponding eigenvalues 
(see (|3.39p and (|3.40p ) as well as using (|4.2j) . the EL equations (|4.3p reduce to the 
conditions 

0C 1L = % 2L = %ir = X 2R mod (deg < 4) , (4.5) 

where 

V) ■= 3 3 T$T^ ¥^ (4.6) 

(for more details see the proof of [7, Lemma 7.1]). 

For all the contributions to the fermionic projector of interest in this paper, it will 
suffice to compute A|A^+| in a perturbation calculation of first or second order. Then 
the complex numbers % nc can be recovered as traces of I n with suitable 2 x 2-matrices, 
as the following lemma shows. 

Lemma 4.1. In a perturbation calculation to first order, there are 2 x 2-matrices %l 
and %r such that 

% nc = Tr c2 (/„ % c ) + (deg < 4) . (4.7) 
In a second order perturbation calculation, one can again arrange (|4.7p . provided that 
the gauge phases v nc in the unperturbed eigenvalues (|3.39|) and (|3.42p must not to be 
taken into account and that the perturbation vanishes on the degenerate subspaces in 
the sense that 

F + (AA) F + = 0. (4.8) 

Proof. In view of ([4.6p . it clearly suffices to show that A|A nc+ | can be written as such 
a trace. Writing 

1 



A|A nc +| — — j" ^(AA nc +) A+ + A + (AA nc + / 

and using (I4.2p . one concludes that it suffices to show that 

AA ncs = Tr c2 (/ n £) (4.9) 
for a suitable 2 x 2-matrix B = B(c, s). 
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The linear perturbation is given by 

AA„ CS = Tr(F ncs AA) . 

As the unperturbed spectral projectors involve a factor I n (see (|3.39p and (|3.42p ). this 
is obviously of the form (|4.9p . 

Using (|4.8p . we have to second order 

AA ncs = T —\ ^( F ncs AAF n , c , { _ s) AA) . (4.10) 

n',c> cs 'Vc'(-s) 

Disregarding the gauge phases u cs in (|3.39p and (I3.42j) . we get 

AA ncs = A _\ Tr(F ncs AAF n , c , ( _ s) AA) 

1 



A., - A_ 



Tt(xcI n F a AAF- 8 AA) 



where in the last line we used the form of the spectral projectors in ()3.39p and ()3.42p 
and carried out the sums over n' and c'. This is again of the form (|4.9p . □ 



Instead of analyzing the conditions (14.5p . we shall always analyze the stronger con- 
ditions 



% L (x, y) = X R (x, y) = c(0 l c a . (4.11 



This requires a detailed explanation, depending on the two cases in (|3.35[) . In case (i), 
when the projectors I n are determined by the chiral gauge potentials, the condi- 
tion (|4.1ip can be understood in two different ways. The first, more physical argument 
is to note that the spectral projectors I n of the matrix product (|3.36p depend on the 
local gauge potentials Al and Ar. In order for these potentials to be dynamical, the 
EL equations should not give algebraic constraints for these potentials (i.e. constraints 
which involve the potentials but not their derivatives). This can be achieved by de- 
manding that the conditions (|4.5p should be satisfied for any choice of the potentials. 
In view of (|4.7p . this implies that (|4.1ip must hold. 

To give the alternative, more mathematical argument, let us assume conversely 
that one of the matrices %l or Xr is not a multiple of the identity matrix. Then the 
perturbation calculation would involve terms mixing the free eigenspaces corresponding 
to Aics and \2cs- More precisely, to first order one would have to diagonalize the 
perturbation on the corresponding degenerate subspace. To second order, the resulting 
contribution to the perturbation calculation would look similar to (I4.10p . but it would 
also involve factors of (Ai cs — \2cs)~^ • In both cases, the perturbed eigenvalues would 
no longer be a power series in the bosonic potentials. Analyzing these non-analytic 
contributions in the EL equations (|4.5p , one finds that they must all vanish identically. 
Working out this argument in more detail, one could even derive (I4.1ip from the EL 
equations. 

In case (ii) in (|3.35p . the projectors I n are isospin-diagonal (I3.43p . so that (I4.5P 
only tests the diagonal elements of JC C . Thus at first sight, (I4.1ip seems a too strong 
condition. However, even in this case the condition (|4.1ip can be justified as follows. 
The left-handed gauge potentials modify the left-handed component of the fermionic 
projector by generalized phase transformations. If the involved gauge potential is 
off-diagonal, it makes an off-diagonal components of P(x, y) diagonal and vice versa. 
As a consequence, satisfying (j4.5|) in the presence of off-diagonal gauge potentials is 
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equivalent to satisfying (|4.5p . We will come back to this argument in more detail in 
Section 

In the next lemma, we express AQ(x,y) in terms of the matrices % c . This lemma 
will be needed in Section [6j 

Lemma 4.2. Under the assumptions of Lemma \4-1\ the kernel AQ(x,y) in fl4.3f> has 

the representation 

AQ(X, V)= l -Yl ^ n c ) InXcU 



2 



where 



n L := % L - i Tr c2 (% L + X R ) l c2 (4.12) 
fane? 7£/j is obtained by the obvious replacements L o R). 

Proof. Using (j4.4p in (|4.3p . only the summands with s = — 1 remain. Applying the 
explicit form of the unperturbed eigenvalues ()3.39p and ()3.40j) . the result follows from 
a straightforward calculation. □ 

The stronger condition (14.11|) is then equivalent to demanding that the relations 

n L (x,y) = = TZ R (x,y) (4.13) 
hold in a weak evaluation on the light cone. 

4.2. The Vacuum. We begin by analyzing the eigenvalues of the closed chain in the 
vacuum. As the fermionic projector is diagonal in the isospin index, we can consider 
the charged sector and the neutrino sector after each other. In the charged sector, the 
eigenvalues can be computed exactly as in [U §5.3]. Using the notation and conventions 
in [7], we obtain 

P(x,y)=ji + i -m 2 YY + mY T$ + (deg < 1) 
3 



A xv =~-:M\ SrJrMT 15 + m 2 YY T^T^ + T^T^ 



3i ■■- I --,(-lW0) ^(0) tf „(-l) 



™y[$t^t^-t^t [0] 



+ m 2 Y 2 T^ + (deg < 2). 

A straightforward calculation shows that the closed chain has two eigenvalues A± , both 
with multiplicity two. They have the form 



A + = 9T ( ( J ) T [ ( 7 1) + m 2 (•••) + (deg < 2) 

[ J !rL_ (4.i4) 

A_ = 9T^T^ + (•••) + (deg < 2), 

where (• • • ) stands for additional terms, whose explicit form will not be needed here 
(for details see 01 eq. (5.3.24)]). 
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In the neutrino sector, by using (|2.36p in the ansatz (jl.9p and (|2.42p . after taking 
the partial trace we obtain 

P( x v ) =®$ T { - 1) + yrt ^ (T { - 1] + <T 2 T (0) ) 

+ l -(m 2 YY T$ + mY T$ + (deg < 1) (4.15) 
3 



XLA xy =- a xlU T[ Q] l) (3 + r rcg if-Jf + r rcg 2]^,) (4.16) 

+ lf(rn 2 YY(T^ + T^) + m 2 Y 2 T^ (4.17) 

+ f my^" 1 ^- ^f^) + (deg < 2) . (4.18) 
The contraction rules ([237]) and (pT39j) yield {(l) 2 = (z + z) (% + zz and thus 

m-z)m-z)=o. 

This shows that the matrix has the eigenvalues z and z. Also applying (|3.34p . the 
eigenvalues of the closed chain are computed by 

^L + = \z T^(3L\ ^+r Icg 5^T^ 2] ) + m 2 (• • • ) 

= 9 T^L\^ + 3 r rcg 6- 2 T$T^ } +m 2 (■■■ ) + (deg < 2) (4.19) 



3 



Ar - - T (_1) 2 (3T (_1) + r T (_1) + r T (0 M + f- • • ^ 

AL - ~ 4 i [o] z \ c,1 [o] + Tre s J [R,o] ^ Tre e - t [H,2]/' + V ; 

= 9 - 3 r reg <T 2 T^^Tgoj + m 2 (• • • ) + (deg < 2) , (4.20) 

where iffl is again given by (|3.9|) . and m 2 (• • • ) denotes the same contributions as 



in (|4.14p with the masses mp replaced by the corresponding neutrino masses fhp. The 
two other eigenvalues are again obtained by complex conjugation (|4.2p . 

The first summands in (|4.19p and (|4.20p are of degree three on the light cone and 
were already analyzed in Section [3l Thus the point of interest here are the summands 
involving 5. Before analyzing them in detail, we point out that they arise for two 
different reasons: The term in (|4.19p is a consequence of the mass expansion for general 
surface states. The term in (|4.20p . on the other hand, corresponds to the last term in 
the contraction rule (|2.39j) . which takes into account the shear of the surface states. 

Let us specify the scaling of the terms involving 5. Recall that the parameter T reg 
scales according to (|3.33p . whereas 5 is only specified by (|2.35p . We want that the 
general surface and shear states make up for the fact that the masses mp of the 
charged fermions are different from the neutrino masses fhp. Therefore, it would be 
natural to impose that the summands involving 5 should have the same scaling as the 
contributions m 2 (• • • ) arising in the standard mass expansion. This gives rise to the 
scaling 

T reg _ 2 

Preg 

and thus 5 ~ m (me) 2 . But 5 can also be chosen smaller. In this case, the terms 
involving 5 in (|4.19p and (|4.20p could dominate the contributions by the standard mass 
expansion. But they do not need to, because their leading contributions may cancel 
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when evaluated weakly on the light cone. With this in mind, we allow for the scaling 

1 Preg 

e<^5< — (me)—. (4.21) 
m 

Assuming this scaling, by choosing the regularization parameters corresponding to the 
factors Tjj2 2] and T^ Q y appropriately, we can arrange that (|4.3p holds. This procedure 
works independent of the masses mp and mp. 

4.3. The Current and Mass Terms. We now come to the analysis of the inter- 
action. More precisely, we want to study the effect of the fermionic wave functions 
in (|2.52p and of the chiral potentials (|3.29[) in the Dirac operator (|2,47p on the EL 
equations to degree four. As in [7, Section 7] we consider the contribution near the 
origin in a Taylor expansion around £ = 0. 

Definition 4.3. The integrand in ()2.32p is said to be of order o(|^| fc ) at the origin 

if the function n is in the class o((|£°| + |£|) fc+L )- Likewise, a contribution to the 
fermionic projector of the form P(x,y) x. rj(x,y) is of the order o(|£| fc ) if rj € 

°((l£°l + l£l) fe+1 - n ). ' 

Before stating the main result, we define the bosonic current jjjR and the Dirac 
current J L/R by 

it/R = 9k A/R ~ aA L/R (4-22) 



(44?) = E *\ m XR/Ll k *f> a} - E ^W^f^ ■ ( 4 -23) 
1=1 1=1 



Note that, due to the dependence on the isospin and generation indices, these currents 
are 6 x 6-matrices. We also point out that for the sake of brevity, in (|4.22p we omitted 
the terms quadratic in the potentials which arise for a non-abelian gauge group. But 
as the form of these quadratic terms is uniquely determined from the well-known 
behavior under gauge transformations, they could be inserted into all our equations in 
an obvious way. Similar to the notation (|2.54p . we denote the partial trace over the 
generation indices by j and J. Moreover, we introduce the 2 x 2-matrix-valued vector 
field Zl by 

t L = J k R K 1+ j k L K 2 +j k R K 3 (4.24) 

- 3m 2 (A k L YY + YYA k A K 4 (4.25) 

+ m 2 (A k L YY + YY A k A K 4 (4.26) 

- 3m 2 [A k R YY - 2YA k L Y + YYA k R ^ K 5 (4.27) 

- 6m 2 (A k L Y Y + Y fAfy K 6 (4.28) 
+ 6m 2 (YA k L Y + YA\ Y) K 7 (4.29) 
+ m 2 {A k L YY + 2YA k R Y + YY A k L ^j K 6 (4.30) 

- m 2 (A k R YY + 2YA\Y + YYA k R \ K 7 , (4.31) 
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where K \ , . . . , Kj are the simple fractions 

3 1 



K- 



K- 



K R 



16vr T (0) 



r (_i) T ( ) r (-i) 



c.c. 



2 — T 



,(0) 



T (o) r (o) r (-i) T (o) 



c.c. 



3 — t: 



,(0) L 



r (-i) r (i) r (-i) T (o) 



1 

4 



1 



1 1 

4 



:C l 



(0) (0) (-1) (0) 
' - 1 [2] >] >] 



(0) L 



(-1) (!) (-!) (0) 
' J [2] J [0] >] 



C.C. 



c.c. 



c.c. 



^6 = — 

12 



#7 



! T r L° ) T r ( n 0) 



T (0) T (-i) 



T (-i) T (o) 



,(0) 



T (o) T (-i) 



T (-i) T (o) 



(-!) (-!) 
1 J f0l J [0l 



T (0) T (0) 
i [l] >] 



T (0) T (0) 



12 



,(0) 



T (o) r (-i) 



T (-i) T (o) 



evaluated weakly on the light cone (|2,32p (and c.c. denotes the complex conjugate). 
Similarly, the matrix Zr is defined by the replacements L f-> R. 

Lemma 4.4. The contribution of the bosonic current (I4.22p and of the Dirac cur- 
rent (|4.23p to the order (deg < 4) + o(|£| -3 ) can be written in the form (|4.7p with 

X L/R = iiu t L/R + (deg < 4) + o(|£~T 3 ) . 

Proof. The perturbation of the eigenvalues is obtained by a perturbation calculation 
to first and second order (see [2J Appendix G] and [TJ Appendix B]). The resulting 
matrix traces are computed most conveniently in the double null spinor frame (f ± ) 
with the methods described in [7, Appendix B]. One finds that AA is diagonal on 
the degenerate subspaces, so that the second order contribution is given by (|4.10p . 
Moreover, the gauge phases v nc in the unperturbed eigenvalues (|3.39p and (|3.42p only 
affect the error term o(|£| -3 ). We conclude that Lemma 14.11 applies, and thus %l 
and %r are well-defined. 

In order to compute "JCl/Ri we need to take into account the following contributions 
to the light-cone expansion of the fermionic projector: 



XL P(x, y)x-^XLt€i / V | 1] f L T {0) 

0,2 | Obi 7*^) 



- XL 



im XL & 



YA R r(°> 
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+ im xl J JO, 1 1 0] (Y(djA? R ) - (djA j L ) y) T« 
+ ^XL$tif[l,0\0]YYA* L T<f s ) 
+ ^XLUif[0,l\0]A i L YYT^ 
+ m 2 X L [ V [1,0\0]YY4 L T^ 

J X 

-m 2 XL [ V [0,0\0]Y4 R YTW 

J X 

+ m 2 XL [ V [0,l\0] 4l YYt(1) 

J X 

(for the derivation see [U Appendix B] and [3j Appendix A]; cf. also Appendix B]). 
A long but straightforward calculation (which we carried out with the help of the of 
the C++ program class_commute and an algorithm implemented in Mathematica) gives 
the resultQ 

We finally mention a rather subtle point in the calculation: According to f|3.39|) 
and (|3,42|) . our unperturbed eigenvalues involve gauge phases and can thus be ex- 
panded in powers of A*^. As a consequence, we must take into account contributions 
of the form (|4.10p where the factors AA involve no gauge potentials, but the unper- 
turbed eigenvalues \ ncs are expanded linearly in A™.^. In this case, the corresponding 
contributions involving no factors of can be identified with contributions to the 
eigenvalues in the vacuum in (|4.14p and (|4.19p , (|4.20[) . Using that the vacuum eigenval- 
ues all have the same absolute value, the contributions linear in A^^ can be simplified 
to obtain the formulas for 3c listed above. Another, somewhat simpler method to get 
the same result is to use that the operator Q is symmetric (|3.4p (see [H Lemma 3.5.1]). 
Thus it suffices to compute the symmetric part (AQ(x,y) + AQ(y,x)*)/2 of the op- 
erator AQ as defined by (|4.3p . This again gives the above formulas for without 
using any relations between the vacuum eigenvalues. □ 

Let us briefly discuss the obtained formula for The summands in (|4.24p involve 
the chiral gauge currents and Dirac currents; they can be understood in analogy to 
the current terms in §7.1 and §7.2]. The contributions (|4.25|) - (|4.3ip are the mass 
terms. They are considerably more complicated than in [71 §7.1]. These complications 
are caused by the fact that we here consider left- and right-handed gauge potentials 
acting on two sectors, involving a mixing of the generations. In order to clarify the 
structure of the mass terms, it is instructive to look at the special case of a U(l) 
vector potential, i.e. Al = An = A-\^2 (with a real vector field A). In this case, the 
terms P~25j) and (j4T26|) cancel each other (note that AY ■ ■ ■ = 3 AY ■■■), and (j4T271) 
vanishes. Similarly, the summand f)4.28j) cancels (|4.30|) . and ()4.29p cancels (|4.3ip . Thus 
the mass terms are zero, in agreement with local gauge invariance. 

2 The C++ program class_commute and its computational output were included as ancillary files to 
the arXiv submission of 7 . The Mathematica worksheets were included as ancillary files to this arXiv 
submission. 
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4.4. The Microlocal Chiral Transformation. The simple fractions K% and K§ 
involve factors which have a logarithmic pole on the light cone. Before working 
out the field equations, we must compensate these logarithmic poles by a suitable 
transformation of the fermionic projector. We again work with a microlocal chiral 
transformation as developed in [7J §7.8— §7.11]. As the generalizations to a system of 
two sectors is not straightforward, we give the necessary constructions step by step. 
Before beginning, we mention for clarity that in the following sections ^4.41 and §4.51 
we will construct contributions to P(x, y) which enter the EL equations to degree four 
only linearly. Therefore, it is obvious that Lemma 14. II again applies. 

As in [7J §7.9] we begin in the homogeneous setting and work in momentum space. 
Then the logarithmic poles on the light cone correspond to a contribution to the 
fermionic projector of the form 

P(k) x ( X L h + xr i> R ) s'(k 2 ) e(-k°) , (4.32) 

where the vector components v J L are Hermitian 2 x 2-matrices acting on the sector 
index. In order to generate the desired contribution (|4.32p , we consider a homogeneous 
transformation of the fermionic projector of the vacuum of the form 

P(k) = U(k) P aux (k) U(k)* (4.33) 

with a multiplication operator in momentum space U{k). With the operator U(k) we 
want to modify the states of vacuum Dirac sea with the aim of generating a contribution 
which can compensate the logarithmic poles. We denote the absolute value of the 
energy of the states by O = \k°\. We are mainly interested in the regime m<!!< e^ 1 
where regularization effects play no role. Therefore, we may disregard the right-handed 
high-energy states and write the vacuum fermionic projector according to (jl.5p . (JTTTJ) 
and (|1.9p . Expanding in the mass, we obtain 

paux = y + mY ) 5(k 2 ) @(-k°) - (Jjt + mY) m 2 Y 2 5'(k 2 ) 0(-k°) + (deg < 0) . (4.34) 

For the transformation U (k) in (|4,33p we take the ansatz 

U(k) = 1 + 4= Z(k) with Z = X L^ lj + XRR j lj, (4.35) 

where IJ and R? are 6 x 6-matrices (not necessarily Hermitian) which act on the gen- 
eration and sector indices. This ansatz can be regarded as the linear Taylor expansion 
of the exponential U = exp(zZ/\/^), giving agreement to [7> §7.9] (in view of the fact 
that the quadratic and higher orders of this Taylor expansion dropped out in [7> §7.9], 
for simplicity we leave them out here). Note that the operator U(k) is in general not 
unitary (for details see Remark 14.81 below) . 
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Applying the transformation (|4.33|) and (|4.35p to (|4.34[) . only the isospin matrices 



are influenced. A short calculation gives 

XL U(ljt + mY)U* = X L{H + mY) (4.36) 

+ j= XL - 0) +^XL (lY - Yt!) (4.37) 

+ -XlW +^XlIY^ (4.38) 

XLU{tf + mY) m 2 Y 2 U* = Xh{tf m 2 YY + m 3 YYY) (4.39) 

2 ■ 3 

7777 f * \ / v *\ 7777 / ' v / N *\ 

+ -j= xl [tYY YYfr ) + -j= xl [tYYY - YYYt ) (4.40) 

+ ^XLkY 2 i* + 7 ^-XLtY 3 K (4.41) 



(and similarly for the right-handed component). Let us discuss the obtained contri- 
butions. Clearly, the terms (14.36H and (|4.39j) are the unperturbed contributions. Gen- 
erally speaking, due to the factor 5(k 2 ) in (14.34 j) . the contributions (|4.37p and (14.381) 
are singular on the light cone and should vanish, whereas the desired logarithmic con- 
tribution (|4.32p must be contained in (|4.40p or (I4.4ip . The terms (|4.37p of order SI - 5 
contribute to the EL equations to degree five on the light cone. Thus in order for them 
to vanish, we need to impose that 

L = = R (4.42) 

LY - YL* = = RY — YR* . (4.43) 

The last summand in (|4.38p does not involve a factor ^ and is even. As a consequence, 
it only enters the EL equations in combination with another factor of m, giving rise 
to a contribution of degree three on the light cone (for details see Lemma B.l]). 
With this in mind, we may disregard the last summand in f|4.38[) . Similarly, the last 
summand in ()4.4ip and the first summand in (|4.40[) are even and can again be omitted. 
In order for the second summand in (|4.40[) to vanish, we demand that 

LYYY - YYYL* = = RYYY - YYYR* . (4.44) 

Then it remains to consider the first summand in (|4.38p and the first summand 
in (|4.4ip . We thus end up with a chiral contribution to the fermionic projector of 
the form 

XL P(k) x ^ xlIU* S(k 2 ) Q(-k°) - ^ xlU Y 2 l* S'(k 2 ) @(-k°) . (4.45) 

Note that the conditions (|4.42p - (|4.44p are linear in L and R, whereas the contribu- 
tion (|4.45p is quadratic. 

Before going on, we remark that at first sight, one might want to replace the condi- 
tions (I4.43|) and (|4.44h by the weaker conditions 

LY - YL* = RY- YR* = ivAk) 1 C 2 

(4.46) 

LYYY - YYYL* = RYYY - YYYR* = iv 3 (k) 1 C 2 

involving two real- valued vector fields v± and v 3 . Namely, as the resulting contribution 
to the fermionic projector acts trivially on the isospin index and is symmetric under 
the replacement L <-> R, it perturbs the eigenvalues of the closed chain in a way that 
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the absolute values of all eigenvalues remain equal, so that the EL equations are still 
satisfied. However, this argument is too simple because the gauge phases must to be 
taken into account. For the contributions in (|4.45p . the methods in [3 §7.11] make it 
possible to arrange that the gauge phases enter in a way which is compatible with the 
EL equations. For the contributions corresponding to (|4.46|) . however, it is impossible 
to arrange that the gauge phases drop out of the EL equations. Hence v\ and 1^3 would 
necessarily enter the EL equations. As the scaling factors in (|4.37p and (|4.40p 

give rise to a different ^-dependence, these contributions to the EL equations would 
have a different scaling behavior in the radius. As a consequence, the EL equations 
would only be satisfied if v± = V3 = 0. 

For clarity, we want to focus our attention to the component of ()4.45p which will 
give the dominant contribution to the EL equations. For the moment, we only moti- 
vate in words how this component is chosen; the detailed justification that the other 
components can really be neglected will be given in the proof of Proposition 14. 61 below. 
In the EL equations, the chiral component ()4.45p is contracted with a factor £. This 
means in momentum space that the main contribution of (|4.45p to the EL equations 
is obtained by contracting with a factor k (this will be justified in detailed in the 
proof of Proposition 14.61 below) . Therefore, we use the anti-commutation relations to 
rewrite (j4.45j) as 

We now contract with k to obtain 

P L [k] := P[(k) k } = 1 (2LiL* k l W - k 2 Dlf) 5(k 2 ) G(-k°) 

2 ( 4 - 47 ) 
- ^- (2Li Y 2 Ljtfk* - k 2 V Y 2 tfj 5'{k 2 ) Q{-k°) . 

As the factor k 2 vanishes on the mass shell, we may omit the resulting terms (for 
details see again the proof of Proposition 14.61 below) . We thus obtain 

p L [k] = 1 L[k] L[k}* 5(k 2 ) e(-k°) - ^ L[k] m 2 Y 2 L[k}* 5'(k 2 ) @(-k°) , (4.48) 

where we set L[k] = Lj{k) W (note that L[k] is a 2 x 6-matrix, and the star simply 
denotes the adjoint of this matrix). The right-handed component is obtained by the 
obvious replacements L — ^ R. 

Let us work out the conditions needed for generating a contribution of the desired 
form (|4.32p . Similar as explained in [7, §7.9], the first summand in (|4.48p necessarily 
gives a contribution to the fermionic projector. For this contribution to drop out of the 
EL equations, we need to impose that it is vectorial and proportional to the identity 
matrix, i.e. 

L[k] L[k}* = R[k] R[k}* = co(fc) l c a (4.49) 

with some constant Co (A;). The second summand in (|4.48[) is of the desired form (|4.32p . 
Keeping in mind that we may again allow for a vector contribution proportional to 
the identity, we get the conditions 

L[k] m 2 Y 2 L[k\* = ^ v L [k] + c 2 (k) l c2 

n (4.50) 
R[k] m 2 Y 2 R[k\* = - v R [k] + c 2 (k) l c2 , 
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where we set VL/n[k] = Vy R (k) kj (and c 2 is another free constant). Our task is to solve 
the quadratic equations (|4.49j) and (|4.50p under the linear constraints (|4.42[) — f|4.44|) . 
Moreover, in order to compute the smooth contribution to the fermionic projector, we 
need to determine the expectation values involving the logarithms of the masses 

L[k] (m 2 Y 2 log(mY)) L[k]* and R[k] (m 2 Y 2 log(mK)) R[k]* . (4.51) 

We next describe a method for treating the quadratic equations (I4.49j) and (I4.50|) 
under the linear constraints (|4.42p (the linear constraints ()4.43p and (|4.44|) will be 
treated afterwards). We first restrict attention to the left-handed component and 
consider the corresponding equations in fj4.42H . (I4.49j) and (|4.50p (the right-handed 
component can be treated similarly). We write the matrix L[k] in components, 



L[k]= 11 12 13 14 15 16 , (4.52) 



hi h% l\3 hi I15 he 
hi h2 I23 hi I25 he 

where the matrix entries l a b are complex numbers. We use the linear relations (I4.42p 
to express the third and sixth columns of the matrices by 

la3 = —lal — la2 , he = _ hi ~ hb { a = lj 2) . (4.53) 

This reduces the number of free parameters to 8 complex parameters, which we com- 
bine to the matrix 

tf=r^ with ^ = (Wo2,Wa5). (4-54) 

We introduce on C 4 the scalar product (., .)o as well as the positive semi-definite inner 
product (., .)2 by 

{ip a ^b)a = (L[k]L[k]*)l and ty a , ^2 = (L[k] m 2 Y 2 L[k]*) a b (4.55) 

(where we implicitly use (14.531) to determine the third and sixth columns of We 
represent these scalar products with signature matrices, 

{4>, 4>)q = (tp, S 4>) C i , (tp, 4>)2 = (tp, S 2 4>)ci ■ 
Expressing (., .) 2 in terms of (., .)o, 

( r 4>A)2 = (iP,Sq 1 s 2 (/))o , 

the resulting linear operator S Q ~ 1 S 2 is symmetric with respect to (.,.)o- Thus by 
diagonalizing the matrix S ~ 1 S 2 , one can construct an eigenvector basis e\, . . . , which 
is orthonormal with respect to (., .)o, i.e. 

(e a , eb)o = Sab , (e a ,eb) 2 = ^ a $ab ■ (4.56) 

As the matrices have real-valued entries, we can choose the eigenvectors e a such that 
all their components are real. Moreover, as the matrices So and S 2 are block-diagonal 
in the isospin index, we may choose the eigenvectors such that e\ and e 2 are have 
isospin up, whereas and have isospin down, i.e. 

ei, e 2 = (*, *, 0, 0) , e 3 , e 4 = (0, 0, *, *) 

(where the star stands for an arbitrary real- valued entry). Finally, we always order 
the eigenvalues and eigenvectors such that 

< \x\ < Hi and /U3 < /U4 . (4-57) 
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Moreover, Writing the vectors ip a in (|4.54p in this eigenvector basis, 

4 

^a = X>2e dl (4.58) 

d=l 

we can express (|4.55p in the simpler form 

4 4 

(L[k\ L[k]*)l = YM< > (i[fc]m 2 y 2 L[A;]*)g = X>^ • (4.59) 

Moreover, the linear condition (|4.42|) is satisfied. 

In order to treat the remaining linear constraints (|4,43p and (|4.44j) . we decompose 
the coefficients in (I4.58P into their real and imaginary parts, 

^tft = 4ft + ih tpi (d = l,...,4). 

Considering the diagonal entries of (|4.43p and ([4.44p shows that 

b\ = b\ = bl = b\ = . 

The off-diagonal entries make it possible to express af , af in terms of a\, a\ and b\ , bf in 
terms of b\, b\, leaving us with the eight real parameters a\,a\, ajj, a| and a\,a\,b\,b\. 

In order to simplify the setting, it is useful to observe that all our constraints are 
invariant if we multiply the rows of the matrix in (|4.54j) by phase factors according 
to 

il>! -> e iv ifa., V>2 -> e - ^ ^2 with ^el. (4.60) 

These transformations only affect the off-diagonal isospin components of the left- 
handed matrix in ()4.50p . With this in mind, we can assume that this matrix has 
real components and can thus be decomposed in terms of Pauli matrices as 

L[k}m 2 Y 2 L[k\* = tl + xa l + za 3 . (4.61) 

Using this in f|4.59j) and evaluating the real part of the off-diagonal elements of (|4.49[) , 
one finds that b\ = = 6|, leaving us with the six real parameters a\, a\, a\, a 2 ,, a|, af. 
With these six parameters, we need to satisfy three quadratic relations in (|4.6ip and 
three quadratic relations in (|4.49j) . This suggests that for given parameters Co and C2 
as well as t, x, z, there should be a discrete (possibly empty) set of solutions. 

In preparation, we analyze the case when all potentials are diagonal in the isospin 
index. 

Example 4.5. (isospin-diagonal potentials) Assume that the parameter x in (|4,6ip 
vanishes. Evaluating the real part of the off-diagonal components of (|4.49[) and (|4.6ip . 
one finds that a\ = = a?,. The diagonal components of (|4.49p and (|4.61 j) give the 
quadratic equations 

M2 - Ml " M2 - Mi 

/ 2\2 ~t + Z + COM4 , 2\2 t-Z-lQUs , , 

a 3 = > K) = • ( 4 - 63 ) 

M4 - A*3 fH~ M3 

For these equations to admit solutions, we need to assume the non-degeneracies 

H2 ± Mi and M3 7^ M4 • 
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Then there are solutions if and only if all the squares are non-negative. In view of our 
sign conventions (|4.57p . we obtain the conditions 

Co Ml < t + z < Co fJ-2 and Co fJ-3 < t — z < Co n± ■ (4.64) 

Provided that these inequalities hold, the matrix entries a\, a\, a 2 and a\ are uniquely 
determined up to signs. For any solution obtained in this way, one can compute the 
logarithmic expectation value ([4.5ip . 

In order to analyze the conditions (|4.64p . we first note that changing the constant C2 
corresponds to adding a constant to the parameter t (see (|4.6ip and (|4.50p ). Hence we 
can always satisfy (|4.64p by choosing Co and C2 sufficiently large, provided that 

Ml < M4 an d /U3 < • (4.65) 

If conversely these conditions are violated, it is impossible to satisfy (I4.64p in the 
case z = 0. The physical meaning of the inequalities (|4.65p will be discussed in 
Remark 14.91 b elow . 

In the next proposition, we use a perturbation argument to show that the inequal- 
ities ()4.65p guarantee the existence of the desired homogeneous transformations even 
if off-diagonal isospin components are present. 

Proposition 4.6. Assume that the parameters fj,±, . . . , /Z4 defined by (|4.56p and (|4.57|) 

satisfy the inequalities (14 .65 p . Then for any choice of the chiral potentials vl and vr 
in (|4.32[) . there is a homogeneous chiral transformation of the form ([4.350 such that 



the transformed fermionic projector ([4.330 is of the form 

Pik) = p(k) + { XL i L + XRi R ) rgj (4.66) 

+ (vectorial) l c a 5{k 2 ) (l + Q(n -1 )) (4.67) 

+ (vectorial) 1 C 2 5'(k 2 )(l + O(frl)) (4.68) 

+ (pseudoscalar or bilinear) Vo, 5'(k 2 )(l + 0(J7 -1 )^ (4.69) 

+ (higher orders in e/\^\) . (4.70) 



Before coming to the proof, we point out that the values of the parameters Co and C2 
are not determined by this proposition. They can be specified similar as in §7.9] by 
choosing the homogeneous transformation such that Co is minimal (see also Section [TJ). 
In order to clarify the dependence on Co and C2, we simply added a subscript c to the 
factor Xjgj . Similar to [TJ eq. (8.3)], this factor can be written in position space as 

T M = 3^ ( bg + m 6(e0) ) + °M ' 
where sl, )S i is a real-valued smooth function which depends on the choice of Co and C2. 
In fact, srpd may even depend on the isospin components of vl and vr; but for ease 
in notation we shall not make this potential dependence explicit. 

Proof of Proposition We first show that for sufficiently large Co and C2 , there are 
solutions of (j4.6ip and of the left equation in f[4.49j) . Evaluating the real part of the 
off-diagonal components of (|4.49[) and ([4.610 . we get linear equations in a\ and a 2 , 
making it possible to express a\ and a 2 in terms of a\, a 2 , a|, a§. These relations do 
not involve Cq nor C2. As a consequence, the diagonal components of ([4.490 and ([4.610 
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give a system of equations, which for large parameters Co and C2 are a perturbation of 
the system (|4.62p and (|4.63[) . Hence for sufficiently large Co and C2, there are solutions 
by the implicit function theorem. 

Repeating the above arguments for the right-handed potentials, we obtain matri- 
ces L[k] and R[k] such that (|4,49p and (|4.50p hold. Moreover, it is clear from our 
constructions that (|4.42|) . (|4.43[) and (|4.44[) are satisfied. It remains to go through all 
the contributions (|4.36[) — f|4.41 [) and to verify that they are of the form (|4.66p ^ (|4.70p . 
Clearly, (|4.36p and (|4.39p combine to the summand P(k) in (I4.66p . The contributions 
in (|4.37p vanish due to (|4.42p and (|4.43[) . The second summand in ([4.38P as well as the 
first summand in (|4.40j) are of the form f|4.69j) . The second summand in (|4.40p vanishes 
in view of (|4.44j) . Hence it really suffices to consider the first summand in (|4.38[) and 
the first summand in (|4.4ip . which were combined earlier in (|4.45[) . 

It remains to justify the contraction with the momentum k, which led us to ana- 
lyze (|4.48p . To this end, we need to consider the derivation of the weak evaluation 
formulas on the light cone in [3J Chapter 4] . More precisely, the expansion of the vector 
component in [U eq. (4.4.6)— (4.4.8)] shows that k and £ are collinear, up to errors of 
the order Moreover, the terms in (I4.47P which involve a factor k 2 are again of 

the order smaller than the terms where the factors k are both contracted to L 
or L*. This explains the error term (14.70p . □ 

We remark that the error term (|4.70p could probably be improved by analyzing those 
components of U (k) which vanish in the contraction V (k) k 3 . Here we shall not enter 
this analysis because errors of the order e/|£| appear anyway when evaluating weakly 
on the light cone (12.32P . 

Exactly as in [7, §7.10], one can use a quasi- homogeneous ansatz to extend the above 
methods to a microlocal chiral transformation of the form 

and one introduced the auxiliary fermionic projector is defined via the Dirac equation 

(U^ 1 )* (i^ - mY) U' 1 P aux = . (4.72) 
This gives the following result. 

Proposition 4.7. Assume that the parameters m, . . . ,[14 defined by (|4.56p and (|4.57|) 

satisfy the inequalities (|4.65p . Then for any choice of the chiral potentials vl and vr 
in (|4.32|) . there is a microlocal chiral transformation of the form (|4.71|) such that the 
transformed fermionic projector P := jjp aux jj* { s j th e form 

P(x, y) = P(x, y) + ( XL i L + Xr1r) (l + 0(\£\/e macro )) (4.73) 

+ (vectorial) 1 C 2 (deg < 2) + (pseudoscalar or bilinear) (deg < 1) (4.74) 

+ (smooth contributions) + (higher orders in e/\£\) ■ (4-75) 

We conclude this section with two remarks. 

Remark 4.8. (Unitarity of U) We now explain why it is would be preferable that 
the operator U in the microlocal transformation were unitary, and how and to which 
extent this can be arranged. We begin with the homogeneous setting (|4.33p and (|4.35p . 
As pointed out after (j4.35|) . the operator U as given by (|4.35p is in general not unitary. 
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However, the following construction makes it possible to replace U by a unitary op- 
erator without effecting out results: We first consider the left-handed matrices L J (k). 
Note that our analysis only involved the partial trace L\k\ of these matrices contracted 
with k. Moreover, by multiplying the columns by a phase (|4.60p we could arrange that 
all the components in (|4,5ip were real. In this situation, a straightforward analysis 
shows that there is indeed a Hermitian 6 x 6-matrix whose partial trace coincides 
with (14.51 f) . By choosing the other components of LP (k) appropriately, one can ar- 
range that the matrices L? (k) are all Hermitian, and (|4.51|) still holds. Similarly, one 
can also arrange that the matrices R 3 (k) are Hermitian. Replacing the ansatz (|4.35|) 
by U = exp(iZ/VU), we get a unitary operator. A straightforward calculation shows 
that expanding the exponential in a Taylor series, the second and higher orders of 
this expansion only effect the error terms in Proposition 14.61 (for a similar calculation 
see [3 §7.9]). 

Having arranged that U is unitary has the advantage that the auxiliary fermionic 
projector defined via the Dirac equation (|4.72p is simply given by P aux = UPU* 
(whereas if U were not unitary, the auxiliary fermionic projector would involve un- 
known smooth correction terms; see the similar discussion for local transformations 
in §7.7]). 

In the microlocal setting (|4.7ip . the transformation U will no longer be unitary, even 
if the used homogeneous transformations U(.,Vl/r) are unitary for every v^/r- Thus it 
seems unavoidable that the fermionic projector defined via the Dirac equation (|4.72l) 
will differ from the operator UPU* by smooth contributions on the light cone (see 
also the discussion after [Zj eq. (7.86)]). But even then it is of advantage to choose 
the homogeneous transformations U(., v^/r) to be unitary, because then the correction 
terms obviously vanish in the limit l ma , C ro — > oo. More precisely, a straightforward 
analysis shows that these correction terms are of the order |£|/Anacro- 

Remark 4.9. (Lower bound on the largest neutrino mass) The inequali- 
ties (|4.64p give constraints for the masses of the fermions, as we now explain. Thinking 
of the interactions of the standard model, we want to be able to treat the case when 
a left-handed but no right-handed gauge field is present. In this case, Co is non- 
zero, but the parameter z vanishes for the right-handed component. In view of our 
sign conventions (|4.57p . the first inequality in (|4.64[) implies that Co > 0. Then the 
inequalities (|4,64p yield the necessary conditions ()4.65p . More precisely, the eigenval- 
ues nx, . . . , are given in terms of the lepton masses by (see also [7J eq. (7.73)]) 

flyZ = — (rh\ + 772-2 + ^3 ~F \J ™l + ^2 + ^3 — ^1 ^2 — ^1 ^1 — "^l ^3 



/i.3/4 = - f mf + m\ + m\ =F y mf + m\ + m\ — m\ rn^ — m| — m\ rrij 

The first inequality in ()4.65|) is satisfied once the mass 7723 of the r-lepton is much larger 
than the neutrino masses, as is the case for present experimental data. However, the 
second inequality in (|4.65p demands that the largest neutrino mass 7713 must be at 
least of the same order of magnitude as m2- In particular, our model does not allow 
for a description of the interactions in the standard model if all neutrino masses are 
too small. 

Before comparing this prediction with experiments, one should clearly take into 
account that we are working here with the naked masses, which differ from the physical 
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masses by the contributions due to the self-interaction (with a natural ultraviolet cutoff 
given by the regularization length e). Moreover, one should consider the possibility of 
heavy and yet unobserved so-called sterile neutrinos. 

4.5. The Shear Contributions. We proceed by analyzing the higher orders in an 
expansion in the chiral gauge potentials. Qualitatively speaking, these higher order 
contributions describe generalized phase transformations of the fermionic projector. 
Our task is to analyze how precisely the gauge phases come up and how they enter 
the EL equations. The most singular contributions to discuss are the error terms 

(vectorial) 1 C 2 (deg = 1) (4.76) 

in Proposition 14.71 If modified by gauge phases, these error terms give rise to the 
so-called shear contributions by the microlocal chiral transformation. In the setting 
of one sector, these shear contributions were analyzed in detail in [TJ §7.11]. As the 
adaptation to the present setting of two sectors is not straightforward, we give the 
construction in detail. 

Recall that the gauge phases enter the fermionic projector to degree two according 
to (|3.30p and fj3 . 31 f) . In order to ensure that the error term (|4.76p drops out of the 
EL equations, it must depend on the gauge phases exactly as (]3.30p . i.e. it must be 
modified by the gauge phases to 



Ul 1 cUl 2 \ fVff 

XL \cU? U?) +XR {o Vg 



x (vectorial) l c a (deg = 1) . (4.77) 



Namely, if (14.77P holds, the corresponding contributions to the closed chain involve the 
gauge phases exactly as in (13.32p . and a straightforward calculation using (13.37j) (as 
well as (|3.9p and (|3.34p ) shows that the eigenvalues of the closed chain all have the same 
absolute value. If conversely (|4.77p is violated, then the eigenvalues of the closed chain 
are not the same, and the EL equations will be violated (at least without imposing 
conditions on the regularization functions). We conclude that the transformation 
law (|4.77p is necessary and sufficient for the EL equations to be satisfied to degree 
five on the light cone. 

In order to arrange (|4.77p . we follow the procedure in §7.11] and write down the 
Dirac equation for the auxiliary fermionic projector 

D P aux = , (4.78) 

where T> is obtained from the Dirac operator with chiral gauge fields (see (|2.47l) 
and ()3.29p ) by performing the nonlocal transformation (|4.72p . 

v -.= (u- 1 )* (i$ x + X l4r + xr4 l - mY ) u- 1 ■ 

For clarity, we begin with the case when U is homogeneous (|4.35|) (but the gauge 
potentials A L j R are clearly varying in space-time; the generalization to a microlocal 
transformation will be carried after (|4,82p below). We decompose T> into its even and 
odd components, 

f = £>odd + Deven , 

where 

£>odd = Xl^Xr + XR^Xl and V cvcn = Xl^Xl + XR^XR ■ 

In §7.11] we flipped the chirality of the gauge fields in T> even . As will become clear 
below, we here need more freedom to modify the gauge potentials in 2) even . To this 
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end, we now simply replace the gauge fields in T> even by new gauge fields A e ^ to be 
determined later, 

E xciu-'r^+xLrr+xRrr-^yp-'xc (4.79) 

c=L/R 

We replace the Dirac equation (|4.78p by 

(2) odd + D^ n ) P aux = 0. 

Exactly as in the proof of Proposition 7.12], one sees that the component ~ Q,^ 1 
of P satisfies the Dirac equation involving the chiral gauge potentials A*£rj£. In view 
of (|4.45p and (|4.49|) . we see that the left-handed contribution of (|4.76p is modified by 
the chiral gauge potentials to 

L[k] Pexp ( - i /VlT)^) W • (4-80) 

Thus similar as in ()3.17|) . gauge phases appear. The difference is that the chirality is 
flipped, and moreover here the new potentials A e ^^ enter. A-priori, these potentials 
can be chosen arbitrarily according to the gauge group (|3. 16|) . 

The basic difficulty is that the matrix L[k] is non-trivial in the generation index 
(see (|4.52p - (|4.54p ). Moreover, the gauge potential Al involves the Maki-Nakagawa- 
Sakata matrix t7 M Ns (see (|3,26p ). Therefore, it is not obvious how ()4.80p can be related 
to ()3.26p . But the following construction shows that for a specific choice of A e £ en the 
connection can be made: We denote the two column vectors of L[k]* by £1,^2 £ C 6 . In 
view of (I4.49p . these vectors are orthogonal. We set ei = ^i/||^i|| and = ^/ll^ll and 
extend these two vectors to an orthonormal basis Ci , . . . , t§ of C 6 . We choose A^ en 
such that in this basis it has the form 

AT a (k,x) = ( if {x) ^ )V *\ , (4.81) 

where we used a block matrix representation in the subspaces (ei, Z2, £3) and (e4, es, Cq). 
Here the potentials are chosen as in (]3.29p . and V € U(3) is an arbitrary unitary 
matrix. We point out that the whole construction depends on the momentum k of the 
homogeneous transformation in (I4.80p . as is made clear by the notation A e ^ en (k, x). 
Substituting the ansatz (I4.8ip in (|4.80p and using that the columns of L[k]* are mul- 
tiples of %\ and C4, we obtain 

L[A;]Pexp(-i J* ' (A^) j Z 3 ) L[k\* 

( U} 1 dU} 2 \ rri1 rril * 6321> / U} 1 dU} 2 \ 

where Ufj as in (|3.30p and d = (t±, Vt\)c?- Choosing V such that d coincides with the 
parameter c in (|3.30p , we recover the transformation law of the left-handed component 
in (|4.77j) . Repeating the above construction for the right-handed component (by flip- 
ping the chirality and replacing L[k] by we obtain precisely the transformation 
law P~77j) . 
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In order to get into the microlocal setting, it is useful to observe that the k- 
dependence of A e ^ en can be described by a unitary transformation, 

A% en (k, x) = W(k) A L {x) W{k)* with W(k) e U(6) . (4.82) 

Interpreting W as a multiplication operator in momentum space and Al as a multi- 
plication operator in position space, we can introduce A e ^ en as the operator product 

^ ven = WA L W* . (4.83) 

We point out that the so-defined potential A^ en is non-local. As the microlocal chiral 
transformation is non-local on the Compton scale, one might expect naively that the 
same should be true for ^4|^ en . However, A e ^ en can be arranged to be localized on 
the much smaller regularization scale e, as the following argument shows: The k- 
dependence of W is determined by the matrix entries of L[k]. The analysis in §4.41 
shows that the matrix entries of L[k] vary in k on the scale of the energy e^ 1 (in 
contrast to the matrix Z, which in view of the factor in (|4.35|) varies on the 

scale to). Taking the Fourier transform, the operator W decays in position space on 
the regularization scale. 

This improved scaling has the positive effect that the error term caused by the 
quasilocal ansatz (j4.83j) is of the order e/£ ma _ cro . Hence the gauge phases enter the 
left-handed component of the error term (|4.76p as 

/ [/n cU 12 \ 

XL ^21 ^22 + O^macro)) (vectorial) 1 C 2 (deg = 1) . 

Carrying out a similar construction for the right-handed component, we obtain the 
following result. 

Proposition 4.10. Introducing the potentials A^^ 1 in the flipped Dirac operator (|4.79p 
according to (|4.83p , the error term (|4.76p in Proposition \4- 7| transforms to 



, Ul l cU L 2 \ (V$ 



[l + 0{e/£ macro )) (vectorial) l c a (deg = 1) 



In this way, we have arranged that the EL equations are satisfied to degree five 
on the light cone. Note that the above construction involves the freedom in choosing 
the basis vectors t2,ts,t5,tQ as well as the unitary matrix V in (|4.8ip . This will be 
analyzed in more detail in §5.11 



4.6. The Energy-Momentum Tensor and the Curvature Terms. Considering 
the contribution of the particle and anti-particle wave functions in (|2,52p at the ori- 
gin x = y gives rise to the Dirac current terms as considered in §4.31 (for details see 
also [7, §7.2]). We now go one order higher in an expansion around the origin £ = 0. 
Setting z = (x + y)/2 and expanding in powers of £ according to 

9(x) = V{z - £/2) = *(z) - \ S s djV{z) + o{\£\) 

*(!/) = V(z + e/2) = + ~ tfdjViz) + o(\i\) 
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we can write the contribution by the particles and anti-particles as 

P(x, V)~-g^ Yl Xc ^ k ~ ^) + °($) + ( even contributions) 



c=L/R 

where 



Vmim = -imE^wy*^*^ +^E^^f a) , (4.84) 

a=l 6=1 

and similar to (|2.54p . the hat denotes the partial trace over the generation index. We 
denote the vectorial component by 

rpkl rpkl _j_ rpkl 

Taking the trace over the generation and isospin indices, we obtain the energy-mo- 
mentum tensor of the particles and anti-particles. 

Lemma 4.11. The tensors T*j R) (gjSD, give the following contribution to the matri- 
ces in (|4.7p . 

Kl/r ~ i k K 8 + (deg < 4) + o(|£~r 2 ) , 
where Ks is the simple fraction 

3 1 



16vr T (o) 



(-1) (0) (-1) 
- 1 rni A rni A rni i ( - ( - 



(Vtofe i/tai Kg differs from K\ on page 37 in that the term —c.c. has been replaced 
by +c.c). 

The obvious idea for compensating these contributions to the EL equations is to 
modify the Lorentzian metric. At first sight, one might want to introduce a metric 
which depends on the isospin index. However, such a dependence cannot occur, as the 
following argument shows: The singular set of the fermionic projector P(x, y) is given 
by the pair of points (x, y) with light-like separation. If the metric depended on the 
isospin components, the singular set would be different in different isospin components. 
Thus the light cone would "split up" into two separate light cones. As a consequence, 
the leading singularities of the closed chain could no longer compensate each other in 
the EL equations, so that the EL equations would be violated to degree five on the 
light cone. 

Strictly speaking, this argument leaves the possibility to introduce a conformal factor 
which depends on the isospin (because a conformal transformation does not affect the 
causal structure). However, as the conformal weight enters the closed chain to degree 
five on the light cone, the EL equations will be satisfied only if the conformal factor is 
independent of isospin. 

The above arguments readily extend to a chiral dependence of the metric: If the 
left- and right-handed component of the fermionic projector would feel a different 
metric, then the singular sets of the left- and right-handed components of the closed 
chain would again be different, thereby violating the EL equations to degree five (for 
a similar argument for an axial gravitational field see the discussion in §9.3]) 
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Following these considerations, we are led to introducing a Lorentzian metric gij. 
Linear perturbations of the metric were studied in [21 Appendix B]. The contributions 
to the fermionic projector involving the curvature tensor were computed by 

P(.x,v)x±Rjk?Z k $T(-V (4.85) 

+ e ^ T(0) + * (deg - 1} + (deg < 1} ' (486) 

where Rjk denotes the Ricci tensor (we only consider the leading contribution in an 
expansion in powers of |£|/^macro)- We refer to (|4.85|) and f|4.86j) as the curvature terms. 
More generally, in [11] Appendix A] the singularity structure of the fermionic projector 
was analyzed on a globally hyperbolic Lorentzian manifold (for details see also |13j). 
Transforming the formulas in [13] to the the coordinate system and gauge used 
in [2], one sees that (|4.85p and (|4,86p also hold non-perturbatively. In particular, the 
results in [11] show that, to the considered degree on the light cone, no quadratic or 
even higher order curvature expressions occur. In what follows, we consider f|4.85 j) 
and (|4.86p as a perturbation of the fermionic projector in Minkowski space. This 
is necessary because at present, the formalism of the continuum limit has only been 
worked out in Minkowski space. Therefore, strictly speaking, the following results 
are perturbative. But after extending the formalism of the continuum limit to curved 
space-time (which seems quite straightforward because the framework of the fermionic 
projector approach is diffeomorphism invariant), our results would immediately carry 
over to a globally hyperbolic Lorentzian manifold. 

Let us analyze how the curvature terms enter the eigenvalues of the closed chain. 
We first consider the case when we strengthen (|4.2ip by assuming that 

e < 5 < — {me) ^ (4.87) 
m 

Preg 

(the case 5 ~ m (me) will be discussed in Section [8]) . The assumption (j4.87j) makes 
it possible to omit the terms ~ m 2 Rij. 

Lemma 4.12. The curvature of the Lorentzian metric gives the following contribution 
to the matrices %l/r in (|4.7|) . 

Xl,Xr ~ 24" 4^ RjU ^ k A °y p °( x > y ) ( 488 ) 

+ T jR ok ee{^ °)* M (4.89) 

+ m 2 R jk (deg = 4) + (deg < 4) + o(||j- 2 ) , 
where Kiq is the following simple fraction of degree four, 

16 ~ 32 I PI I PI I W°] P] t°] WO] J ' 

(and Po(x,y) and A® y denote the vacuum fermionic projector and the closed chain of 
the vacuum, respectively). 

Proof. The contribution (|4.85p multiplies the fermionic projector of the vacuum by a 
scalar factor. Thus it can be combined with the vacuum fermionic projector Pq to the 
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expression 

c xy P (x, y) with c xy := 1 + ^ R jk ^£ fc . (4.91) 

Hence the closed chain and the eigenvalues are simply multiplied by a common pref- 
actor, 

a _ r 2 4 o \ - r 2 A o 

■"■xy — ^ X y -"-xy ) ^ncs — ^ X y ^ncs • 

As a consequence, the contribution (|4.85[) can be absorbed in (|4.88 jl . 

The summand (|4.86p is a bit more involved, and we treat it in the t- formalism. The 
closed chain is computed by 

A.y-^Rjki^iT^T^ (4.92) 
+ ^R j kith k T^T$. (4.93) 



Similar as explained for the chiral contribution after (|2.70p . the eigenvalues A nc + are 
only perturbed by (|4.92p . More precisely, 



A j_~ — R-u f j f k T (0) T M) 



and the other eigenvalues are obtained by complex conjugation (|4.2j) . In particular, 
one sees that the eigenvalues are perturbed only by a common prefactor. Combining 
the perturbation with the eigenvalues of the vacuum, we obtain 



A 



nc+ 



(i + 1 R jk ee) xl 



In view of (|4.2p . this relation also holds for the eigenvalues A nc _. We conclude 
that ()4.86p can again be absorbed into (|4.88p . A short calculation using (|4.6p shows 
that the contributions so far combine precisely to (|4.88p . 

It remains to consider the effects of shear and of the general surface terms. The 
shear contribution is described by a homogeneous transformation of the spinors which 
is localized on the scale e (for details see Appendix[X]). Since this transformation does 
not effect the macroscopic prefactor c xy in (I4.9ip . the eigenvalues are again changed 
only by a common prefactor. Hence (|4.85p drops out of the EL equations for the shear 
states. The contributions (]4.92p and (I4.93p . on the other hand, do not involve t, and 
are thus absent for the shear states. We conclude that also (|4.86|) drops out of the EL 
equations for the shear states. 

We finally consider the general surface states. As (I4.92P is a smooth factor times 
the vacuum fermionic projector, the Ricci tensor again drops out of the EL equations. 
For the remaining term (I4.93p . the replacement rule (I2.5ip yields the contribution of 
the general mass expansion 
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As a consequence, 



[i?,0] J 





\ _ 3 o , «>& T ^g T (l) T (-l) Tv J T 1 

IT - 27 R «>* !l£i 7W , r(0) T (l) \ _J_ | T 

L ~ 32 ifc ^ 5 2 V [i? ' 0] [0] [0] / 1 [ 



(-1) |2 / 1 





Similarly, 



giving the result. □ 

4.7. Scalar /Pseudoscalar Potentials, the Higgs Field. As explained in §8.5], 
the Higgs potential of the standard model can be identified with suitable scalar /pseudo- 
scalar potentials in the Dirac equation. As shown in [7J Lemma B.l], the contributions 
by the pseudoscalar potentials to the fermionic projector drop out of the EL equations. 
The scalar potentials, on the other hand, contribute to the EL equations to degree three 
on the light cone. As the detailed computations are rather involved, we postpone the 
analysis of these contributions to a future publication. 

5. Structural Contributions to the Euler-Lagrange Equations 

In this section, we analyze additional contributions to the EL equations to degree 
four on the light cone. These contributions will not enter the field equations, but they 
are nevertheless important because they give constraints for the form of the admissible 
gauge fields and thus determine the structure of the interaction. For this reason, we 
call them structural contributions. 

5.1. The Bilinear Logarithmic Terms. We now return to the logarithmic singular- 
ities on the light cone. In £ j4.3l we computed the corresponding contributions to the EL 
equations to the order o(|£| -3 ) at the origin. In §4.41 we succeeded in compensating the 
logarithmic singularities by a microlocal chiral transformation. The remaining ques- 
tion is how the logarithmic singularities behave in the next order in a Taylor expansion 
around £ = 0. It turns out that the analysis of this question yields constraints for the 
form of the admissible gauge fields, as is made precise by the following proposition. 

Proposition 5.1. Assume that the parameter C2 in (|4.50p is sufficiently large and that 
the chiral potentials in (|3.29p satisfy the conditions 

A\ x - Aft = ±(Af - A% for all space-time points . (5.1) 

Moreover, in case (i) in f)3.35|) we assume that the MNS matrix and the mass matrix 
satisfy the relation 

/ Um NS \ yy = YY ( - ^ Um NS \ 2) 
Umns / \Umns / 

Then one can arrange by a suitable choice of the basis ei, . . . , and the unitary ma- 
trix V in (|4.8ip that the contributions to the EL equations ~ |£| -3 log |£| vanish. 
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If conversely (|5.ip does not hold and if we do not assume any relations between the 
regularization parameters, then the EL equations of order |£|~ 3 log|£| are necessarily 
violated at some space-time point. 

The importance of this proposition is that it poses a further constraint on the form of 
the chiral gauge potentials. 

The remainder of this section is devoted to the proof of Proposition 15.11 Generally 
speaking, our task is to analyze how the gauge phases enter the logarithmic singularities 
of the fermionic projector. We begin with the logarithmic current term 

Xi P aux (x, y) x -2 xl f [0, | 1] f L H T (1) , 

J X 

which gives rise to the last summand in (|4.24p (and similarly for the right-handed 
component; for details see [3 eq. (B.17)~(B.18)] or [H Appendix B]). According to the 
general rules for inserting ordered exponentials (see [3J Definition 2.9] or [H Defini- 
tion 2.5.5]), the gauge potentials enter the logarithmic current term according to 

— 2 xl [ V [0,0\l] Pe-*-/?^*-*)* jl(*)7iPe- i £ ,ji L(v-*)i r« (5.3) 

J X 

(where Pe = Pexp again denotes the ordered exponential ()3. 18|) ) . Performing a Taylor 
expansion around £ = gives 

XLP™ x (x,y) x - 1 XL jj. (£±^) 7i rW (5.4) 

+ 1 xl (4& An + An 4o) + ( 5 - 5 ) 

(note that in (|5.5p it plays no role if the functions are evaluated at x or y because the 
difference can be combined with the error term). 

We arranged by the microlocal chiral transformation that the logarithmic singu- 
larity of (|5.4p is compensated by the second summand in (|4.73|) . Since both (|5.4p 
and (|4.7ip involve the argument (x + y)/2, the logarithmic singularities compensate 
each other even if x and y are far apart (up to the error terms as specified in (]5.5p and 
Proposition 14. 7|) . Thus it suffices to analyze how the gauge phases enter (|4.73p . To this 
end, we adapt the method introduced after (|4.80p to the matrix products in (|4.50p . 
Beginning with the left-handed component, the square of the mass matrix is modified 
by the gauge phases similar to (|5.3p and (|5.4j> . (|5.5p by 

XL m 2y2 ^ Xl f V p e -ij:Al(z- X ) k m 2 Y 2p e -ifj;A^v-z) l dz 
J x 

= xl m 2 Y 2 - l - XL m 2 (A{£j Y 2 + Y 2 4^) + o(\£\) 

(for details see O Section 2 and Appendix A]). When using this transformation law 
in (|4.50p . we need to take into account that, similar to 114. 80 j) . the chiral gauge po- 
tentials A L / R must be replaced by A c ^y L n . Thus we need to compute the expectation 
values 

XLP{x,y) x -~m 2 L[k](AT D [t]Y 2 + Y 2 AT n [ti})L[k}* , (5.6) 

where the square bracket again denotes a contraction, ^ft 611 ^] = (A|^ eri )fc £ fc . 

Again choosing the basis t±, . . . ,Zq, the potential ^4|^ en is of the form (|4.8ip . Now 
we must treat the diagonal and the off-diagonal elements of Al separately. Obviously, 
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the diagonal entries in (|4.8ip map the eigenvectors Z\ and to each other. Hence (|5.6j) 
gives rise to the anti-commutator 

XLP(x,y)^- 1 - {(^ J, { ^,L[k] m 2 Y 2 L[k]^ 

For the off-diagonal elements of Al, the matrices V and V* in (|4.81|) make the situation 
more complicated. For example, the lower left matrix entry in (I4.8ip maps ei to a non- 
trivial linear combination of e4,e5,C6, i.e. for any 6 x 6-matrix B, 

L\k\ B (° °) L\kV -( 114112 (iW + BlV? + BlV " ] ° 

Similarly, the upper right matrix entry in (|4.8ip maps to a non-trivial linear combina- 
tion of ti, t%, 63. As a consequence, the off-diagonal elements of Al yield a contribution 
to (|5.6p of the general form 



XLP(x,y)^Af[i]G{k)+Af[i\G{ky with G = (^12 jLV (5-8) 

Here the 2 x 2-matrix G{k) depends on C2 and vi,\k\ as well as on the choice of the basis 
vectors £2 , £3 , es , e 6 an d the matrix V in (|4.8ip . Counting the number of free parame- 
ters, one sees that G{k) can be chosen arbitrarily, up to inequality constraints which 
come from the fact that V must be unitary and that the entries of the matrix m 2 Y 2 in 
the basis (ci, . . . , ee) cannot be too large due to the Schwarz inequality. These inequal- 
ity constraints can always be satisfied by suitably increasing the parameter C2- For this 
reason, we can treat G(k) as an arbitrary matrix involving three free real parameters. 

The right-handed component of the microlocal chiral transformation can be treated 
similarly. The only difference is that the right-handed gauge potentials are already 
diagonal in view of (I3.29p . Thus we obtain in analogy to (15. 7ft 

XR P(x,y) x -- { j^S , \ v R [k] + c 2 (k) 1 C X (5.9) 

whereas (I5.7P has no correspondence in the right-handed component. 

Comparing (|5.3p with (I5.7p . (|5.8p and (15. 9p . one sees that the transformation laws 
are the same for the diagonal elements of Al and Ar. For the off-diagonal elements 
of Al, we can always choose G(k) such as to get agreement with (|5.3p . We conclude 
that by a suitable choice of G(k) we can arrange that the transformation law (I5.3P 
agrees with (|5,7p . (|5.8p and (|5.9p . As a consequence, the logarithmic poles of the 
current terms are compensated by the microlocal chiral transformation, even taking 
into account the gauge phases to the order o(|£|). 

We next consider the logarithmic mass terms 

XL P aux (x, y) x m 2 xl /\ I 0] YY4 L T« (5.10) 

J X 

[0,0\0]Y4 r YTM (5.11) 
+ m 2 xl T[0, 1 1 0] 4 L YY T« , (5.12) 

J X 



rv 

m 2 X L / [u,U|Ujr# R 
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which give rise to the summand f|4.27|) (for details see [71 eq. (B.29)-(B.31)] or [U 
Appendix B]). Here the gauge phases enter somewhat differently, as we now describe. 

Lemma 5.2. Contracting the logarithmic mass terms (|5.1U|) - (|5.12p with a factor £ 
and including the gauge phases, we obtain 

l -Ti{iixLP^{x,y)) 



im 2 



9 

m 



A{( Zl ) YY - 2YA R (z 2 ) Y + YYA? L (z 3 j) £ T« (5.13) 



((Ai^iAl^YY + 2{A{^)YY{AU k ) + YY ) (5-14) 



+ ^ Y{A^){A R ^ k ) Y rW + (||] 2 ) , (5.15) 

where 

3x + y x + y x + 3y 

*1 = , Z 2 = — > *3 = ■ 

The right-handed component is obtained by the obvious replacements L f-> R. 

Proof. Following the method of O proof of Theorem 2.10], we first choose a special 
gauge and then use the behavior of the fermionic projector under chiral gauge trans- 
formations. More precisely, by a chiral gauge transformation we can arrange that Al 
and A R vanish identically along the line segment xy. In the new gauge, the mass 
matrix Y is no longer constant, but it is to be replaced by dynamical mass matri- 
ces Y L / R {x) (see (3j eq. (2.8)] or [4, eq. (2.5.9)]). Performing the light-cone expansion 
in this gauge, a straightforward calculation yields 

- Tr {i( XL P aux (x, y)) X m 2 (Y L {x) Y R (y) - J (Y L Y R )(z) dz) T« . 

Transforming back to the original gauge amounts to inserting ordered exponentials 
according to the rules in [3, Definition 2.9] or [4, Definition 2.5.5]. We thus obtain 

- Tr {it xl P anx (x, y)) x m 2 Y Pe -< % A ^ Y 
ry (5-16) 

m 2 / Pe" 4 £ A ^ {z ~ x) > YY Pe~ l £ A ^ { ^ dz T« . 



Expanding in powers of £ gives the result. □ 

The contribution (15.13P is the mass term which we already encountered in (j4J27j . In 
contrast to (|5.4p and (|4.7ip , the term (|5.13p does not only depend on the variable (x + 
y)/2. However, taking the partial trace, for a diagonal potential Al we obtain in view 
of (f3T29l) that 

(4(*0 YY + YY A{(z 2 )) b = (a l ( Zi ) + Ai(z 2 )) YY 

= 2AU^±y)^YY + o{\^. 



2 

This makes it possible to compensate the logarithmic singularity of f)5 . 13|) by the 
second term in (|4,73p . up to the specified error terms. For the off-diagonal potentials, 
the situation is more complicated and depends on the two cases in ()3.35p . If we are 
in case (i) and (|5.2p is satisfied, then (|5.17p also holds for the off-diagonal potentials. 
As a consequence, the logarithmic singularity of (|5.13p can again be compensated 
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by the second term in (|4.73p . However, if we do not impose (|5.2p . then it seems 
impossible to compensate the off-diagonal logarithmic terms by a microlocal chiral 
transformation (|4.73p . If we assume instead that we are in case (ii) in (|3.35p . then the 
spectral projectors I n are diagonal (|3.43p . so that off-diagonal potentials are irrelevant 
as they do not enter the EL equations (|4.5p . We conclude that we can compensate the 
logarithmic singularities of (|5.13p in case (i) under the additional assumptions (|5.2p . 
and in case (ii) without any additional assumptions. 

The terms (|5.14p and (|5.15p . on the other hand, are quadratic in the chiral gauge 
potentials. Analyzing whether these terms are compatible with the transformation 
law (|5.7p . (|5.8p and (|5.9p of the microlocal chiral transformation gives the following 
result. 

Lemma 5.3. Consider the component of the fermionic projector which involves a 
bilinear tensor field and has a logarithmic pole on the light cone, 

P(x,y) x (xLhf&y)^ + XRffifavhitj) T {1) (5.18) 

(where h % ^ R is a smooth tensor field acting as a 2 x 2-matrix on the isospin index). 
If (|5.ip holds and C2 is sufficiently large, then one can arrange by a suitable choice of 
the basis ti, ■ ■ ■ , te and of the unitary matrix V in (|4.8ip that 

XLl?i + XRh i j t = Wl& (5.19) 

(where h 13 is a suitable tensor field which acts trivially on the isospin index). If con- 
versely (|5.ip does not hold, then (15.191) is necessarily violated at some space-time point. 

Proof. We first analyze the right-handed component. If (|5.13p is transformed according 
to (|5.9p . we could argue just as for the logarithmic current terms to conclude that the 
contribution of the form (|5.18p vanishes. Therefore, it suffices to consider the terms 
obtained by subtracting from (|5,13p - (|5.15p the term (|5,13p transformed according 
to (|5.9p . This gives the second order terms 

B := ^ [a r [^] 2 YY + 2A R [$ YYA R [$ + YYA R [tf) T« 

- ^ (Ar[£\ YA l [£] Y - YA L [t? Y + YA L [t] Y A R [i]) T« . 
We decompose Al into its diagonal and off-diagonal elements, denoted by 

Am=Ai[$+A° L [t]. 

A straightforward calculation using the identity 

A L \e = A d L \e + Am 2 + {^uiim 

gives 

B = ^ Y {[A R \i\ - A d L [i]f + Am 2 ) y T {1) (5.20) 

~Y [A R [i] - A d L [£\,Am}Y T« . (5.21) 

Clearly, the matrix A£[£] 2 is a multiple of the identity matrix. The matrix (Ar[£] — 
Af^) 2 , on the other hand, is a multiple of the identity matrix if and only if (|5,ip 
holds. The anti-commutator in (|5,2ip vanishes on the diagonal. It vanishes provided 
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that (|5.ip holds. We conclude that the contribution (|5.20j) and (|5.2ip acts trivially on 
the isospin index if and only if (|5.ip holds. In this case, 

B = ^ Y - Ailt]) 2 + A° L [tf) Y T« . (5.22) 

It remains to show that under the assumption (|5.ip . we can arrange that the cor- 
responding left-handed contribution, which we again denote by B, is also a multiple 
of the identity matrix, and that it coincides with (|5.22p . In order to compute the 
relevant left-handed contribution, we subtract from (|5.14p and (|5.15p the term (|5.13p 
transformed according to (|5.7p . 

B = ^{Am, {A L [Z] YY - 2YA R [i) Y + YYA L [£\) } T« 

- ^ {A L [tf YY + 2A L [$ YYA L [S] + YYA L [£] 2 ) T (1) + ^ YA R [tf Y T« . 

In view of ([5.8)1 , we can add contributions which involve A 12 or A 2 ^ ■ A short calculation 
shows that in this way, we can arrange that B is again of the form (|5.22j) . □ 

The next lemma gives the connection to the EL equations. 

Lemma 5.4. The contributions to the EL equations ~ |£|~ 3 log|£| vanish if and only 
if the condition (|5.19]1 holds. 

Proof A direct computation shows that the terms of the form ([5.18)1 contribute to the 
EL equations of the order |£| -3 log |£| unless ([5.19P holds. Therefore, our task is to show 
that it is impossible to compensate a term of the form ([5.18)1 by a generalized microlocal 
chiral transformation. It clearly suffices to consider the homogeneous setting in the 
high-frequency limit as introduced in [7, §7.9]. Transforming to momentum space, the 
contribution (15. 18)1 corresponds to the distribution 

fhijk j S"(k 2 )e(-k°) . (5.23) 

Having only three generations to our disposal, such a contribution would necessarily 
give rise to error terms of the form 

Arf hij k j 5'{k 2 ) 0(-k°) or hij k j 5(k 2 ) Q(-k°) . 

These error terms as large as the shear contributions by local axial transformation 
as analyzed in §7.8], causing problems in the EL equations (for details see [3 
§7.8 and Appendix C]). Instead of going through these arguments again, we here 
rule out (I5.23P with the following alternative consideration: In order to generate the 
contribution (|5,23p . at least one of the Dirac seas would have to be perturbed by a 
contribution with the scaling 

— A y hij k j 8(k 2 - m 2 a ) e(-k°) . 

Due to the factor k 3 , this perturbation is by a scaling factor Q larger than the pertur- 
bations considered in [3 §7.9]. Thus one would have to consider a transformation of 
the form (cf. eq. (7.58)]) 

U = exp (iZ) with Z = 0(0°) . 
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This transformation does not decay in and thus cannot be treated perturbatively. 
Treating it non-perturbatively, the resulting shear contributions violate the EL equa- 
tions. □ 

Combining Lemmas 15.31 and 15.41 gives Proposition 15.11 



5.2. The Field Tensor Terms. We now come to the analysis of the contributions 
to the fermionic projector 

xl p(x, v)~\xl$£ if 7,7, r(°) - xl & jjO, 1 1 0] F l l l3 r(°) 

= \xl J\2a-l)^Fi^ j T^+ i - XL J\ ijkl F»ei 5 l l T^, (5.24) 

which we refer to as the field tensor terms (see [2j Appendix A], |3, Appendix A] and [7J 
Appendix B]; note that here we only consider the phase-free contributions, to which 
gauge phases can be inserted according to the rules in [3]). In [TJ, the field tensor 
terms were disregarded because they vanish when the Dirac matrices are contracted 
with outer factors £. Now we will analyze the field tensor terms in the t- formalism 
introduced in H2.7\ This will give additional constraints for the form of the admissible 
gauge fields (see relation (I5.37D below). 

In this section, the corrections in r reg are essential. It is most convenient to keep the 
terms involving r reg in all computations. We assume that we evaluate weakly for such 
a small vector £ that we are in case (ii) in (|3.35[) (this will be discussed in Section [7|). 
It then suffices to consider the sector-diagonal elements of the closed chain. Moreover, 
by restricting attention to the first or second isospin component, we can compute the 
spectral decomposition of the closed chain in the neutrino sector (n = 1) and the chiral 
sector (n = 2). separately. For a uniform notation, we introduce the notation 



MP 



r£ } if n = 2 



with Lo ^ as given by (j3.9|) . Then the unperturbed eigenvalues are given by 



A nL _ = 9 T^M^ , X nR _ = 9 M^T$ . 
Moreover, using the calculations 



. M (o) r (-i) |r (_1) | 

XL P(x, y) = 3i X L$ [ °' : = 3i xl t 



Kr- T3, no- it \ Mn ^1 rr(0 

^— -| X rP(x, y) = 3ixr$ — (o)]- T [0] 
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in (|4.3p . we can write the EL equations as 

/ 1 \ \T { ~ 1] \ 

A\X nL .\ A|A nV _| Mf = (5.25) 

^ n',c> ' \Mn | 

A| W| " 4 £ A IV C '-| TgJ = . (5.26) 

n',c' ' I [01 I 



Note that in the case r reg = 0, these equations reduce to our earlier conditions (14.51) 
and gSJ). 

Our task is to analyze how (|5.24p influences the eigenvalues A nc _ of the closed chain. 



Lemma 5.5. The field tensor terms (|5.24|) contribute to the eigenvalues A nc _ by 

XnL - x I T (2a _ 1} Trca ( /n ^ 6 1}) j ) (5,27) 

+ I J V Tr c2 (j n ei iW i k (t^) 1 ) T^mF + (deg < 2) (5.28) 
W * ~ £(2a - 1) Tr c2 (l n F| | 4 (^J Mf 7^ (5.29) 
- | jP Tr c2 (j B ^ F% f M^T® + (deg < 2) . (5.30) 



Proof. We first consider the effect of a left-handed field on the left-handed eigenval- 
ues. Every summand in (|5.24p involves a factor £T'°). As the factor /,( ) gives no 
contribution (see (|2.66p ). we regularize (|5.24p in the t- formalism by 

XL P(z, y)~\xL jj2a - 1) i 7i + \xi.f *m F l L j f 7 V Zjjjf (5.31) 

(where the hat again denotes the partial trace). For computing the effect on the 
eigenvalues, we first multiply by the vacuum fermionic projector Po(y,x) to form the 
closed chain. Then we multiply by powers of the vacuum chain (|2.69p and take the 
trace. Since the number of factors i in (|2.69p always equals the number of factors £, and 
taking into account that (|5.3ip vanishes when contracted with a factor £, we conclude 
that the factor Po(y,x) must contain a factor i. In view of (]2.66p . this means that 

we only need to take into account the contribution Po(y,x) X — 3i /[gi -^[op We ^ us 
obtain 

XL A xy x | XL Jj2a - 1) ii Ft 7i iff ^ 

~ 4 y j l 7 7 M ^ pi ■ 

Since the last Dirac factor involves t, this contribution vanishes when multiplied by 
the first summand in (|2.69p . Thus our field tensor term only influences the eigen- 
value X n L-- A short calculation gives (|5.27p and (|5.28p . Similarly, a right-handed field 
only influences the corresponding right-handed eigenvalues by (|5.29p and (|5.30p . The 
result follows by linearity. □ 

Before going on, we remark that the above contributions do not appear in the stan- 
dard formalism of the continuum limit, where all factors £ which are contracted to 
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macroscopic functions are treated as outer factors. In order to get back to the stan- 
dard formalism, one can simply impose that i*L£V = 0. However, this procedure, 
which was implicitly used in [7], is not quite convincing because it only works if the 
regularization is adapted locally to the field tensor. If we want to construct a regular- 
ization which is admissible for any field tensor (which should of course satisfy the field 
equations), then the contributions by Lemma 15.51 must be taken into account. 

Corollary 5.6. Introducing the macroscopic functions 

a n l/r = j f(2a - 1) Tr c2 (in^Ji (5.32) 

the absolute values of the eigenvalues are perturbed by the field tensor terms (I5.24p 
according to 



\\\ T - \~— ' In r T (0) T (-1) — T (-1) T (0) 

{±\\iL-\ - (— jy: \ a nL i rnl l m +a nL l M l r 
l T [0] I 



A|A nfi „| = T^-fc {a nR M^Mi- X) +a^M^Mi 0) ) 
Proof. Writing the result of Lemma 15.51 as 



AA„ L _ = 2a nL T^M^ , AA ni? „ = 2a nR M^T$ , 

we obtain 

A|A " i -' = ^km^i a ^ T ^ M ^) - wwM^^) ■ 

|-t[ ] M n | I [0] I 

The calculation for A|A n #_| is analogous. □ 

After these preparations, we are ready to analyze the EL equations f|5.25j) and f|5.26|) . 
We begin with the = 0. Then we can set M$ = T®, giving the condi- 



tions (14.51). where now 



| T (-l) 



<Y _ A|X I 1 B 1 A/f(°) - n T^T^T^V 4-75— t^T^-VtW 
•A-nc — t±\*nL-\ I (o) | n ~ nc [0] 1 [0] 1 [0] " nc 1 [0] 1 [0] 1 [0] ' 

n I 

This formula can be simplified further with the integration-by-parts rules. Namely, 
applying (j2.40|) . we obtain 



u - V I 1 rnl J rnl J rnl I [0] [0] 1 [0] + J [0] 1 [0] J [0] 



Using this relation, we conclude that 



or _ f 0n — \ T (0) T (-1) T (0) _ _/o„C„ ^_|_qW„ l\ T (») T (-1) T (0) 

•A-nc — I ■zUnc u nc I -i jqi -t j-Qj J- jqj — I -rtc^u nC/ ) -f- OLLiLyu nc j \ J- jqj -t jgj ± tqi . 

If any non-trivial gauge field is present, the four macroscopic functions a nc will not all 
be the same (note that even for a vectorial field which is acts trivially on the isospin 
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index, the contribution term f)5.33j) has opposite sign for a n L and a n n). This implies 
that (|4.5p can be satisfied only if we impose the regularization condition 



T (-l) 

>] 


T (0) 


|r(0) 
\ L [0] 


| T (0)| 

l>] 1 


>] - 


T (0) 
>] 



' \|\ 1 -"] 1 -/'"" _ "J 1 rp(0) ( rpWrpi-l) , nn(— 1) rp{0)\ 

- ^|AiL_|— -ttt- i, n] TTZTrfW fOl I ""^ fOl fOl +O n Lir 01 i, n] . (D.di)J 



which must hold in a weak evaluation on the light cone. 

In order to compute the effect of r reg , we first note that the perturbations A|A2 C — | 
do not involve obvious from Corollary \5.6\i . Moreover, the contribution 

of these eigenvalues to (|5.25p and f|5.26|) for n = 2 is independent of r rcg . In view 
of ([5.34p . these contributions drop out of the EL equations. Next, the eigenvalue 
contributes to f|5.25j) and (|5.26[) for n = 2 by 

1 

~ 4 

This is in general non-zero. Thus in order to allow for left-handed gauge fields in 
the neutrino sector, we need to impose additional conditions on the regularization 
functions. The simplest method is to impose that 

r£J 0] = (deg < 1) • 0((me) Plcs ) pointwise . (5.36) 

Then = Ijj?, and flOg) again vanishes as a consequence of (|5.35[) . We remark 
that (|5.36p could be replaced by a finite number of equations to be satisfied in a 
weak evaluation on the light cone. But as these equations are rather involved, we 
here prefer the stronger pointwise condition ()5.36p . We also note that, in contrast to 
the condition ()5.34p for the regularization of ordinary Dirac seas, the relation ()5.36j) 
imposes a constraint only on the right-handed high-energy states. 

It remains to consider the terms involving T^^. These are A|Ai#_| as well as the 

factor \M ( n 1] \ in (^261) in case n = 1. Collecting all the corresponding contributions 
to the EL equations, we get a finite number of equations to be satisfied in a weak 
evaluation on the light cone. Again, we could satisfy all these equations by imposing 
suitable conditions on the regularization. However, these additional conditions would 
basically imply that T^ R i = vanishes, meaning that there are no non-trivial regular- 
ization effects. For this reason, our strategy is not to impose any more regularization 
conditions. Then the EL equations (|5,25p and (|5.26p are satisfied if and only if there 
is no right-handed gauge field in the neutrino sector and if the vectorial component is 
trace- free, 

Tr(hF R ) = and Tr(F L + F R ) = , (5.37) 



because only under these conditions all the equations in involving T^ R ^ or T^ R ^ vanish. 



6. Projection on the Dynamical Degrees of Freedom 

6.1. Preparatory Considerations. Let us briefly review our general strategy for 
deriving the effective field equations. Our starting point is the fermionic projector of 
the vacuum, being composed of solutions of the free Dirac equation (cf. (|2.46[) ) 



(i0 - mY)^ = 
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As explained in £14. 3[ a wave function ^ gives rise to a contribution to the fermionic 
projector of the form (cf. (|2.52j) ) 

- , (6.1) 

which enters the EL equations. Our method for satisfying the EL equations is to 
introduce a suitable potential into the Dirac equation (see (|2.47p ) 

(i^ + !B-m7)f = 0. (6.2) 

After compensating the resulting logarithmic singularities by a microlocal chiral trans- 
formation (see §4.4p . we can hope that the remaining smooth contributions to the 
fermionic projector by the Yang-Mills current can compensate the contribution by the 
Dirac current (|6.ip . This general procedure is worked out for an axial potential in [7] . 

In the present setting of two sectors, the situation is more involved because the po- 
tential 23 in (16. 2p must satisfy the additional constraints (see (13.290 and the conditions 
arising from the structural contributions in Section [5]) . In order to explain the basic 
difficulty, let us consider the simplified situation where the potential must be vectorial 
and isospin diagonal; more precisely, 

S = 4a 3 (6.3) 

with a vector field A. This constraint for the form of the potential implies that the 
YM current must also be of a special form. As a consequence, we cannot expect to 
compensate an arbitrary contribution of the form (16. 2D . The best we can hope for is 
to compensate the vectorial and isospin diagonal component of (|6.2p . suggesting the 
resulting field equation 

d jk A k - OAj = c^y/* (6.4) 

with a coupling constant c. However, this choice of field equation is not unique. More 
generally, one could try to compensate another component of (|6.2p . leading to the field 
equation 

d jk A k - OAj = ^{ca) (6.5) 

with a fixed vector c G M 3 . In general terms, one can say that, due to the constraints, 
we can satisfy the EL equations only "in the direction of" the admissible variations 
of the Dirac operator. This situation bears some some similarity with the Lagrange 
multiplier method as used in minimization problems with constraints. But for the 
Lagrange multiplier method, one needs to know what the "orthogonal complement" of 
the admissible variations are. For our variational principle, however, there is no notion 
of an "orthogonal complement". In order to give this notion a meaning, one would 
have to compute how the action changes in the direction of non-admissible variations. 
It is not clear how this computation could be carried out, because it seems to involve 
unknown regularization details. Even if such computations could be carried out, they 
would certainly be long and tedious and seem impractical. 

As a way out, additional input is obtained by demanding that the resulting system 
of equations should be consistent in the sense that it admits non-trivial solutions. 
In order to explain the method in the above simplified example (|6.3p . we take the 
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divergence of (|6.5p and use (|6.2p together with (|6,3p to obtain 
= d^(ca)-fj^ + H(ca)~jjdj^ 

= i (W - m)^(ca)jj^ - ~ty(co)(i$ - m)*J (6.6) 

= —i (^jfLa 3 (ca)jj^ - ~$>(ca)4& 3i &) ■ 
The resulting equation admits non-trivial solutions only if 

a 3 (ca) = (ca) a 3 . 

This is our consistency condition. It reduces the freedom to choosing the field equations 
within the two-parameter family c = cl + da 3 . In order to understand the significance 
of such consistency conditions, one should keep in mind that more complicated systems 
of equations give rise to many more consistency conditions. For instance, additional 
consistency conditions are obtained by demanding that the total energy-momentum 
tensor be divergence free. More complicated consistency conditions arise if off-diagonal 
potentials and the MNS mixing are taken into account. Therefore, it is to be expected 
that consistency conditions drastically reduce the freedom in formulating the field 
equations. 

Since it is not know how to systematically derive all consistency conditions, our 
method is to try to recover the field equations as the Euler-Lagrange equations corre- 
sponding to an underlying "effective Lagrangian." We now explain this method in the 
above simplified example, whereas the general construction will be given in the next 
section. We first note that the Dirac equation (|6.2p can be obtained from the Dirac 
action action 

JU{i$ + 'Z-mY)y d 4 x (6.7) 

by varying Vl/. In order to derive the field equations, we choose a kinetic term for the 
bosonic potential, 

J (djA k ) (&A k ) d 4 x , (6.8) 

which is chosen such that varying A gives the left side of (|6.5p . The "effective action" 
is now chosen as a linear combination of (|6.7p and (|6.8p . As an overall constant is 
irrelevant, this method only involves one free constant. The resulting field equation 
is necessarily of the form (|6.4p . Working with the effective action has the advantage 
that the conserved quantities obtained via Noether's theorem immediately guarantee 
consistency of the equations (for example, (|6.6p is a consequence of current conserva- 
tion). 

6.2. Projection on the Dynamical Degrees of Freedom. Quantities like currents 
and fields take values in the Hermitian 6 x 6-matrices and have a left- and right-handed 
component. Thus, taking the direct sum of the two chiral components, it is useful to 
introduce the real vector space 

6 6 := Symm(C 6 ) Symm(C 6 ) , (6.9) 

where Symm(C 6 ) denotes the Hermitian 6 x 6-matrices. For example, the Dirac cur- 
rent (|4.23p can be regarded as an element of &q, 

d ■= (Jl,Jr) g ©e 

(here we disregard the tensor indices, which will be included later in a straightfor- 
ward way). Similarly, the Yang-Mills current (|4.22p can be regarded as an element 
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of 66- However, it can take values only in a subspace of ©6, as we now make precise. 
We denote the local gauge group corresponding to the admissible gauge potentials 
by S C U(6)l x U(6)i? and refer to it as the dynamical gauge group (recall that in §3.21 
we found the group (|3.28p . and that taking into account the additional constraints 
encountered in Section [5l the dynamical gauge group is a proper subgroup of (|3.28p ). 
The dynamical gauge potentials are elements of the corresponding Lie algebra q = T e S, 
which can be identified with a subspace of &q, the so-called dynamical subspace. The 
dynamical potentials and corresponding Yang- Mills currents take values in the dynam- 
ical subspace, 

A := (A L ,A R ) G and (JlJr) € . (6.10) 
In our context, the usual Dirac Lagrangian takes the form 

^Dirac = + B - mY)V . (6.11) 

By varying the Dirac wave functions, one gets the Dirac equation (|2.47p . Conversely, 
the Dirac equation (12. 47ft forces us to choose the Dirac Lagrangian according to (16. lip . 
The point is that by varying the Yang-Mills potential and the gravitational field, the 
Dirac Lagrangian also determines how the Dirac spinors couple to the bosonic fields. 
Our goal is to rewrite certain components of (|3.3p in a way which is compatible with 
the coupling of the Dirac spinors as determined by the structure of the Dirac action. 
First, we need to modify (I6.1ip in order to build in that the fermionic projector involves 
the partial trace (|2.53p . implying that the coupling of the fermions to the bosonic fields 
is necessarily described after taking the partial trace. To this end, we replace (|6.1ip 
by 

£ D irac = Tr C 6 (t + % - m7)*) , (6.12) 

where f denotes the matrix which takes the partial trace over the generations and acts 
trivially on the isospin index, 




t = 1 1 1 <g> 1 C 2 



Varying the wave function in (I6.12|) yields the equation (ity + "B — my)f |i which 
is weaker than the original Dirac equation because the matrix f is singular (more 
precisely, we only obtain that the sum of the Dirac equations for the three generations 
must be satisfied). Thus it is not possible to derive the Dirac equations from the 
effective action. But the EL equations of the effective action are compatible with the 
Dirac equation. This ensures that our system of field equations and Dirac equations 
is consistent, and that all the conservation laws obtained from Noether's theorem are 
respected. 

In order to restrict attention to the essence of (|6.12p . we first note that taking the 
partial trace gives rise to an operator 

~ : 6 6 -> 6 2 C 6 6 , 

where in the last inclusion we regard a symmetric 2 x 2-matrix as a 6 x 6-matrix 
which acts trivially on the generation index (and &2 denotes similar to (|6.9p the chiral 
Hermitian 2 x 2-matrices). The coupling of the fermionic wave functions to the chiral 
potentials in (|6.12p is described by the bilinear form 

(.|.) : 6 2 x -> R : (2, A) ^ Tr c6 (J L A L + J R A R ) , (6.13) 
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where Jur is the Dirac current (|4.23p . We refer to (.|.) as the fermion-boson pairing. 
Next, we introduce the following subspaces of ©2, 

3 = q and 91 = {<J | (3\A) = VA £ 0} , 

giving rise to the direct sum decomposition 

6 2 = 3®m. 

The space 3, referred to as the interacting subspace, has the interpretation as the 
components of the Dirac current (always after taking the partial trace) which couple 
to the bosonic field. The non-interacting subspace 9T, on the other hand, does not take 
part in the interaction. 

We next analyze how the potentials and currents enter the EL equations (|6.11|) . 
As explained in §4.11 we work with the stronger conditions (|4.13j) . We consider the 
contributions to degree four on the light cone after compensating the logarithmic 
poles and evaluating weakly on the light cone. Writing the left- and right-handed 
components again as a direct sum, we have symbolically 

(ll L (x,y),K R (x,y)) x + $ B {d jk A k - OA,) + $m(^)) , (6.14) 

where the linear operators <£j, &b and <3?Af describe how the fermionic and bosonic 
currents as well as the mass terms enter the matrices IZl and TZr, 



We want to deduce from (|6.14p a set of equations in which the current d enters only 
in terms of the fermion-boson pairing (|6.13p . The obvious idea is to multiply (I6.14|) 
by the inverse of <3? j and then take the fermion-boson pairing with A 3 . Unfortunately, 
this method does not work because $j is in general not invertible. But the method 
can be realized if we multiply instead by a matrix Hj, which is uniquely determined by 
imposing that the non-interacting subspace should not enter the effective Lagrangian. 
More precisely, we make the following assumptions: 

(a) The mapping <3?j|a : 3 — > 62 is injective. 

(b) *j(3) n $j(9t) = {0}. 

(c) The sets &b(b) and <&m(s) are contained in <i> 7(62). 
Under these assumptions, we can introduce a mapping 

5j:$j(6 2 )^3c6 2 

as follows. Suppose that v £ $7(62). According to (b), we can uniquely decompose v 
as 

v = + vrf\ with G &j(3) and v<yi £ $ jffi) ■ 

Using (a), there is a unique vector u £ 3 with Qj(u) = v%. Setting Hj(v) = u defines 
the desired linear mapping. It has the properties 

S j|*j(3) = (*jb) _1 and 2j| $j{gr j) = . 

Next, according to (c), the vector (j6.14D lies in the image of 4>j. Hence we may 
multiply (I6.14p by Hj. Taking the fermion-boson pairing with a potential A' € Q gives 
the necessary conditions 

((Zjo$ J )(3 j ) + (Zjo$ B )(d jk A k -nA j ) + (Ejo$ M )(A j ) A')=0 VA'eg. (6.15) 
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These conditions can be further simplified: If 3 lies in the interacting subspace, 
then (Hj o = 3- On the other hand, if 3 lies in the non-interacting subspace, 

then (Hj o <3?j)(£f vanishes. In the latter case, the fermion-boson pairing (3|.A') = 
also vanishes. Hence we may replace f)6. 15j) by 

3j + (Sj o $ B ) (d jk A k - DAj) + (Sj o $ M ) (a,-) | a'} = o wi' e g . (6.16) 

With the above construction, we have reduced the EL equations (|6.1ip to equations 
corresponding to the dynamical degrees of freedom as described by the subspace g C 
G>6- We refer to this procedure as the projection on the dynamical degrees of freedom. 

The above construction can be adapted in a straightforward way to the gravitational 
field: We first write the contribution of the energy-momentum and the Ricci tensor to 
the EL equations in analogy to (|6.14p symbolically as 

{K L (x,y),K R (x,y)) x (* T (f jk ) + <& curv (i^)) ^* , (6.17) 
where &t and & CUIV are linear mappings 

$ T : ©2 -> ©2 , $curv : R ©2 ■ 

We introduce 3 curv = <(1, 1)> C ©2 and set 9I curv = {3 | 1)) = 0}, giving rise 

to the direct sum decomposition 

©2 *^curv © ^H*curv • 

Similar to the above conditions (a)-(c), we make the following assumptions. 

(d) The mapping <&Tb cur v : 3 CU rv -> ©2 is injective. 

(e) $T(3 curv ) n $ T (^curv) = {0}. 

(f) $ curv (M) C $t(6 2 ). 

Under these assumption, there is a unique linear operator St : 3>t(©2) — ^ 3 curv with 
the properties 

H T|$ T pcurv) = (^Tbcurv)" 1 aIld S T | $ T (Otcurv) = ■ 

Multiplying ()6.17j) by St and taking the fermion-boson pairing with (1, 1), we obtain 
the necessary condition 

(f jk + (St o $ curv ) (R jk ) | (l, 1)> = . (6.18) 

6.3. The Effective Lagrangian. In the equation ()6.16j) we arranged that the Dirac 
current couples to the bosonic potentials just as required from the Dirac action (|6.1ip . 
Thus in order for this equation to be of variational form, we merely need to assume that 
the terms involving the bosonic current and the mass term are variational derivatives 
of a corresponding Lagrangian. A simple sufficient condition, which will always be 
satisfied in our applications, is that the corresponding quadratic forms are symmetric. 
We thus make the following assumption: 

(g) The bilinear forms 

((Ej o $ B )(.) | . } : x g -> R and <(Sj o | .) : g x g -> E (6.19) 

are symmetric. 
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Then we can introduce Yang-Mills Lagrangian by 

£ym = \ ((Ej o $ B )(d jk A k - UAj) | A j ) + 1 ((Hj o $ M )(A) \ A) 

By adding irrelevant divergence terms, we can write this Lagrangian in the equivalent 
but more familiar form 

£ym = -~ <(Hj o ^ B )(d j A k - d k Aj) | (&A k - d k Ai)) + i ((5j o $ M )(.A) | A> . 
Adding the Dirac Lagrangian (|6.1ip gives the effective action 

S cS = / (£ D irac + £ Y m) d 4 x . (6.20) 

J M 

Varying ^ gives the Dirac equation (|2.47p . Varying the Yang-Mills potential gives 
precisely the equations (|6.16p with A £ g. 

In order to incorporate the field equation for the gravitational field into the effective 
Lagrangian, one replaces the integration measure in (|6.20p by y det g d A x and the 
Dirac operator as usual by that in curved space-time. Moreover, one adds the Ein- 
stein-Hilbert action 

5 EH = / ((H r o$ curv )(i2)| (1,1)) V - det gd A x, (6.21) 
JM 

where R denotes scalar curvature. One should keep in mind that this Lagrangian is 
determined only up to terms which contribute to the EL equations (I3.2p to degree three 
or lower. In particular, if the Ricci tensor is a multiple of the metric, the term Rj k ^^ k 
in Lemma 14.121 is of degree one on the light cone, giving rise to a contribution which 
can be absorbed in the error term. In other words, to the considered degree four on 
the light cone, the Ricci tensor is determined only up to multiples of the metric. This 
gives precisely the freedom to add the cosmological Lagrangian 

f 2A i : l4 

/ — v/ — det g a x 
JM k 

for an arbitrary value of the cosmological constant A. In principle, the cosmological 
constant could be determined in our approach by evaluating the EL equations to degree 
three on the light cone. But this analysis goes beyond the scope of the present work. 

We finally point out that the above method determines the effective Lagrangian and 
fixes the coupling constants (except for the cosmological constant). Thus the struc- 
ture of the interaction is described completely by the underlying EL equations (|3.3p 
corresponding to the causal action principle (11.31) . 

6.4. General Remarks. The projection on the dynamical degrees of freedom gives 
a method to extract from the EL equations (|3.3p a system of field equations which 
are of variational form. The wish to rewrite the field equations as the critical points 
of an effective action is also motivated by the fact that we describe the interaction 
by the underlying causal action principle (|1.3p . and it seems natural to assume that 
the variational form of the equations should be preserved when taking the continuum 
limit. However, at present there is no mathematical method to get directly from the 
action (|1.2[) to the above effective action. Our methods make it necessary to first go 
to the EL equations (|3.3p corresponding to (jl.2p . Analyzing these EL equations in the 
continuum limit gives relations in which the variational form is no longer apparent. By 
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demanding that the equations should be again of variational form, we found a unique 
way of rewriting the field equations as coming from an effective Lagrangian. 

Although this procedure seems sufficient for all practical purposes, it is not quite 
satisfying from the theoretical point of view. There are two conceivable methods to 
improve the situation, which we mention as possible directions of future research: 

(A) In order to prove that the EL equations ()3.3|) are again of variational form in 
the continuum limit, one could analyze systematically how the non-dynamical 
degrees of freedom contribute to the EL equations (examples for such non- 
dynamical degrees of freedom are the gauge potentials which were ruled out 
in §3.21 and §31 but there are many other non-dynamical variations which have 
not been considered so far). Similar as in the Lagrange multiplier method, 
it suffices to satisfy the EL equations in the directions transverse to the non- 
dynamical contributions. This analysis, which might make it necessary to 
impose certain conditions on the structure of the regularization, could give a 
detailed justification of the evaluations (|6.16p and (|6.18p . 

(B) It would be desirable to have a mathematical method which would allow to 
take some kind of "continuum limit" of the action (jl.2p to obtain the effective 
action directly (i.e. without analyzing the EL equations). The general difficulty 
is that the causal action principle (jl.3p involves constraints (in particular the 
condition that P should be a projector of fixed rank) which have no obvious 
correspondence in the continuum limit. This difficulty is also reflected in a quite 
different structure of the action principles: In (|1.3f) the action is minimized. 
For the effective action, a minimization makes no sense because the effective 
action is unbounded from below. Instead, one merely considers critical points 
of the effective action. This general difference must be taken into account in 
any approach towards relating (11.21) directly to an effective action. 

7. The Field Equations for Chiral Gauge Fields 

We now use the methods of Section [6] to compute the effective action for the coupling 
of the Dirac field to the gauge fields. In order to determine the dynamical subspace q C 
©6 (defined before (|6.I0p ). we first recall that in §3.21 we derived the admissible local 
gauge group (|3.28p together with the representation of the gauge potentials (|3.29[) . 
In Section we obtained further restrictions for the gauge potentials. Namely, the 
analysis of the bilinear logarithmic terms in §5.11 revealed that the diagonal elements 
must satisfy the constraint (|5.ip . The field tensor terms in §5.371 on the other hand, 
gave us the two linear constraints (|5.37p for the field tensor, which due to gauge 
symmetry we can also regard as constraints for the potentials. Putting these conditions 
together, we conclude that the dynamical gauge potentials must be of one of the two 
alternative forms 



The potentials of the form (|7.ip do do not form a Lie algebra. This means that the 
structure of (|7,ip is not preserved under local gauge transformations corresponding 
to the potentials of the form (|7.ip . As this seems to be inconsistent, we disregard 




or 
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this case. We thus restrict attention to the remaining case (|7.2|) . where q is the Lie 
algebra su(2), which acts on the left-handed component of the spinors and involves a 
MNS mixing matrix. Hence 

3 = {(av,0) | ve M 3 } and 01 = {(a 1, B) \ a G R, B G Symm(C 2 )} . 

Using the results of Section [H it is straightforward to compute <J>j to degree four 
on the light cone. Applying Lemma 14.41 and ()4.24j) . we obtain 

and similarly for the right-handed component. Using this identity in (|4,12p gives 

K L x ii k {j\ - ~ Tr c2 (J| + J k R ) 1^) K X . 
Using this formula in (|6.14p . we obtain 

*j((Jl, Jr)) = ((4, h) - \ Tr c2 (J L + J R ) (1, 1)) K x + (deg < 4) . (7.3) 

A short calculation shows that the conditions (a) and (b) in ^6.21 are satisfied. Fur- 
thermore, since the image of $j consists of all trace- free matrices in 62, it is obvious 
from (|4.12p that the condition (c) is also satisfied. Moreover, in view of (|7.3p . the 
mapping Hj is simply given by 

1 3 

Ej((B L ,B R )) = -—J2^aB R ) (<r Q ,0) . 

Next, the bilinear forms in ()6.19p are symmetric, and thus the condition (g) in §6.31 
also holds. We conclude that the field equations reduce to the three equations 

= Tr(a a % R ) , a = 1,2,3. 

Here %r contains the contributions by the left-handed Dirac and bosonic currents of 
Lemma 14.41 with the logarithmic poles compensated by the microlocal chiral trans- 
formation according to Proposition 14.71 We state the general structure of these field 
equations and discuss them afterwards. 

Theorem 7.1. Writing the S\J(2)-gauge potentials in components 

At= l -Tr{a«A L ) 

(and similarly for the currents), the field equations read 

jt -M 2 a At = c a Jt + (/ [0] * 3lT + (/ [2] * A L ) a , 

where jl and Jl are the currents ()4.22p and (j4.23|) . respectively. The mass parame- 
ters M a and the coupling constants c a satisfy the relations 

Mi = M 2 and a = c 2 . (7.4) 

Finally, the distributions fm\ and f^] are convolution kernels. 

Note that the relations fj7.4|) follow immediately from the symmetry under transfor- 
mations of the two isospin components according to 

(e^ \ T 
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Qualitatively speaking, this theorem can be understood similar to the results in [7]. 
Also, the calculations use exactly the same methods. In particular, the convolution 
kernels fm\ and /[ 2 i are computed and interpreted just as in [?J §8.1 and §8.2]. In view 
of these similarities, we here omit the detailed computations and only point to two 
steps in the computations which are not quite straightforward. First, as mentioned 
after Proposition 14. 6\ the constants Co and c 2 are not determined by this proposition. 
Following the strategy used in 7.9], we can fix these constants by minimizing Co- 
Thus we choose the microfocal chiral transformation in such a way that the vecto- 
rial contribution (|4.49j) to the fermionic projector is as small as possible. Using this 
method, for a given regularization one can also compute the coupling constants and 
the masses similar as in §8.6]. 

The second step which requires an explanation concerns the computation of the 
coupling constants and bosonic masses for a given regularization in the spirit of [7, §8.6]. 
Here one must distinguish the two cases (i) and (ii) in (|3.35p . Which of these cases 
applies depends crucially on the choice of the parameter p reg in (|3,33p . In particular, 
by choosing p reg sufficiently small (and thus the parameter r reg in ()2.43p sufficiently 
large), we can arrange that we are in case (ii). In order to keep the setting as general 
as possible, we deliberately left open which of the cases should be physically relevant. 
We found that all our computations up to and including Section H] apply in the same 
way in both cases. In the analysis of the bilinear logarithmic terms in §5.1} however, 
our constructions apply in case (i) only under the additional assumption (|5.2p . The 
analysis of the field tensor terms in §5.21 was carried out only in case (ii) (and at 
present it is unclear how the results could be extended to case (i)). This gives a strong 
indication that the physically relevant scaling should indeed be described by case (ii). 
This scaling can be realized by choosing the parameter p reg in (I3.33P sufficiently small. 
Thus in a physical model, the the parameter r reg in (|2.43p should be chosen sufficiently 
large. 

Arranging in this way that we are in case (ii), it remains to justify the transition 
from the EL equations (|4.5p to the stronger conditions (]4.1ip . We already indicated 
an argument before Lemma 14.21 We are now in the position to make this argument 
precise: Recall that in case (ii), the spectral projectors I n are isospin-diagonal (|3.43p . 
The perturbation of these spectral projectors by the gauge phases leads to a finite 
hierarchy of equations to be satisfied in a weak evaluation on the light cone. With this 
in mind, it suffices to satisfy (|4.5p with I n according to ()3.43p . But clearly, we must 
take into account that the gauge phases enter the matrices % nc , as we now explain. We 
begin with the Dirac current terms. As the left-handed component of a wave function 
is modified by the gauge phases in the obvious way by 



the gauge potential enters the Dirac current term as described by the replacement 



In this way, the off-diagonal components of the Dirac current enter the diagonal matrix 
entries of % n and thus the EL equations (|4.5p . Since the gauge currents have the same 
behavior under gauge transformations (see (15. 3p ). their off-diagonal elements enter the 
EL equations in the same way. For the mass terms, there is the complication that they 
have a different behavior under gauge transformations (for the logarithmic terms, this 
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was studied in (|5.16p . whereas for the contributions of the second order perturbation 
calculation, the dependence on the gauge phases can be read off from the formulas 
given in [H Def. 7.2.1]). This different behavior under gauge transformation does not 
cause problems for the logarithmic poles, because we saw in £|5,ll that the logarithmic 
poles on the light cone can be arranged to vanish. Thus the only effect of the different 
gauge behavior of the mass terms is that it modifies the values of the bosonic mass 
corresponding to the off-diagonal gauge potentials. The easiest method to describe this 
effect quantitatively is to again work with the EL equations (|4.1ip . but to modify the 
off-diagonal matrix elements of %l and % R by multiplying the contributions (|4.25[) 
(I4.3ip with numerical factors which take into account the linear behavior under off- 
diagonal left-handed gauge transformations. It is planned to work out the masses and 
coupling constants for a specific example of an admissible regularization in a separate 
publication. 



8. The Einstein Equations 

We come to the computation of the coupling constant in the Einstein-Hilbert ac- 
tion (|6.2ip . We begin by computing the mapping <J>y as defined by (I6.17p . Similar 
to (|7.3p , a short calculation using Lemma 14.111 gives 

$T(T L ,f R ) = ((f fl , f L ) - Tr c2 (f L + f R ) (1, 1)) K 8 + (deg < 4) . 

This mapping vanishes on 3 curv , so that condition (d) in §6.21 is violated. The way 
around this problem is to compute the next higher order in T reg : 

Lemma 8.1. The matrix &t defined by (|6.17p has the representation 

$t(T l ,T r ) = ((T R ,f L ) - ~ Tr c2 (f L + f R ) (1, 1)) Kg (l + 0((me)^ 

+ r re ,Tr c2 (f i + f /? )((- 1 ^.(jj _° 1 ))^n 
+ 0((me) 2pres ) (deg = 4) + (deg < 4) , 
where K\\ is the following simple fraction of degree four, 

1 1 



K 



ii 



16tt t (0) 



(-1) (o) T (-i) 
1 [R,0] 1 [0] J [0] "^ LX - 



U) 



Proof. Follows by a straightforward calculation similar as explained in the proof of 
Lemma 14.41 □ 



To explain the structure of the above result, we point out that the matrix $^ is trace- 
free. This can be understood immediately from the fact that (]4.3p vanishes if the 
terms A|A racs | are all equal. 

Combining Lemma 14.121 with Lemma 14.21 and comparing with (]6.17p , one sees that 
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A straightforward calculation shows that $t\3 cutv as computed in Lemma 18. II has the 
properties (d) and (e) in §6.21 Moreover, a short calculation yields that the condi- 
tion (f) is also satisfied, and that 

H T o$ curv (l) = ^-(l,l). 

Taking the fermion-boson pairing with (1,1) and comparing with (|6.2ip . we obtain 
the following result. 

Theorem 8.2. A ssuttic that the paTQineters S and Preg 

satisfy the scaling (14.870 . Then 

effective action involves the Einstein- Hilbert part 

£ EH = - (R + 2A) + R ■ (m 2 (me)- pres ) , (8.2) 

where R is scalar curvature, k is the gravitational constant and A the cosmological 
constant. The gravitational constant is given by 

-^16 

where K\\ and Kiq are the simple fractions (|8.ip and (|4.9Q[) . both evaluated weakly on 
the light cone (|2.32p . 

We point out that the cosmological constant remains undetermined; see the explana- 
tion after (I6.2ip . But our result excludes corrections to the Einstein-Hilbert action of 
higher order in the curvature tensor. Note that the simple fractions K\\ and K\§ are 
both of degree four, and thus their quotient is of the order one. Hence 

k ~ 5 2 . 

This means that the Planck length is to be identified with the length scale 5 describing 
the shear and general surface states (see (I2.34p and (|2.39ft ). 

We next explain how this theorem could be extended to the case 

S~ — ( m£ )¥ , (8.3) 
m 

In this case, the terms ~ m 

2 RjkC j C k in Lemma EG] are of the same order as those ~ 
Tieg/^ 2 RjkC J £ k an d must be taken into account. They can be obtained by a straightfor- 
ward computation. The statement of Theorem 18.21 will remain the same, except that 
the form of K\q will of course be modified, and that the error term in (|8.2p disappears. 
The only structural difference is that (I4.90P will then involve factors 2™ , which have 
logarithmic poles on the light cone. It does not seem possible to compensate these log- 
arithmic poles by a microlocal transformation. Therefore, in order for the logarithmic 
poles to drop out of the EL equations, one must impose that 

3 3 

a=l a=l 

This constraint for the neutrino masses can be understood similar as in Remark 14.91 
Working out the detailed computations seems an interesting project for the future. 

We finally note that for completeness, one should also compute how the energy- 
momentum tensor of the gauge fields enters the fermionic projector and verify that its 
effect on the EL equations to degree four is compatible with its coupling to the Einstein 
equations as given by the effective action. Since these computations are rather lengthy, 
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I \lfJ ~ I lb I U\ fb — lib 11 lib ^ U 

p£ W = n„,„2, - ~ (I^l + W (A.2) 



we postpone them to a future publication. These computations would conclude the 
complete analysis of the EL equations to degree four on the light cone up to errors of 
the order 

Q(x,y) = (deg = 4) • o(|£T 2 ) + (deg < 4) . 

Appendix A. The Regularized Causal Perturbation Theory with 

Neutrinos 

A.l. The General Setting. For clarity, we begin with a single Dirac sea (i.e. with 
one direct summand of (|2.42p or (|2.55p ). Thus without regularization, the vacuum is 
described as the product of the Fourier integral (jl.7p with a chiral asymmetry matrix, 

P = Xt with t = P m and X = 1, \L or xr , (A.l) 

under the constraint that X = 1 if m > 0. We again denote the regularization by 
an index e. We always assume that the regularization is homogeneous, so that P e is 
a multiplication operator in momentum space, which sometimes we denote for clarity 
by .P . If m > 0, we assume that the regularization satisfies all the conditions in [U 
Chapter 4]; see also the compilation in (TJ Section 3]. In the case m = 0, we relax the 
conditions on the shear and allow for general surface states, as explained in Section [2T3l 
In the low-energy regime, P £ should still be of the form (jA.ip . i.e. 

(# + m) 5(k 2 - m 2 ) if m > 

Xtf5{k 2 ) ifm = 

However, in the high-energy regime, P e will no longer satisfy the Dirac equation. 
But in preparation of the perturbation expansion, we need to associate the states P e 
to eigenstates of the Dirac operator (not necessarily to the eigenvalue m). To this 
end, we introduce two operators V^hift and T4hear with the following properties. The 
operator V^hift has the purpose of changing the momentum of states such that general 
surface states (as in Figure [Q (B)) are mapped onto the mass cone, i.e. 

:= ^(«8hift(*0) , (A.3) 

where v shift : M — > M is a diffeomorphism. The operator Vshean on the other hand, 
is a unitary multiplication operator in momentum space, which has the purpose of 
introducing the shear of the surface states (i.e. it should map the states in Figure H](B) 
to those in Figured (C)), 

(VshearV') (k) = V sheai (k) ip(k) with V^hear^) unitary. 

These operators are to be chosen such that the operator P £ defined by 

P £ = W Khift P £ Khift Khcar (A.4) 

is of the following form, 

pe {k)= \d{k)(^ + m{k)l)5{k 2 -m(k) 2 ) if m > 

| d{k)Xk'5{k 2 ) ifm = 

with X as in (TAT]) . Thus in the massive case, P £ should be composed of Dirac 
eigenstates corresponding to an energy-dependent mass m(k) > 0, and it should have 
vector-scalar structure. In the massless case, we demand that k 2 = 0, so that the states 
of P e are all neutral. The ansatz (|A.5P is partly a matter of convenience, and partly 
a requirement needed for the perturbation expansion (see Proposition I A. II below). 
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Moreover, we assume for convenience that P e is only composed of states of negative 
energy, 

P £ (k) = if k 2 < or k° > 0. (A.6) 

In view of (|A.2|) . it is easiest to assume that V^hift and V^hear are the identity in the 
low-energy regime, i.e. 

V r shear(^) = 1 and v shiit {k) = k if |fc°| + \k\ < . (A.7) 

Then, by comparing ()A.2p with (|A.5p . one finds that 

d(k) = l and m(k) = m if |Jfe°| + \k\ < e' 1 . (A.8) 

The required regularization of P £ (x,y) on the scale e is implemented by demanding 
that 

d(k) decays on the scale + \k\ ~ e _1 . (A. 9) 

In view of their behavior in the low-energy regime, it is natural to assume that the 
functions in ()A.7p and ()A.8p should be smooth in momentum space and that their 
derivatives scale in powers of the regularization length, i.e. 

\Vld(k)\ ~ e^\d(k)\ , |V^m(£;)| ~ e^\m(k)\ 

|V^ shcar (A;)| ~ e^\V sheai (k)\ , |V> shift (A0| ~ E W K hift (A;)| . (A ' 10) 

Clearly, the above conditions do not uniquely determine the function d and the oper- 
ators Kshift and V^hear- But we shall see that the results of our analysis will be inde- 
pendent of the choice of these operators. We remark that the transformation Kjhear is 
analogous to the transformations XJ\ considered in [H Appendix D] (see [H eq. (D.22)]), 
except that here we consider only one unitary transformation. 

The last construction immediately generalizes to a system of Dirac seas. Namely, 
suppose that without regularization, the auxiliary fermionic projector of the vacuum 
is a direct sum of Dirac seas (see for example (I2.42|) or fl2.55f) ) , 

^max 

paux = x ^ 

1=1 

Then we introduce P £ simply by taking the direct sum of the corresponding regularized 
seas 

^max ^max ^max 

paux ._ ^ pe ^ V shi{t := ^ VfOffi > Khoar : = @ Kshear ■ 

i=l l=\ 1=1 

Setting 

paux i/ t/ pe \t— 1 \t— 1 

~~ ''shear ''shift r " s hift ''shear ' 

the operator P £ satisfies the Dirac equation 

{tf -mY(k))P £ {k) =0, (A.ll) 
where the mass matrix is given by (cf. (|2.44j) or (I2.57P ) 

(max 

mY(k) = rri£ . 
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In the low-energy regime, we know furthermore that 

Vshearik) = 1 and v shi{t (k) = k ) 

} if |/c u | + \k\ < e , 
P aux (/c) = Xt J 

where X and t are given as in (|2.45|) . Clearly, the regularity assumptions (jA.lOp are 
imposed similarly to p aux . Finally, we need to specify what we mean by saying that 
two Dirac seas are regularized in the same way. The difficulty is that, as mentioned 
above, different choices of d, V s h eaiT and t^wft ma Y gi ve r i se to the same regularization 
effects. In order to keep the situation reasonably simple, we use the convention that if 
we want two Dirac seas to show the same regularization effects, we choose the corre- 
sponding functions d as well as V^hear and t^hift to be exactly the same. If conversely 
two Dirac seas should show different regularization effects, we already choose the cor- 
responding functions d to be different. Then we can say that two Dirac seas labeled 
by a and b are regularized in the same way if d a = df,. In this case, our convention 
is that also (T^heax)a = (Kjhear)& and (V r s hift)a = (V s H{t)b- This notion gives rise to an 
equivalence relation on the Dirac seas. In the formalism of Section f2.6|. the equivalence 
classes will be labeled by the parameters r 4 reg (see ()2.43p and (|2.58p ). 

A. 2. Formal Introduction of the Interaction. Now the interaction can be intro- 
duced most conveniently by using the unitary perturbation flow \10\ Section 5]. In 
order not to get confused with the mass matrix, we introduce an additional spectral 
parameter /x into the free Dirac equation, which in momentum space reads 

(Jf-mY(k) -//!)#(£) = 0. 

For this Dirac equation, we can introduce the spectral projectors p, the causal funda- 
mental solutions k and the symmetric Green's functions s can be introduced just as 
in [H §2.2], if only in the formulas in momentum space we replace m by mY{k). For 
clarity, we denote the dependence on fi by an subscript +fj, (this notation was used 
similarly in [U §2.6]; see also [H §C.3] for an additional "modified mass scaling", which 
we will for simplicity not consider here). We describe the interaction by inserting an 
operator 2 into the Dirac operator, 

V = i$ + 2 - mY(k) . 

After adding the subscript +/u to all factors p, k or s in the operator products in [TO] 
Section 5], we obtain an operator U which associates to every solution ip of the free 
Dirac equation (ity — mY — = a corresponding solution rp of the interacting 

Dirac equation (ity + 2 — mY — \i l)ip = 0, 

17(8) : ip^ip. 

The operator U is uniquely defined in terms of a formal power series in 2. Taking fi 
as a free parameter, in (TOl Section 5] the operator U is shown to be unitary with 
respect to the indefinite inner product (|2.3p . We now use U to unitarily transform all 
the Dirac states contained in the operator P £ and set 



^ aux = W ^shift C/(8) P s UCB)- 1 V s ^ t . (A.12) 



This construction uniquely defines the regularized auxiliary fermionic projector with 
interaction P aux in terms of a formal power expansion in 2. The fermionic projector 
is then obtained by taking the partial trace (see (|2.53p or (|2.56p ). 
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A. 3. Compatibility Conditions for the Interaction. In order to derive the struc- 
ture of the admissible 23, we first consider a perturbation calculation to first order and 
assume that 23 is a multiplication operator in position space having the form of a plane 
wave of momentum g, 

B(x) = r B q e~ iqx . (A.13) 

In this case (cf. @J eq. (D.14)]), 

/oo 
dp{s + ^p + ^P e + P E p + ^s +fl ) . 
-oo 

Using a matrix notation in the direct sums with indices a, b € {1, ... , ^ m ax}, we obtain 
in momentum space (for the notation see [U Chapter 2] or |10j ) 

/oo 
d(i \s M (k + q) (2,)g Pmb+M (£0 CP £ )£(fc) 
-°° (A.14) 
+ (ps)»(fc + g)p ma+M (fc + g) (B,)g 



This equation was already considered in [2, Section 3] and [U Appendix D]. However, 
here we analyze the situation more systematically and in a more general context, 
pointing out the partial results which were obtained previously. 

For clarity, we analyze (|A.14p step by step, beginning with the diagonal elements. 
For ease in notation, we assume that "B q has only one non-trivial component, which is 
on the diagonal, 

(H q )t = 6 at S M B (A. 15) 

with £ G {1, . . . , £ max } and B a matrix acting on Dirac spinors. Shifting the integration 
variable according to mi + p — > fi, we obtain (cf. [H eq. (D.15)]) 

(AP £ Y e (k + q,k) 

dpi [s^k + q)B p^k) (P e )[{k) + (P £ )l(k + q) Pll (k + q)B s^k)} 



— oo 
oo 



pp 

(k + p 2 

PP 



d » ^ \ jr. , ^ .2 (# + i + y) B W + A*) Kk 2 - fJ 2 ) (P £ ) £ t(k) 



+ (P%(k + g) 5((k + g) 2 — (tf + 4 + p) B (Ifi + fj) 



f 

j — t 



k 2 — fi 2 
PP r 

2A; t 2 {(# + + m) b + m) 5(^ 2 - z^ 2 ) (P £ )I(A;) 



2% + g 2 

- (P £ )|(A; + g) 5((k + g) 2 - M 2 ) (0 + ^ + M ) B (0 + //)}, 

where in the last step we used that the argument of the 5-distribution vanishes. Car- 
rying out the //-integration gives (cf. [U eq. (D.15)]) 

{AP%(k + q,k) = - ^ 2 { ((# + rf) B + B k) (P%{k) 

4kq + Zq v\ J (A. 16) 



Here the principal part has poles if 2kq + q 2 = 0, leading to a potential divergence 
of AP 6 . In order to explain the nature of this divergence, we first point out that if 23 
had been chosen to be a smooth function with rapid decay, then AP would have been 
finite (see the proof of [U Lemma 2.2.2]). Thus the potential divergence is related to the 
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fact that the plane wave in (|A.13P does not decay at infinity. A more detailed picture is 
obtained by performing the light-cone expansion (see [2] and [H Appendix F]). Then 
one can introduce the notion that (|A.16|) is causal if its light-cone expansion only 
involves integrals along a line segment xy. Since such integrals are uniformly bounded, 
it follows immediately that all contributions to the light-cone expansion are finite for 
all q. If conversely (|A.16P diverges, then the analysis in [4, Appendix F] reveals that 
individual contributions to the light-cone expansion do diverge, so that unbounded 
line integrals must appear (see also the explicit light-cone expansions in pQ). In this 
way, one gets a connection between the boundedness of (|A.16P and the causality of the 
light- cone expansion. 

Unbounded line integrals lead to contributions to the EL equations whose scaling 
behavior in the radius is different from all other contributions. Therefore, the EL 
equations are satisfied only if all unbounded line integrals drop out. The easiest way 
to arrange this is to demand that the fermionic projector itself should not involve any 
unbounded line integrals. This is our motivation for imposing that 

(AP) £ (k + q, k) should be bounded locally uniformly in q. (A. 17) 

Let us analyze this boundedness condition for (|A.16j) . Since the denominator 
in (IA.16D vanishes at q — > 0, we clearly get the necessary condition that the curly 
brackets must vanish at q = 0, 

[{Jf,B},P £ (k)} =0. (A.18) 
Using (1A.5I) together with the identity 

[{#, b} , = [k 2 , b] + m - W = , 

we find that (|A.18h is automatically satisfied in the case X = 1. The situation is 
more interesting if a chiral asymmetry is present. If for example X = xl, we get the 
condition 

[{l j t,B},XLk']=0. 

This condition is again trivial if B is odd (meaning that {£>,7 5 } = 0). However, if B 
is even, we conclude that 

0=y{{#,£},W=7 5 W 

(where in the last step we used that k 2 = in view of (|A.5j) ). As k is any state on the 
lower mass shell, this rules out that B is a bilinear potential, leaving us with a scalar 
or a pseudoscalar potential. In order to rule out these potentials, we next choose a 
vector q with qk = 0, set q = eq and consider (|A.16p in the limit e — > 0. Then the 
denominator in (lA.16j) diverges like e~ 2 , so that the curly brackets must tend to zero 
even ~ e 2 . 

((# + ei)B + Bk^jP £ (k) - P £ {k + £<?)((# + e|)B + = 0(e 2 ) . (A.19) 

Using that P £ (k) is left-handed and that B is even, we find that the first summand 
in ()A.19p is right-handed, whereas the second summand is left-handed. Hence both 
summand must vanish separately, and thus 

= + |)B + Bfy) P £ {k) = q'Bk' d{k) 5{k 2 ) . 
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This condition implies that B must vanish. We conclude that if X = xl, only odd 
potentials may occur. We can write this result more generally as 



■BX = X*% . 



(A.20) 

We have thus derived the causality compatibility condition (|2.49p from our boundedness 
condition (|A.17[) . This derivation is an alternative to the method in [4j §2.3], where 
the same condition was introduced by the requirement that it should be possible to 
commute the chiral asymmetry matrix through the perturbation expansion. 

So far, we considered (|A.17|) in the limit q — > 0. We now analyze this condition for 
general q. Using (|A.5|) and (|A.20|) . a short calculation gives 

Akq + 2q 2 

we we set m = rri£. If d(k) = d(k + q), the transformations 

f 5' (k 2 - m 2 + r(2kq + q 2 )) dr = — —5(k 2 -m 2 +T(2kq + q 2 ))dT 

Jo 2kq + q z J dr 

1 



(AP £ )f (fc + q,k) = -X (d(k) 5(k 2 - m 2 ) - d(k + q) 5({k + q) 2 - m 2 )) . 



5((k + q) 2 - m 2 ) - 5(k 2 - m 2 )) 



X { lkafi 2 6{k2 ~ ^ ( d{k + q) ~ d(k) ) ■ (A ' 21) 



2kq + q 2 

show that AP £ is indeed a bounded distribution for any q. Thus it remains to be 
concerned about the contribution if 5(k) / 5{k + q), 

4kq + 2q 2 

Unless in the trivial case B = 0, this contribution is infinite at the poles of the de- 
nominator. We conclude that in order to comply with the condition (|A.17|) . we must 
impose that the weight function d(k) in (|A.5j) is constant on the mass shell k 2 = m(k) 2 . 
This is indeed the case in the low-energy regime (|A.8p . However, in the high-energy 
region, the function d{k) is in general not a constant (and indeed, assuming that d(k) 
is constant would be in contradiction to (|A.9|) ). Our way out of this problem is 
to observe that (|A.2ip implies that the light-cone expansion of (|A.16P is in general 
not causal, in the sense that it involves unbounded line integrals. However, using 
that d(k + q) — d(k) ~ <?|Vd|, the scalings q ~ ^ m i C ro an d (jA.lOj) show that these 
non-causal contributions to the light-cone expansion are of 

higher order in e/^macro • 

(A.22) 

This consideration shows that the perturbation expansion will give rise to error terms 
of higher order in ^/-^macro* In what follows, we will always neglect such error terms. 
If this is done, the above assumptions are consistent and in agreement with (1A.17|1 . 
provided that the causality compatibility condition (|A.2Q[) holds. 

Before moving on to potentials which mix different Dirac seas, we remark that the 
above arguments can also be used to derive constraints for the possible form of P £ (k), 
thus partly justifying our ansatz lAJjl 

Proposition A.l (Possible form of P £ ). Suppose that AP £ as given by (fA~T6i) satisfi es 
the condition (|A.17p . We renounce the assumptions on P £ (see (|A.5|) . (|A.6[) . (|A.8P 
and (jA.lOp ), but we assume instead that the admissible interaction includes chiral or 
axial potentials. Then for every k, there are complex coefficients a,b,c,d£ C such that 

P £ (k) = al + ibj 5 + ck' + dj 5 k' . (A.23) 



THE CONTINUUM LIMIT OF A FERMION SYSTEM INVOLVING NEUTRINOS 



81 



Proof. We again analyze the necessary condition (|A.18|) . This condition is clearly 
satisfied if B is a vector potential. Thus by linearity, we may assume that B = >y 5 4 is 
a vector potential. It follows that 

o=[{#, 7 5 4},p £ (fc)] = [ 7 5 [4^},P%k)). 

Decomposing P £ (k) into its even and odd components, by linearity we can again 
consider these components after each other. If P £ (k) is odd, i.e. 

we obtain 

= 7 5 {[4,^], J P £ (A;)} =Aj 5 ((ku)4- {Au)^) + 4((kv) 4 - (Av) If) . 

Since A is arbitrary, it follows that u and v must be multiples of k. 
If P £ (k) is even, we obtain the condition 

[\m,p £ (k)\ =o. 

This condition is obviously satisfied if P £ (k) is a scalar or a pseudoscalar. Writing 
the bilinear component of P £ (k) in the form F^rf^ with an anti-symmetric tensor 
field F, we get the condition 

= [[4, , F ii7 V] = 2F lJ k i [4, y] - 2F ij A i 7 J ] • 

Since ^4 is arbitrary, it follows that F = 0, concluding the proof. □ 

The step from ()A.23P to our the stronger assumption (|A,2|) could be justified by the 
assumption that the image of P £ should be negative definite or neutral, and further- 
more by assuming that without chiral asymmetry the matrix 7 s is absent, whereas the 
chiral asymmetry is then introduced simply by multiplying with \L or XR- We finally 
remark that in [U Appendix D], the condition (|A.18|) is analyzed for B a scalar poten- 
tial to conclude that P £ (k) should commute with the Dirac operator (see [U eq. (D.16) 
and eq. (D.17)]). This is consistent with our ansatz (|A.5j) . which is even a solution of 
the Dirac equation (jA.llj) . However, we here preferred to avoid working with scalar 
potentials, which do not seem crucial for physically realistic models. 

We next consider instead of (|A.15j) a general potential "B q which may have off- 
diagonal terms in the direct summands. Then ()A.14j) can be evaluated similar as in 
the computation after (IA.15j) . but the calculation is a bit more complicated. Therefore, 
we first compute the integral of the first summand in (jA. 14j) , 

dfi s ma+fl (k + q) (S,)g Pmb+IM (k) (P £ ) b b (k) 



d[i e(m b + fj,) lu ; - , 2 S(k 2 - (m b + fi) 2 ) 



PP 

(k + q) 2 - (m a + 
x(k' + 4 + m a + ti) (S,)g + m b + n) (P £ ) b b (k) 



— — (P E )$(k) , (A.24) 

[i=±\k\—mt > 



where we set 



fga _ 1 pp ( W + 4 + m a + ii) (gg)g (lf + m b + ri { 
b 2 \2kq + q 2 - (m 2 a - m 2 ) - 2^(m a - m b ) } 
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and \k\ = \/W (note that, in view of our assumption (|A.6p . the factor (P e )^(k) guar- 
antees that the above expression vanishes if k 2 < 0). Treating the second summand 
in ()A,14p similarly, we obtain 

tra E lJ %$^*i- £ ® ? ^t? ■ (A ' 26) 

/j,=±\k+q\—m a jtt=±|fe|— mi 

This formula is rather involved, but fortunately we do not need to enter a detailed 
analysis. It suffices to observe that (|A.25|) has poles in q, which lead to singularities 
of (IA,24p . Thus the only way to satisfy the condition (|A.17j) is to arrange that contri- 
butions of the first and second expression on the right of (|A.26j) cancel each other. In 
view of CO]), the first expression involves d a (k + q), whereas in the second expression 
the term db(k) appears. This shows that in order to get the required cancellations, 
the functions d a and <4 must coincide. Using the notion introduced on page [771 we 
conclude that ¥> may describe an interaction of Dirac seas only if they are regularized 
in the same way. An interaction of Dirac seas with different regularization, however, 
is prohibited by the causality condition for the light-cone expansion. For brevity, we 
also say that the interaction must be regularity compatible. 

A. 4. The Causal Perturbation Expansion with Regularization. We are now 

ready to perform the causal perturbation expansion. In |10|. Section 5], the unitary 
perturbation flow is introduced in terms of an operator product expansion. Replacing 
the Green's functions and fundamental solutions in this expansion by the corresponding 
operators of the free Dirac equation (|A.11|) . we can write the operator U("B) P £ 
as a series of operator products of the form 

Z, '.= C\ H> • • • S Cp H P e H> Cp+i S • • • S Cfc , 

where the factors G\ are the Green's functions or fundamental solutions of the free 
Dirac equation (|A.ll|h The operators Ci are diagonal in momentum space, whereas 
the potential "B varies on the macroscopic scale and thus changes the momentum only 
on the scale ^ macro . Thus all the factors C\ will be evaluated at the same momentum p, 
up to errors of the order ^ macro . We refer to this momentum p, determined only up 
to summands of the order ^ macro as the considered momentum scale. In view of the 
regularity assumptions on the functions d and m in (jA.lOj) . we may replace them by 
the constants d(p) and m(p), making an error of the order (1A.22|) . This evaluation 
of the regularization functions is referred to as the fixing of the momentum scale, and 
we indicate it symbolically by | S caiep- Since "B is regularity compatible, we may then 
commute the constant matrix d to the left. Moreover, we can apply the causality 
compatibility condition (|A.20|) together with the form of X in (|A.5[) and (|A.2[) to also 
commute the chiral asymmetry matrix X to the left. We thus obtain the expansion 

ump £ um~\ calep = Y: E ^xdc^c,,^ •••£c fe+M L Jep (A.27) 



k=0 a=0 



+ (higher orders in £/ 



-macro; 7 



where we set 

^max ^max 

I = $If and d=0d € (p), 
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and c a are combinatorial factors. Here the combinatorics of the operator products 
coincides precisely with that of the causal perturbation expansion for the fermionic 
projector in |10t Theorem 4.1]. 

A. 5. The Behavior under Gauge Transformations. In order to analyze the be- 
havior of the above expansion under [/(l)-gauge transformations, we consider the case 
of a pure gauge potential, i.e. B = with a real-valued function A. Then the gauge 
invariance of the causal perturbation expansion yields 

17(B) P £ ^(B)- 1 !^ p = (e iA P £ e- iA )\ SC3le p + (higher orders in e/Wo) • (A.28) 

According to (|A.12[) . we obtain p aux by applying with Vghift and V^hear- The transfor- 
mation Kshift is a subtle point which requires a detailed explanation. We first consider 
its action on a multiplication operator in momentum space M{k). Then, according to 
the definition (|A.3|) . 

(V^tMV^ik) = (MV s ^){v ahm (k)) 

(ushift(fe)) ip(k) (A. 29) 

so that the transformation again yields a multiplication operator, but with a trans- 
formed argument. To derive the transformation law for multiplication operators in 
position space, we first let / = e~ iqx be the operator of multiplication by a plane 
wave. Then 

(V shift fV s ^)(k) = (fV s ^){v shiit (k)) 

= (KhiftV^^W^) -q)= ^~J ft (vsmft(fc) - ?))• 

This can be simplified further if we assume that the momentum q ~ ^ m a C ro ^ s macro- 
scopic. Namely, the scaling of the function u s hift in (|A.10|) allows us to expand in a 
Taylor series in q, 

(y shi{t fV sh l h ij)(k) =tp(k- Dv^l^^ q^j + (higher orders in e/l macro ) . 

Thus V^hift / Khift * s again a multiplication operator in position space, but now corre- 
sponding to the new momentum 

L(k)q with L(k) := Dv;^ it \ Vshift{k) . 

Again in view of the regularity assumptions (jA.lOp , when fixing the momentum scale 
we may replace the argument k by p, i.e. 

Khift Khiftlscale P = e~^ p ^ x + (higher orders in e/Wo) • 

Using the relation (L(p)q) x = qL(p)*x, we can rewrite this transformation law simply 
as a linear transformation of the space-time coordinates. Then the transformation 
law generalizes by linearity to a general multiplication operator by a function / which 
varies on the macroscopic scale, i.e. 

Khift f(x) KiMtLaOe P = f( L (p)* x ) + (higher orders in e/£ macIO ) . (A.30) 

With (IA.29P and (jA.30j) . we can transform (lA.28j) to the required form (|A.12j) . Using 
that V^hear commutes with scalar and macroscopic multiplication operators (again up 
to higher orders in e/£ ma _ cro ), we obtain in view of (IA.4j) 

p£ (*>y)Laie P = e lA(L(p) * x) P £ (x,y) e-^ L ^ + (higher orders in e/Wo) . (A.31) 
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Except for the factors L(p)*, this formula describes the usual behavior of the fermionic 
projector under gauge transformations. In particular, if we do not consider general 
surface states and V^hift = lj then our perturbation expansion is gauge invariant. How- 
ever, if we consider general surface states described by a non-trivial operator V^hift) then 
the matrix L will in in general not be the identity, and the transformation law (|A.3ip 
violates gauge invariance. 

Our method to recover gauge invariance is to replace the gauge potential A by a 
more general operator A, which in momentum space has the form 

A(v+ q -,p- q -) :=A{L{pT\), (A.32) 

where A is the Fourier transform of the classical potential to be used when no regular- 
ization is present. In the case A = qDA of a pure gauge field and fixing the momentum 
scale, we then find that A coincides with the multiplication operator -A(((L(p)*) _1 a:)), 
just compensating the factors L(p)* in (|A.31j) . In view of the regularity assump- 
tions (lA.lOj) . the matrix L scales in powers of the regularization length. Thus A is 
a nonlocal operator, but only on the microscopic scale e. On the macroscopic scale, 
however, it coincides with the classical local potential. We also point out that the 
compatibility conditions worked out in Section IA.3I under the assumption that 23 is a 
multiplication operator are valid just as well for the nonlocal potential (1A.32|) . because 
after fixing the momentum scale, A reduces to a multiplication operator, so that our 
previous considerations again apply. 

A. 6. The Regularized Light-Cone Expansion. We are now in the position to 
perform the light-cone expansion. Our starting point is the operator product expan- 
sion (|A.27j) . We choose 23 to be the Fourier transform of a general multiplication or 
differential operator and introduce 13 in analogy to (|A.32p as a non-local operator. 
After fixing the momentum scale, this operator reduces again to a multiplication or 
differential operator. Then the light-cone expansion can be performed exactly as de- 
scribed in [2] and [H §2.5]. Finally, one can transform the obtained formulas with V^hift 
and Kjhear (again using the rules (|A.29P and (|A.30|) ). Since the resulting line integrals 
do not depend on the momentum scale, the regularization only affects the factors T^ n \ 
The condition that 23 should be regularity compatible can be described by the param- 
eters T 4 reg . We thus obtain the formalism described in §2.51 and §2.61 

We finally compare our constructions with those in [2J Appendix D]. Clearly, the 
constructions here are much more general because they apply to any order in pertur- 
bation theory and may involve a chiral asymmetry. Moreover, the momentum shift 
operator V^hlft makes it possible to describe general surface states, and also we allow 
for a large shear of the surface states (whereas in [U Appendix D] we always assumed 
the shear to be small). Nevertheless, the basic idea that in order to preserve the 
gauge invariance in the presence of a regularization, one should replace the classical 
potentials by operators which are nonlocal on the microscopic scale already appeared 
in [H Appendix D] (see the explanation after [U eq. (D.26)]). Thus [4, Appendix D] 
can be regarded as a technical and conceptual preparation which is superseded by the 
constructions given here. 
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Notation Index 



P(x, y) - kernel of fermionic projector, 

m 

A xy - corresponding closed chain, [2] 

/ZfAjy] - causal Lagrangian, [2] 

| . | - spectral weight, [2] 

S[P] - causal action, [2] 

T[P] - constraint, [2] 

(M, (.,.)) - Minkowski space, E] 

- auxiliary action, [2] 
fi - Lagrange multiplier, E] 
P N (x,y) - vacuum fermionic projector 

in neutrino sector, [3] 
P (x,y) - vacuum fermionic projector 

in charged sector, [3] 
mp - masses of charged fermions, E] 
P m - vacuum Dirac sea of mass to, [3] 
mp - neutrino masses, [3] 
mY - mass matrix, [19] 
C^mns ~ Maki-Nakagawa-Sakata (MNS) 

matrix, H [28] 
M a - bosonic mass, H 
c a - coupling constant, H 
Rjk - Ricci tensor, H [52] 
R - scalar curvature, H EH [74] 
Tjk - energy- momentum tensor, d [51] 
A - cosmological constant, [4] [69] [74] 
k - gravitational constant, [1 [71 
pe _ re g U l ar i zec l fermionic projector, [5] 
£ - vector y — x, E] 

Q{x,y) - operator in EL equations, [5] 
El 

<.|.> - inner product on wave func- 
tions, [6] 

Aj - eigenvalues of the closed chain, [10] 
T m 2(x,y) - Fourier transform of lower 
mass shell, [T5l 
- regularized term of mass expan- 
sion, [T6l 
T|™| - ordinary shear term, [TBI 
c reg - regularization parameter, [16] 
deg - degree on light cone, [16] 
5 - length scale of shear and general 
surface states, [T7] 



T\ R p i - describes mass expansion of gen- 
eral surface states, [T7I 
Tj"' p j - describes shear states, [J7J 

Zo™' ) - factor for contraction rules, [T7] 
P aux - auxiliary fermionic projector, [18] 
X - chiral asymmetry matrix, [18] 
r reg ~~ dimensionless parameter for high- 
energy states, CEl [30] 
t - distribution composed of vacuum 

Dirac seas, [19] 
t - corresponding object with interac- 
tion, HH] 

23 - operator in Dirac equation, [191 EHl 
M 

v , A - denote the partial trace, [20] 
£ - real lightlike vector in t-formalism, 

EU 

P^ - lightlike component of vacuum fer- 
mionic projector in t-formalism, 

EU 

- vector describing regularization in 
i- formalism, [22] 
AnL ~ eigenvalues of closed chain, [241 
Fncs ~ corresponding spectral projectors, 

El 

T {n) _ T (n) T (n) ,„ 

?7 C - unitary matrix involving gauge pha- 
ses^ 

Pexp, Pe - ordered exponential, E7J [55] 



A L - 


- left-handed gauge potential, [29] 
[56] 




right-handed potential in neutrino 




sector, E9J EH [56] 




right-handed potential in charged 




sector, EH [51 [56] 


Preg " 


- determines scaling r reg = (me) Pies , 


[30] 


^nc — 


eigenvalues of matrix involving pha- 




ses, [31] 




corresponding spectral projectors, 




EU 


A±- 


eigenvalues of closed chain in vac- 




uum, EU 



m 



F. FINSTER 



F± - corresponding spectral projectors, 

ED 

AQ(x,y) - operator Q to degree four 
on light cone, [33] 

% nc - matrices entering the EL equa- 
tions to degree four, [33] 

% c - matrices entering the EL equations 
to degree four, [33] 

1Z C - matrices entering the EL equa- 
tions to degree four, [35] 

°(l£| fc ) _ order at the origin. |3"T1 

j c - bosonic current. [371 

J c - Dirac current, [37] 

2c ~ matrix composed of current and 
mass terms, [37] 

U(k) - homogeneous transformation . HU1 

Q - absolute value of energy, 00] 

Z(k) - generator of homogeneous trans- 
formation, 0o] 

L(k), R(k) - chiral components of Z(k), 

M 

Sq, S*2 - signature matrices, 03] 
U(x, y) - microlocal chiral transforma- 
tion, 06] 

D even , D dd ~~ even and odd components 

of Dirac operator, 08] 
feven _ even component with flipped 

chirality, 1351 
ei, . . . , C6 - orthonormal basis, l49l 
MjP - short notation for factors 

or Lf Q ],W\ 
Symm(C n ) - Hermitian n x re-matrices, 

m\ 

& n = Symm(C 6 ) Symm(C 6 ) - left- 
and right-handed matrices, [65] 
- dynamical subspace, [66] 
^Dirac _ Dirac Lagrangian, [66] 
f - matrix taking the partial trace, [66] 
(.|.) - fermion-boson pairing. [671 
3, 9t - interacting and non-interacting 
subspace, [67] [71] 

_ describe coupling to EL 
equations, [67] 
Hj - inverts $ j on 3, [67] 

^curv _ describe coupling to EL equa- 
tions, [68] 



3 CU rv> 9t C urv ~ interacting and non- inter- 
acting subspace for gravity, [68] 

£ym ~ Yang-Mills Lagrangian, [69] 

S e g - effective action, [69] 

5eh _ Einstein-Hilbert action, [69j [71] 

Kshift, Kihear - operators in the regular- 
ized causal perturbation theory, 

ITS 
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Subject Index 



/.-formalism, |2"T1 

action 

auxiliary, [2] 
causal, [2] 
Dirac, IBB1 
effective, ES M 
Einstein-Hilbert, ESI El 
Yang-Mills, EH 

bilinear logarithmic term, [54] 

causal perturbation theory, [T9l 
causality compatibility condition , [T9l [801 
chiral asymmetry matrix, [18] 
chiral cancellation, [5] [10] 
contraction rule, [T7J 
cosmological constant, [69j [Til 
curvature term, [52l 

degree, [TBI 

dynamical subspace of ©6 , EBJ 

Einstein equations, H 
energy-momentum term, [51] 
evaluation on the light cone 

pointwise, [30] 

weak, [TBI 

fermion-boson pairing, [BTJ 
fermionic projector 

auxiliary, [18] 

of the vacuum, [3] 
field tensor term, [B01 

gauge field equations, EJJ 
gauge group 

dynamical, IBB1 
gauge potential 

dynamical, IBB1 
general surface state, [10] 
gravitational constant. [Til 

integration-by-parts rule, [TTJ 
interacting subspace of ©2, [671 

Lagrangian 



causal, [2] 
Dirac, [BU] 
effective, ESI M 
Einstein-Hilbert, EH El 
Yang-Mills, [69] 
light-cone expansion, [T9l 

Maki-Nakagawa-Sakata (MNS) matrix, 

mass expansion, EU [161 E3 EH 
mass matrix, [19] 
massive state, [T3l 

neutrino 

chiral, El El 

massive, [3] 
non- interacting subspace of 69. [671 
null state, [TBI 

order 

at the origin, [37] 

projection on the dynamical degrees of 
freedom, [68l 

regularization 

macroscopic away from light cone, 

El EH 

naive, El E] 

spherically symmetric, O [TOl 
regularization parameter, [TBI 

sector 

charged, [3] 

neutrino, [3] 
shear contribution 

by the microlocal chiral transfor- 
mation, EH] 
shear state, E2 
simple fraction, E7J 
singular mass limit, E] 
state stability, El [Til 

transformation of the fermionic projec- 
tor 

microlocal chiral, EfJ 
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